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INTRODUCTION 


This  document  contains  reproductions  of  technical  papers  on  some 
of  the  most  recent  research  results  on  aircraft  structures  from  the 
NACA  Laboratories.  These  papers  vere  presented  by  members  of  the 
staff  of  the  NACA  Laboratories  at  the  NACA  conference  held  at  the 
Langley  Aeronautical  Laboratory  March  7>  19^1>  The  purpose  of  this 
ccxiference  vas  to  convey  to  those  Involved  In  the  study  of  aircraft 
structures  these  recent  research  results  aM  to  provide  those  attending 
an  opportunity  for  discussion  of  the  results. 

The  papers  In  this  document  are  In  the  same  form  In  \dilch  they 
were  presented  at  the  conference  so  that  distribution  of  them  ml^t 
be  pxxmpt.  The  original  presentation  and  this  record  are  considered 
as  ccsqplementary  to^  rather  than  as  substitutes  for,  the  Ca^ttee's 
system  of  cosqplete  and  formal  reports. 

A  list  of  the  conferees  Is  Included. 
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SOME  CONSIllERATIONS  OF  WING-FLEXIBILITY  EFFECTS 
ON  RESPONSE  DUE  TO  GUSTS 
By  Eldon  E.  Kordes  and  John  C.  Houbolt 
Langley  Aeronautical  Laboratory 


There  is  continued  concern  as  to  the  role  that  structural  flexi¬ 
bility  plays  in  airplane  reaction  due  to  gusts  and  in  gust  measurements. 
This  paper  is  intended  to  give  an  evaluation  of  the  ii^portance  of  wing 
bending  flexibility. 

Analyses  that  atten^t  to  deal  thoroughly  with  all  the  structural 
flexibilities  are  -rather  involved  and  are  not  well  suited  for  evaluation 
studies  which  require  the  calculation  of  a  great  many  individual  cases. 
Therefore,  as  a  basic  stepping  point,  the  problem  was  idealized  somewhat 
as  follows.  Airplane  motion  was  assumed  to  be  represented  two  degrees 
of  freedom:  vertical  motion  of  the  airplane  and  fundamental  wing 
bending.  Two-dimensional  unsteady-lift  theoiy  was  used  and  it  was 
assumed  that  the  idealized  airplane  encounters  prescribed  gust  varia¬ 
tions  in  the  direction  of  flight.  The  thought  behind  representing  air¬ 
plane  motion  ly  two  degrees  of  freedom  was  to  see  whether  a  method  suit¬ 
able  for  making  trend  studies  could  be  found  and  to  see  if  some  basic 
parameters  related  to  flexibility  effects  could  be  established. 

A  classical  treatment  of  the  idealized  system  led  to  response 
equations  vdiich  si;iggested  certain  natural  parameters  and  these  are  shown 
in  figiire  1.  The  purpose  of  this  figure  is  to  show  the  nature  of  the 
parameters  since  they  denote  conditions  for  dynamic  similarity;  airplanes 
having  similar  values  of  these  parameters  would  have  similar  responses 
(that  is,  have  similar  values  for  gust-response  factor  as  well  as  for 
certain  other  nondimensional  quantities  relating  to  the  effects  of  the 
gust  on  the  airplane).  The  first  parameter  is  associated  with 
vertical  motion  of  the  airplane.  It  is  very  similar  to  the  relative 
density  or  so-called  mass  parameter  that  has  been  referred  to  in  other 
gust  studies.  In  this  expression  W  is  the  gross  weight  of  the  air¬ 
plane,  n  is  the  value  of  one-half  the  theoretical  lift-curve  slope, 
p  is  air  density,  g  is  the  acceleration  of  gravity,  Cq  is  a  refer¬ 
ence  chord,  and  S  is  wing  area.  Because  wing  bending  is  included,  a 
similar  parameter  is  also  found  and  might  be  considered  as  the  mass 

parameter  associated  with  the  foindamental  mode.  The  only  difference 
between  this  parameter  and  is  that  W  is  replaced  by  W^,  where 

is  defined  by  the  integral  shown,  which  involves  the  span^se  weight 
distribution  w  and  the  fundamental-mode  curve  y^^  normalized  to  unit 
tip  deflection.  The  next  parameter  is  X  which  is  given  in  terms  of 
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the  natural  frequency  of  the  fundamental  mode,  the  reference  chord,  and 
the  forward  velocity  of  the  airplane.  This  parameter  is  very  similar 
to  the  reduced  frequency  parameter  used  in  flutter.  The  fourth  param^ 
eter  H  represents  the  distance  in  chords  from  the  beginning  of  an 
assumed  single  gust  to  its  peak  vertical-velocity  value.  The  next  paraan 

eter  —  is  associated  with  the  bending  moment  due  to  lift  and  with  the 

alleviation  moment  associated  with  mass  inertia  effects,  is  the 

distance  from  the  center  line  of  the  airplane  to  the  center  of  lift,  and 
x„j  is  the  distance  from  the  center  line  to  the  center  of  gravity  of  the 
semispan  mass. 

It  may  be  of  interest  to  see  how  the  two-degree-of -freedom  analysis 
which  led  to  these  parameters  agrees  with  a  more  complete  analysis 
method. 

In  figure  2  the  bending  stress  developed  near  the  root  of  a  partic¬ 
ular  twin-engine  transport  airplane  flying  into  a  10-foot-per-second 
sharp-edge  gust  is  shown.  The  results  obtained  by  a  recurrence-matrix 
method  which  involves  the  use  of  influence  coefficients  at  six  semispan 
stations  (reference  l)  and  the  two-degree-of -freedom  analysis  agree 
rather  well  over  the  time  interval  shown,  which  includes  the  maximum 
stress  developed.  These  results  lend  a  measure  of  confidence  in  the  use 
of  the  simplified  analysis  for  making  trend  studies. 

In  the  evaluation  of  flexibility  effects  a  possible  approach  would 
be  to  investigate  what  effect  individual  variations  in  the  basic  param¬ 
eters  would  have  on  airplane  response.  This  approach  has  the  objection, 
however,  that  the  physical  significance  may  not  at  once  be  evident. 
Therefore,  it  may  be  more  instructive  to  consider  what  happens  when 
operating  conditions  for  a  typical  airplane  are  varied,  such  as  a  change 
in  speed  or  a  change  in  mass  distribution,  even  though  a  change  in  one 
operating  condition  may  cause  a  change  in  several  of  the  parameters. 

Thus  the  present  evaluation  is  made  ty  considering  four  specific  oper¬ 
ating  effects  and  these  are  listed  in  figure  3.  The  first  effect  deals 
with  gust  gradient  distance.  The  second  phase  deals  with  gust  shape  and 
for  the  present  investigation  three  types  have  been  considered;  they  are 
a  sine  gust,  a  sine  squared  gust  (which  may  also  be  called  a 

1  -  cos  gust),  and  a  triangular  gust.  The  present  investigation  is 
limited  to  single  gusts,  but  it  is  recognized  that  gust  sequence  or  pene¬ 
tration  of  successive  gusts  may  be  important.  The  third  effect  involves 
the  static  loading  conditions.  The  case  of  a  twin-engine  airplane  with 
fuel  tanks  outboard  of  the  engines  is  chosen  and  three  loading  conditions 
are  considered.  The  first  is  that  with  a  full  load  in  the  fuselage  and 
wing  fuel  tanks  nearly  empty.  The  second  loading  condition  is  with  half 
fuselage  load  and  about  half  fuel  load,  and  the  third  loading  condition 
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is  that  vrith  fuselage  nearly  e]]q>ty  and  with  tanks  full.  The  total  air¬ 
plane  weight  for  each  of  these  three  conditions  is  approximately  the  same 
so  that  the  effect  of.  load  distribution  could  be  studded.  The  fourth 
effect  considered  might  be  called  a  speed  effect;  it  deals  with  changing 
the  speed  of  the  airplane  and  is  not  meant  here  as  a  compressibility 
effect.  In  making  use  of  the  parameters  shown  in  figure  1  the 
product  XH  can  be  formulated,  which  can  be  shown  to  be  equal  to  p 


times  the  period  ratio  where  Tq 

^1 


is  the  time  required  to  reach 


maximium  gust  velocity  and  T-j^  is  one-fourth  the  natural  period  of  the 
fundamental  mode . 


This  ratio  is  formulated  because  an  analogous  period  ratio  has 
been  found  to  be  very  significant  in  response  calculations  of  elastic 
structures  under  impulsive  loadings.  The  ratio  may  be  changed  by 
changing  X,  which  is  a  function  of  V,  or  by  changing  H.  This  speed 
effect  may  thus  be  indicated  to  some  extent  by  phase  1  indicated  in 
figure  3,  which  deals  with  H,  but  to  obtain  a  direct  measure  of  speed 
itself,  all  conditions  are  held  fixed  and  only  speed  is  varied. 


To  evalviate  the  effects  outlined  in  figtire  3,  case  history  studies 
were  made  on  a  twin-engine  airplane  having  a  gross  weight  of  about 
U0,000  pounds.  The  speed  used  was  approximately  2^0  miles  per  hour, 
unless  otherwise  noted. 


Consider  the  effect  of  phase  1  -  that  of  varying  gust  gradient  dis¬ 
tance.  Fig\ire  shows  the  results  of  encoiintering  gusts  of  varying  gust 
gradient  distances,  v^ere  the  gust  was  considered  to  be  of  the  sine 
type.  The  upper  portion  of  this  figure  shows  the  variation  of  a  bending- 
moment  parameter  Cfi  with  H.  This  nondimensional  bending-moment 

parameter  can  be  used  to  calculate  the  net  incremental  moment  developed 
due  to  the  gust  from  the  equation  shown  in  the  figure,  where  k  is 
readily  calculated  from  the  lift-curw  slope,  air  density,  wing  area,  and 
distance  to  center  of  lift,  7  is  forward  velocity,  and  U  is  peak  gust 
velocity,  and  of  course  Cjj  is  picked  from  the  curve.  Values  are  shown 
of  Cjj  that  would  be  obtained  if  the  airplane  were  considered  rigid  and 

if  the  airplane  were  considered  flexible.  Thus,  the  difference  in  the 
two  curves  may  be  considered  as  a  dynamic  overshoot  effect.  The  ratio 
of  the  ordinates  of  the  curves  would  lead  to  a  response-factor  curve 
which  is  shown  in  the  lower  portion  of  the  figure.  It  can  be  seen  that 
for  gradient  distances  of  about  U  the  dynamic  overshoot  is  approximately 
20  to  25  percent,  and  for  gust  gradient  distances  of  about  10  or  greater 
the  overshoot  is  only  about  1  percent  or  less.  The  results  indicate 
that  for  the  larger  gradient  distances  the  particular  airplane  considered 
acts  essentially  as  a  rigid  body  and  that  dynamic  overshoot  is  only 
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Inqportant  in  this  range  of  distance.  These  points  will  be  discussed  In 
■ojre  detail  subseqiMntly. 

The  Influence  that  gust  shape  has  on  response  factor  Is  shown  In 
figure  5*  In  general  the  trend  Is  very  similar  for  the  three  gust  shapes 
considered  and  the  response  factor  Is  very  nearly  1  In  all  cases  for  H 
of  10  or  greater.  These  results  are  for  the  second  loading  condition 
shown  In  figure  3,  but  similar  agreement  would  be  fo\md  for  the  other 
two  conditions.  From  this  It  can  be  seen  that  the  effect  of  gust  shape 
Is  not  quite  as  critical  as  has  sometimes  been  supposed,  at  least  for 
single  gusts.  Since  shape  does  not  appear  to  be  critical,  the  sine  gust 
was  arbitrarily  chosen  for  subsequent  calculations.  . 

The  next  phase  considers  the  effect  of  loading  conditions  on  air¬ 
plane  response.  The  variation  of  the  response  factor  with  gust  gradient 
distance  for  the  three  loading  conditions  considered  Is  shown  In 
figure  6.  It  should  be  noticed  that  as  the  load  Is  removed  from  the 
center  of  the  airplane  and  Is  Increased  In  the  outboard  portion,  the 
effect  Is  to  Increase  the  response  factor  over  a  wider  range  of  gradient 
distance.  Thus,  the  effect  of  Increased  spanwlse  mass  Is  to  Increase 
the  dynamic  response  effect.  Considering  the  response  factor  alone, 
however,  may  be  misleading,  as  can  be  seen  by  considering  the  lower  half 
of  the  chart  where  Cm  Is  plotted  against  H  for  each  of  the  three 
loading  conditions.  These  curves  are  of  primary  Interest  to  the  designer 
since  they  lead  to  the  net  Incremental  bending  moment  that  would  be 
developed.  The  condition  of  full  fuselage  load  leads  to  the  largest 
value  of  Cm  and,  hence,  bending  moment,  and  there  Is  an  appreciable 

reduction  In  Cm  as  the  load  goes  outboard  on  the  span.  Thus,  the 
loading  conditions  where  load  Is  concentrated  near  the  center  line  Is 
the  critical  one,  as  would  be  expected.  The  fact  that  the  response 
factor  Is  high  for  the  outboard  loading  condition  Is,  frcmi  the  design 
standpoint,  not  Important  in  this  case. 

It  Is  of  Interest  to  consider  the  stresses  that  might  be  expected 
to  occur  during  flight.  To  do  this  It  is  worthwhile  to  review  some  of 
the  results  that  have  been  obtained  by  the  Langley  Oust  Loads  Branch 
from  aix*planes  In  thunderstorm  operation.  In  figure  7  Is  shown  a  shot¬ 
gun  pattern  of  the  variation  of  peak  gust  velocity  with  gtist  gradient 
distance  (see  reference  2),  where  H  is  given  In  feet  Instead  of  chords. 
These  data  were  gathered  from  analysis  of  time-history  records  and  were 
deduced  by  assuming  that  the  grists  encountered  were  of  a  triangular 
nature.  Assume  that  It  Is  desired  to  calculate  the  maximum  stresses 
that  would  be  developed  by  the  gusts  Indicated  here.  The  curved  line 
from  the  origin  to  U  ■  50  feet  per  second  represents  essentially  the 
envelope  curve  of  the  data  up  to  a  gradient  distance  of  about  90  feet, 
and  a  nearly  horizontal  line  Is  assumed  to  apply  thereafter.  If  given 
gradient  distances  are  considered,  gusts  which  lead  to  «*y<mnin  stresses 
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would  naturally  be  those  having  the  greatest  vertical  velocity,  or  those 
which  fall  along 'the  envelope  curve.  To  determine  the  stresses  produced 
by  the  gusts  represented  by  the  envelope  CTirve,  the  first  loading  con¬ 
dition  of  the  example  airplane  was  chosen  since  it  leads  to  the  critical 
bending  moment.  The  bending  stress  that  develops  near  the  root  of  the 
wing  as  it  varies  with  gust  gradient  distance  is  shown  in  figure  8.  The 
gra^ent  distance  is  given  here  in  both  feet  and  chord  units.  The 
results  for  the  elastic  airplane  and  the  airplane  assiuned  rigid  are  shown, 
and  these  results  tend  to  show  rather  clearly  the  effect  of  dynamic 
response.  A  rather  interesting  effect  can  be  noticed  as  gust  gradient 
distance  approaches  10  chords,  which  is  in  the  neighborhood  of  the 
nominal  value  considered  in  design.  In  this  region  maximum  stresses  are 
reached  but  dynamic  overshoot  is  very  small.  It  might  be  concluded  that, 
for  this  example,  a  rigid-body  analysis  would  be  quite  adequate  for 
design  purposes. 

It  would  be  instructive  to  f-lnd  out  how  increased  speeds  might 
affect  this  conclusion.  That  is,  consider  the  gust  which  led  to  the 
maximum  stresses  shown  in  figure  8  and  vary  the  speed  so  that  the  period 
ratio  Tq/T^  becomes  the  same  as  the  period  ratio  that  is  found  in  the 

region  of  the  figure  where  H  is  between  3  and  6  chords  and  dynamic 
overshoot  is  large.  The  first  loading  condition  was  used  and  the  air¬ 
plane  flown  with  increasing  speeds  through  sine  gusts  of  10-chord  gra¬ 
dient  distance.  The  results  are  shown  in  figure  9  where  the  dynamic 
response  factor  at  two  altitudes  is  plotted  as  a  function  of  both 
velocity  and  period  ratio  Tq/T^^.  Considering  first  the  response  factor 

at  sea-level,  it  is  seen  that  it  is  very  nearly  one  over  the  entire 
velocity  range  considered.  It  was  noted  earlier  that  when  the  period 
ratio  is  varied  by  varying  H  alone,  the  response  factor  might  change 
by  as  much  as  29  percent.  Thus  the  period  ratio  in  gust  treatment  is 
much  more  coiiq}lex  than  in  the  case  of  elastic  structures  \inder  iii9}ulslve 
loading.  The  reason  is  that  a  change  in  forward  speed  not  only  changes 
the  period  ratio  but  changes  the  aerodynamic  damping  characteristics  as 
well. 


The  results  presented  for  an  altitude  of  29,000  feet  show  an 
increase  in  the  response  factor  over  the  velocity  range  considered.  Here 
again,  considering  the  response  factor  alone  may  be  misleading;  the^ 
incremental  bending  moment  at  hOO  miles  per  hour  and  29,000  feet  would 
actually  be  about  0.6  of  the  bending  moment  at  sea  level  for  the  same 
velocity. 

The  final  results  presented  are  for  the  effect  of  flexibility  on 
accelerations.  Some  acceleration  measurements  have  been  made  on  a 
Hartln  202  airplane  in  flight  through 'rough  air  (reference  3)  and  it 
was  of  interest  to  see  whether  the  resvilts  obtained  co\ild  be  predicted 
by  theoretical  response  calculations.  To  make  the  calc\xlatlons,  the 
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gusts  shown  In  figure  7  were  again  considered.  The  results  are  shown 
in  figure  10  v^re  the  acceleration  developed  at  the  airplane  center 
line  is  plotted  against  the  acceleration  at  the  nodal  points  of  the 
fundaaental  node,  this  latter  being  used  as  a  oeasvire  of  average  airplane 
acceleration. 

According  to  theoretical  calcxilatlons  all  gusts  falling  in  the 
region  where  H  is  between  23  and  90  feet  (see  fig.  7)  would  establish 
a  band  formed  by  the  c\irve  and  the  line  shown  in  figure  10.  The  data 
obtained  in  flight  are  shown  by  the  plotted  points  in  the  lower  portion 
and  it  may  be  noticed  that  when  considered  as  a  group  the  flight  data 
merge  rather  well  into  the  band  established  theoretically. 

To  siamnarlze,  a  simplified  two-degree~of -freedom  analysis  has  been 
applied  to  a  twin-engine  transport  airplane  in  order  to  examine  the 
effect  of  wing  flexibility  on  gust  response.  Conditions  varied  were 
gust  gradient  distance ,  gust  shape,  static  load  distribution,  and  forward 
velocity.  In  the  course  of  the  investigation,  parameters  felt  to  be  of 
general  significance  for  gust  response  were  established.  Although  the 
calcxilatlons  were  made  for  a  particular  airplane,  the  trends  established 
should  be  similar  for  other  airplanes  of  the  same  general  class.  Further 
stu4y  must  be  made  of  the  effect  of  penetration  of  successive  gusts, 
spanwise  gust  variation,  compressibility,  and  probability  of  gust  occur¬ 
rence.  Also,  it  is  planned  that  the  same  general  approach  will  be  used 
to  investigate  aircraft  of  other  classes. 
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Figure  1.-  Basic  parameters  for  airplane  response  to  gusts 
iresultlng  from  two-degree-of-fireedom  analysis. 
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Figure  2.-  Cosgiarlson  of  tvo-degree>of- freedom  and  rectirrence-matrlx 
solutions  for  response  to  10- foot-per- second  sharp-R^BR  gusts. 
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Figure  5»-  Response  fiactor  for  three  gust  shapes  as  a  function 

of  gust  gradient  distance. 


Plgore  6.-  Response  factor  and  hendlng-moment  paraaster  Cff  for 
^uree  loading  conditions  as  a  function  of  gust  gradient  distance 


Figure  7.-  Gust  velocity  as  a  function  of  gust  gradient  distance* 

(See  reference  2.) 


Figure  8.-  Maxlanm  stress  as  a  function  of  gust  gradient  distance 


CCnnSBNTIAL 


RESPOIME 

nUiTOR 


V.MPH 

Is 

T| 


Figure  9*-  Response  factor  as  a  function  of  velocity  and 
period  ratio  Tq/Ti . 
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Figure  10.-  Acceleration  at  center  of  gravity  as  a  function 
of  nodal-point  acceleration. 
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OSCILLATING  PRESSURES  ON  FUSELAGES  AND 
VTENGS  DUE  TO  PROPELLERS 
Ry  Leslie  W.  Lassiter  and  Harvey  H.  Hubbard 
Langley  Aeronautical  Laboratory 


A  rotating  propeller  is  an  intense  source  of  oscillating  press\ires. 
At  large  distances  from  the  propeller  these  pressures  are  recognized  as 
merely  unpleasant  noise,  but  at  points  in  the  immediate  vicinity  of  the 
propeller  the  oscillating  pressures  are  of  such  magnitude  that  they  are 
capable  of  exciting  destructive  vibrations  in  nearby  parts  of  the  air¬ 
plane  structure.  Such  vibrations  are  excited  at  the  blade-passage 
frequency  of  the  propeller  or  at  some  integral  multiple  of  that  fre¬ 
quency.  The  fatigue  problem  associated  with  vibrations  of  this  nature 
is  a  serious  one,  since,  with  excitation  frequencies  of  the  order  of 
100  cycles  per  second,  the  structure  is  subjected  to  the  stress  cycle 
many  times  dxiring  a  single  flight.  Fatigue  of  metal  parts,  then,  is 
greatly  accelerated. 

Several  instances  of  this  type  of  vibration  problem  have  been 
reported.  There  have  been  instances  where  panels  failed  in  the  fuselage 
wall  near  the  propeller  tip  and,  in  some  cases,  wing  panels  in  the 
vicinity  of  the  propeller  have  also  failed.  The  present  paper  is  con¬ 
cerned  with  an  evaluation  of  the  oscillating  pressures  that  occur  in 
both  phases  of  the  problem.  The  discussion  is  concerned  first  with  the 
fuselage  wall  problem,  the  natxire  of  which  is  illustrated  schematically 
in  figure  1.  The  vibration  amplitudes  along  the  fuselage  wall  are  shown 
in  an  exaggerated  manner  by  the  shaded  panel  blur.  It  is  seen  that  the 
greatest  amplitudes  occur  at  about  l/U  radius  ahead  of  and  behind  the 
plane  of  rotation. 

It  is,  of  course,  desirable  from  the  designer's  standpoint  that  he 
be  able  to  predict  the  magnitudes  of  the  oscillating  pressures  so  that 
he  can  calculate  the  vibration  amplitudes  and  associated  stresses  in 
the  structure.  An  experimental  and  theoretical  investigation,  therefore, 
was  conducted  to  determine  the  significance  of  the  various  parameters 
involved  in  the  problem.  This  investigation  was  primarily  concerned 
with  pressures  in  the  region  near  the  propeller  tips  and  the  type  of 
data  obtained  is  illustrated  in  figxxre  2.  This  figure  presents  a  sample 
free-space  pressure  distribution  (at  zero  forward  speed)  of  the  fvinda- 
mental  Fourier  component  of  pressure  near  the  propeller  tip,  as  obtained 
along  the  line  parallel  to  the  axis  of  rotation  at  a  distance  of 
1/6  radius  from  the  tip.  The  root-mean-square  values  of  the  pressvire 
are  plotted  in  pounds  per  square  foot  as  a  fxinction  of  distance  in  terms 
of  propeller  radii  from  the  rotational  plane.  The  maxim\im  values  occur 
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at  distances  of  l/h  radius  ahead  of  and  behind  the  propeller  plane,  and 
the  largest  pressures  are  contained  in  a  region  about  1  radius  wide  in 
the  vicinity  of  the  propeller  plane. 

For  a  propeller  operating  in  free  space,  the  pressure  was  foimd  to 
be  a  function  of  propeller  dimensions,  torque,  thrust,  tip  Mach  number, 
and  the  point  in  space  at  which  the  measurement  was  made.  The  presence 
of  a  reflecting  surface  such  as  the  fuselage  causes  the  pressure  at  its 
surface  to  be  increased  above  the  free-space  value  at  the  same  point. 

The  amount  of  this  increase  is  a  function  of  the  fuselage  shape;  a  flat 
surface  causes  a  doubling,  whereas  a  curved  surface  increases  the  pres¬ 
sure  a  lesser  amount. 

A  theory  which  is  useful  in  predicting  these  oscillating  pressiu*es 
at  any  point  in  space  near  the  propeller  was  developed  and  is  presented 
in  reference  1.  The  agreement  of  the  theory  with  experiment  may  be  seen 
ty  comparison  of  the  curves  in  figure  2.  Good  agreement  was  found  for 
the  lower  harmonica  up  to  a  tip  Mach  number  of  1.00  which  was  the  limit 
of  this  series  of  tests. 

In  addition  to  the  problem  of  propeller-excited  vibrations  in  the 
fuselage  wall,  similar  difficulties  have  been  encountered  in  the  wing 
structure  of  pusher  configurations  like. that  shown  in  figure  3.  The 
shaded  portion  of  the  wing  indicates  the  region  where  trouble  occurred. 

In  this  area  fatigue  failures  of  the  secondary  structure  of  the  wing, 
the  wing  panels,  and  the  trailing-edge  structure  have  been  noted.  They 
have  been  observed  to  occur  not  only  in  the  region  immediately  ahead  of 
the  propeller,  but  also  in  the  region  beyond  the  propeller  tip.  Recently, 
the  work  on  fuselage  pressures  near  the  tip  of  the  blade  has  been 
extended  to  include  the  area  ahead  of  the  propeller  where  a  wing  might 
be  located.  In  this  latter  investigation  measurements  were  made  along 
lines  parallel  to  B-B  in  fig\ire  3,  in  addition  to  the  previous  measure¬ 
ments  along  lines  parallel  to  A-A.  The  experiments  were  conducted  over 
a  tip  Mach  number  range  extending  to  1.20  and  the  data  were  correlated 
with  theory  at  several  different  tip  Mach  numbers.  The  results  of  these 
tests  are  illustrated  in  figures  U  to  8. 

Figure  4  presents  the  free-space  root-mean- square  pressure  of  the 
fundamental  Fourier  component  as  a  function  of  distance  from  the  pro¬ 
peller  axis.  Pressure  units  are  po\mds  per  square  foot.  These  parti¬ 
cular  measurements  were  made  along  the  line  1/U  radius  ahead  of  the  pro¬ 
peller  at  tip  Mach  numbers  of  0.90  and  1,20.  The  maximum  value  of  the 
fundamental  Fourier  component  is  seen  to  occur  at  the  0.5  station  for 
both  tip  Mach  numbers,  but  the  magnitude  at  1.20  is  approximately  twice 
that  at  a  tip  Mach  number  of  0.90. 

Similar  graphs  of  the  distribution  of  harmonics  higher  than  the 
first  show  a  tendency  for  the  maximum  of  the  curve  to  move  progressively 
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outward  toward  the  blade  tip.  At  about  the  sixth  harmonic  the  maximum 
is  found  to  occur  even  beyond  the  tip. 

The  dashed  cttrve  at  the  top  of  figure  U  is  the  distribution  calcu¬ 
lated  ty  the  theory  for  a  tip  Hach  number  of  1.20.  It  is  seen  that, 
even  at  this  tip  Hach  nuinber,  the  theoxy  agrees  reasonably  well  with 
the  measured  distribution.  It  should  be  pointed  out  that,  for  conven¬ 
ience  in  calculating  these  pressures,  the  torque  and  thrust  forces 
acting  on  the  blade  are  assumed  to  be  concentrated  at  some  effective 
radius . 

The  effective  radius  then  is  a  function  of  the  blade  loading  and 
the  manner  in  vdilch  the  forces' at  each  blade  element  contribute  to  the 
free-space  press\ires  at  a  point  in  space  for  a  given  harmonic.  The 
value  of  the  effective  radius  thus  may  differ  for  the  various  harmonics 
and  for  different  blade  loading. 

For  the  fiindamental  frequen<y  the  effect  of  increasing  distance 
from  the  propeller  is  beneficial  in  decreasing  the  free-space  pressures. 
This  attenuation  effect  is  Illustrated  in  figure  5*  The  curve  on  the 
right  indicates  the  attenuation  with  distance  from  the  tip  in  the 
z  direction,  whereas  the  attenuation  with  distance  in  the  y  direction  is 
shown  b7  the  curve  at  the  left.  The  pressure  amplitudes  of  the  funda¬ 
mental  component  are  thus  seen  to  be  attenuated  rapidly  as  distance  is 
Increased.  There  is  some  evidence,  however,  to  indicate  that  the  effect 
of  distance  may  be  less  beneficial  in  the  case  of  the  higher-order 
harmonics  which  have  spatial  distributions  different  from  that  of  the 
fundamental  Foiirier  component. 

As  was  stated,  the  curves  shown  in  fig^lres  2  to  5  relate  only  to 
the  fmidamental  component  of  pressure.  It  must  be  remembered  that  the 
total-pressure  variation  Is  generally  composed  of  many  of  these  Fourier 
conqponents,  all  Integral  multiples  of  the  blade-passage  frequency.  This 
is  especially  true  near  the  blade  tip  because  of  the  tendency  for  the 
higher  harmonics  to  peak  in  that  region;  therefore,  the  total-pressure 
variation  with  time  is  not  by  any  means  a  sinusoidal  one.  This  varia¬ 
tion  is  Illustrated  in  figure  6,  which  shows  a  sample  wave  form  of  the 
total  pressure  in  free  space. 

This  wave  form  represents  a  time  history  of  the  pressure  at  a  point 
in  space  during  one  revolution  of  a  two-blade  propeller.  Each  pair  of 
peaks,  positive  and  negative,  corresponds  to  a  blade  passage.  Positive 
pressure  is  measured  upward  from  the  reference  line;  negative  pressure 
is  downward.  Although  the  particular  wave  form  shown  was  obtained  from 
a  free-space  measurement,  its  general  shape  applies  also  to  wave  forms 
observed  at  a  reflecting  surface.  Neither  the  free-space  nor  surface- 
pressure  wave  forms  are  necessarily  symmetrical  about  the  reference 
pressure  as  this  one  is;  in  general  they  are  not  symmetrical. 
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When  a  wing  is  introduced  into  the  oscillating  pressure  field,  it 
is  acted  upon  by  the  pressure  in  two  ways.  First,  the  phase  and  n^ni- 
tude  differences  between  pulses  impinging  on  the  top  and  bottom  s\irfaces 
of  a  wing  cause  an  alternating  differential  pressure  to  exist  across  the 
wing.  Near  the  trailing  edge,  where  the  structure  is  relatively  thin, 
this  differential  pressure  tends  to  flex  the  structure.  Figure  7  illus¬ 
trates  the  variation  of  differential  pressxire  across  a  thin  wing  as  a 
function  of  distance  from  the  propeller  axis  and  gives  a  con^arison  with 
free-space  peak-to-'  ik  values  at  the  same  location.  The  wing  surface 
used  consisted  of  a  thin  flat  board  without  camber.  The  values  plotted 
in  this  figure  are  the  peak-to-peak  values  of  the  total  pressure  obtained 
from  wave  fonns  similar  to  that  of  figure  6.  It  may  be  that  the  total 
pressure  is  not  of  as  much  importance  as  its  individual  components,  but 
it  is  thought  that  discussion  of  the  total  pressure  is  Justified  hy  the 
large  differences  in  magnitude  and  distribution  existing  between  it  and 
its  conqponents. 

It  is  seen  that  the  pressure  differential  across  a  thin  wing  is 
much  greater  than  the  free-space  pressure  at  the  same  point.  The  two 
have  similar  spatial  distributions,  however,  with  the  maximum  pressure 
occurring  in  the  region  near  the  blade  tip  and  slightly  past  the  tip. 

It  will  be  remembered  that  the  fvindamental  component  of  pressxire  peaked 
at  the  0.5  station j  however,  the  higher  harmonics,  of  which  many  are 
present,  peak  further  out  along  the  blade  to  give  the  resultant  total 
differential  pressvire  distribution  shown  in  the  figure.  It  is  this  dif¬ 
ferential  pressxire  which  may  have  been  responsible  for  structural 
failures  near  the  wing  trailing  edge,  since  it  imposes  oscillating  loads 
which  contribute  to  accelerated  fatigue. 

Another  factor  to  be  considered  in  conjunction  with  the  wing 
structure  is  the  pressure  acting  on  each  surface  individually.  It  is 
this  pressure  which  affects  the  wing  panels  and  their  associated  members. 
There  is  a  large  difference  in  the  magnitudes  of  these  panel  pressures 
for  a  given  wing  surface,  depending  upon  whether  the  propeller  blade 
approaches  the  surface  or  recedes  from  it.  Figure  8  illustrates  the 
variation  of  the  surface  pressure  with  distance  along  the  wing  for  both 
top  and  bottom  surfaces  as  the  blade  passes.  The  values  plotted  here 
are  obtained  from  the  positive  and  negative  amplitudes  of  the  total- 
pressure  wave  forms,  as  illustrated  in  figure  6.  The  curves  at  the  top 
relate  to  the  surface  approached  by  the  blade  (in  this  illustration  the 
top  surface);  those  at  the  bottom  of  the  chart  relate  to  the  opposite 
STorface,  from  which  the  blade  is  receding.  The  peak  pressures  in  both 
cases  are  plotted  in  units  of  pounds  per  square  foot,  con^ression  or 
suction. 

From  the  curves  at  the  top,  it  is  seen  that  where  the  surface  is 
approached  by  the  blade,  the  pressure  variation  with  time  is  essentially 
symmetrical  about  the  reference  or  atmospheric  pressure.  Wing  panels  on 
that  surface,  then,  are  subjected  alternately  to  large  compression  and 
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suction  loads.  Both  positive  and  negative  magnitudes  are  greatest 
opposite  the  propeller  tip,  and  they  remain  quite  high  even  beyond  the 
tip.  It  will  be  remembered  that  the  trouble  area  indicated  in  figure  3 
extended  past  the  tip  also.  If  the  magnitudes  on  the  surface  from  which 
the  blade  receded  are  cooqpared  with  those  cn  the  opposite  surface,  it  is 
seen  that  the  sxirface  approached  by  the  blade  is  subjected  to  greater 
pressures  and  may  be  expected  to  prove  most  critical  structurally.  It 
is  in^ortant  to  note  that  the  direction  of  propeller  rotation  is  signifi¬ 
cant  since  that  determines  which  surface  is  subjected  to  greater  pres¬ 
sures.  Reversal  of  propeller  rotation  in  the  diagram  of  figure  8  would 
result  in  making  the  bottom  surface  most  critical. 

In  summary:  two  phases  of  the  propeller-excited  vibration  problem 
have  been  presented.  One  phase  pertains  to  fuselage  pressures  in  the 
region  near  the  propeller  tip;  the  other  pertains  to  pressTires  on  the 
wing  structure  near  the  propeller.  It  has  been  shown  in  this  preliminary 
investigation  that  these  oscillating  pressures  may  be  calculated  and 
hence  the  order  of  magnitude  of  the  problem  can  be  defined. 
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Figure  3*-  Observed  area  of  fatigue  failures  In  pusher 

configurations. 


Flgnre  1^.-  Free-space  pressure  distributions  of  the  fundsMntal 
Fourier  cosqponsnt  ahead  of  a  propeller. 
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Figure  Free-epace  pressures  of  the  fundaaental  Fourier 
conponent  as  a  function  of  distance  front  the  propeller. 
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Figure  7.*  Total-preaaure  dlatrlbutlon  for  free  space  coo^ared  to 
differential  preaaurea  acroaa  a  thin  vlng. 


Flfore  8.-  Distribution  of  total  pressure  at  the  surfaces  of  a 
thin  vlsg  during  a  propeller  blade  passage. 


CCMFIDSITIAL 


COHFirOTlAL 


27 


WING-ANALTSIB  PR0BLK1I8 


CORFIOIRTIAL 


COHFIDEIRrUL 


29 


DBFUCTIQIIS  OF  SOLID  N  OR  W  ASD  SWEPT  WINGS 
Bor  Oeorse  V.  Zender  and  Richard  R.  Heldenfels 
Langley  Aeronautical  Laboratory 


Thin  high-epeed  aircraft  vings  nay  be  swept  and  have  solid  or  nearly 
solid  construction.  Swept  wings  of  fairly  large  aspect  ratio  (aspect 
ratios  at  least  greater  than  2)  resesible  ordinary  cantilever  beaas  out> 
board  of  the  triangular  region.  This  suggests  that  ordinary  beam  theory 
with  slight  modification  nay  be  used  to  calculate  the  defomatlons . 

An  eqieriaental  study  is  therefore  being  undertaken  to  determine  the 
extent  to  which  bean  thecnry  may  be  used  for  calculating  the  defomatlons 
of  solid  wings  of  swept  and  N  or  W  plan  forms.  The  first  models  tested 
in  this  program  are  shown  in  figure  1.  The  models  had  a  rectangular 
cross  section,  4^^  sweep  angle,  were  of  248-T4  aluminum  alloy,  and  were 
full -span  models  with  a  center  section.  The  center  section  was  clamped 
between  two  tui^sten  carbide  blocks  and  loads  were  applied  to  both  halves 
of  each  model  in  order  to  minimize  the  effect  of  Jig  defomatlons. 
Systematic  deflection  data  were  taken  on  one  side  only,  but  two  points 
were  checked  on  the  other  side. 

Figure  2  shows  the  deformations  of  the  swept  model  subjected  to 
vertical  tip  loads.  The  left-hand  figure  shows  the  vertical  deflection 
of  the  center  line  and  the  right-hand  figure  shows  the  twist  of  the  wing. 
Two  types  of  twist  are  shown.  The  first  is  the  twist  about  the  center 
line  which  night  be  called  the  structural  twist  and  is  designated  by  the 
aynhol  6.  The  other  is  the  twist  measured  in  a  plane  parallel  to  the 
root  which  is  equivalent  to  the  angle  of  attack  and  which  is  designated 
by  the  symbol  a.  The  angle  of  attack  a  contains  components  of  the 
twist  6  and  the  bending  slope.  Figure  2  shows  three  theoretical  curves 
which  Illustrate  the  results  that  can  be  obtained  from  ordinary  beam 
theory,  using  various  assumptions  for  the  effective  root.  The  root 
designated  by  the  length  Li  measured  adong  the  center  line  results  in 
too  little  deflection}  the  distance  Lg,  which  is  the  length  corresponding 
to  the  Intersection  of  the  center  lizie  and  the  clamped  root,  results  in 
too  much  deflection.  These  two  effective  roots  are  the  extremes  which 
ml^it  result  from  an  arbitrary  selection.  The  figure  also  shows  that  an 
effective  root  one-third  the  distance  trcn  Li  to  L2,  designated  by  L, 
gives  very  good  agrement  with  the  test  data  for  deflections.  It  is 
noticed,  however,  that  the  use  of  L  slightly  overestimates  the  angle- 
of -attack  change.  This  results  Arom  using  ordinary  beam  theory  which 
neglects  the  small  amount  of  structural  twist  produced  by  bending  loads, 
as  shown  by  the  square  test  points  in  the  right-hand  figure.  The 
component  of  this  twist  parallel  to  the  root  offsets  slightly  the  angle- 
of -attack  change  due  to  bending  of  the  model. 
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Many  tests  were  made  under  vertical  loads  to  determine  the  variation 
of  the  effective  root  location  with  the  point  of  load  application.  Fig¬ 
ure  3  suonmurlzes  the  results  of  these  tests.  The  location  of  the  effec¬ 
tive  root  a/b  Is  plotted  as  a  function  of  the  ratio  l/c.  The  distance 
a  Is  measured  from  the  rear  corner  of  the  sweptback  wing  to  the  effective 
root,  and  b  Is  the  distance  from  the  rear  corner  to  a  root  at  the  Inter¬ 
section  of  the  center  lines  of  the  outer  part  and  center  section.  The 
distance  along  the  center  line  to  the  point  of  application  of  the  load 
Is  designated  by  Z,  and  c  Is  the  structural  chord.  The  circles  give 
the  test  results  and  the  solid  line  has  been  faired  through  the  test 
points,  note  that  as  Z/c  becomes  larm,  that  Is,  bending  moment 
predcAlnates  over  shear  at  the  root,  a/b  approaches  a  constant  value. 

The  two  dashed  lines  are  estimated  locations  of  the  effective  root  for 
sweep  angles  of  30°  and  60°.  These  estimates  are  based  on  test  data 
obtained  from  reference  1.  This  figure  may  also  be  used  for  a  distributed 
loading  and  the  effective  root  obtained  considering  Z  as  the  distance 
to  the  centroid  of  the  loading. 

A  similar  Investigation  was  made  to  determine  the  effective  root 
auisoclated  with  a  torque  about  the  center  line.  The  twists  produced  by 
a  tip  torque  about  the  center  line  are  sho%m  In  figure  4.  It  was  found 

that  the  effective  root  coxild  be  taken  at  L  ^  that  Is,  ^  «  1,  and 

the  twist  e  calculated  from  elmaentary  torsion  theory.  The  solid  line 
shows  the  twists  calculated  with  this  effective  root.  The  use  of 
different  effective  roots  for  bending  and  torsion  Is  Inconvenient,  but  In 
most  wing  deformation  analyses  the  torques  are  small  compared  with  the 
bending  moment  and  very  little  error  Is  Incurred  by  using  the  effective 
root  for  bending  for  the  torques.  In  all  the  subsequent  comparisons 
between  experiment  and  calculations,  the  calculations  were  made  on  this 
basis. 

Figure  ^  shows  the  deformations  produced  by  a  couple  applied  to  the 
tip  of  the  swept  model.  The  plane  of  the  couple  Is  parallel  to  the  root} 
thus,  effectively  It  Is  a  combination  of  a  torque  and  a  bending  moment. 

TZie  deformations  consist  of  deflections  and  rotations  a  and  0.  Calcu¬ 
lations  are  In  good  agreement  with  the  test  data.  Effective  root  location 

a  1 

was  deteralned  by  taking  g  *  which  Is  the  asymptotic  value  shown  In 
figure  3* 

Figure  6  shows  similar  results  when  the  couple  Is  applied  at  the  half 
span.  Note  that  outboard  of  the  applied  load,  the  rotations  a  and  6 
are  constant  and  the  deflection  varies  linearly  as  would  be  expected  for 
rigid  body  movements .  Cedculatlons  are  again' In  good  agreement  with 
tests. 
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The  preceding  methods  for  the  deformation  analysis  of  swept  wings 
can  be  applied  to  M  or  W  wings  If  the  Inner  and  outer  portions  are 
treated  as  Individual  swept  wings.  Figure  7  shows  the  results  for  a 
vertical  load  applied  at  the  tip  of  a  V  model.  Again  the  left-hand  fig¬ 
ure  shows  the  vertical  deflection  of  the  center  line  and  the  right-hand 
figure  shows  the  twist.  The  plane  of  the  structural  twist  d  Is 
different  for  each  part  of  the  wing}  thus^  there  Is  a  discontinuity  In 
the  curve  at  the  midspan  station.  Fairly  good  agreement  Is  obtained 
between  calculations  and  tests,  hut  not  quite  as  good  as  for  the  swept 
case.  The  procedure  for  analyzing  the  Inner  portion  of  the  model  Is  the 
same  sis  that  previously  described  for  the  swept  model.  The  deformations 
of  the  outer  part  are  obtained  from  superposition  of  rigid  body  motions 
due  to  deformations  of  the  Inner  part  and  deformations  of  the  outer  part 
due  to  the  loads  on  the  outer  part.  The  latter  was  calculated  using  an 
effective  root  determined  from  figure  3*  Note  that  for  this  model,  the 
Inner  part  Is  under  combined  bending  and  torsion  which  causes  a  decrease 
In  the  angle  of  attack  near  the  root,  but  farther  outboard  the  angle  of 
attack  Increases. 

Figure  8  shows  similar  results  for  a  vertical  load  applied  at  the 
half  span  of  a  W  model.  The  location  of  the  effective  root  a/b  was 
k/lO  for  the  effective  root  at  the  inboeurd  end  of  the  model.  No  effective 
root  was  necessary  for  the  outer  panel  since  no  loads  were  applied  in 
this  region.  The  calculated  deflection  and  the  rotations  a  and  9  are 
constant  In  the  outer  part  as  would  be  expected  from  rigid  body  movements. 
The  actual  deflection  vcLries  slig^itly  because  of  coupling  between  bending 
amd  twisting.  Note  that  the  angle  of  attack  is  negative  at  all  stations 
along  the  span.  Since  the  tip  and  half-span  loads  have  opposite  effects 
on  the  angle -of -at tack  change,  the  structural  designer  has  some  control 
over  the  angle-of- attack  changes  of  M  or  W  wings.  For  swept  wings, 
however,  all  vertical  shear  loads  cause  a  decrease  in  the  angle  of  attack 
and  no  such  control  Is  possible. 

Figure  9  shows  results  for  a  couple  applied  at  the  tip  parallel  to 
the  plane  of  the  root.  Both  parts  of  the  model  are  under  combined  bending 
and  torsion  which  cause  both  deflections  and  rotations.  At  the  midspan 
station,  where  the  direction  of  the  loads  change,  the  deflection  shows  a 
change  In  slope  and  9  has  a  discontinuity. 

Figure  10  shows  similar  results  when  a  couple  Is  applied  at  the  half 
span.  In  the  outer  part,  which  Is  not  loaded,  a  and  9  are  constant 
and  the  deflection  varies  llne€u*ly. 

In  conclusion  it  has  been  shown  that  the  deformations  of  swept 
and  M  or  W  wings  can  be  calculated  by  elementary  beam  theory  with  effec¬ 
tive  roots.  The  accuracy  has  been  demonstrated  for  45°  specimens  and 
several  types  of  arbitrary  loads.  The  material  presented  is  the  initial 
phase  of  an  investigation  which  will  ultimately  include  the  effective 
root  location  for  other  sweep  angles,  plan  forms,  and  cross  sections. 
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Figure  1.-  Svept  and  M  or  W  models. 


2.-  Deflections  and  tvlsts  of  svept  model  for  vertical 

tip  load. 
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Figure  3.-  Effect  of  load  position  on  location  of  effective 

root. 
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Figure  4.-  Twist  of  swept  model  for  tip  torque  about  center  line 
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Figture  5.-  Deflections  and  twists  of  swept  model  for  tip  couple. 


Figure  6.-  Deflections  and  twists  of  swept  model  for  half-span  couple 
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Flgtore  7»-  Deflections  and  twists  of  W  model  for  vertical  tip  load 
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Figure  8.-  Deflections  and  twists  of  W  model  for  vertical 

half- span  load. 


CONFIMRTIAL 


CCHfiiaBTlAL 


Figure  9»-  Deflections  and  twists  of  W  model  for  tip  couple 


Figure  10. •  Deflections  and  twists  of  W  andel  for  half-span  couple 
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DEFLECTIONS  AND  STRESSES  OF  SOLID  DELTA  WINGS 
^  Nanuel  Stein  and  J.  Edward  Anderson 
Langley  Aeronautical  Laboratory 


The  previous  paper  by  George  W.  Zender  and  Richard  R.  Heldenfels 
has  discussed  the  application  of  beam  theory  to  thin  solid  wings  of  high 
aspect  ratio.  For  thin  solid  wings  of  low  aspect  ratio  it  is  no  longer 
adequate  to  use  beam  theory  to  study  structural  deformation.  Wings  of 
this  type  are  more  nearly  plates  than  beams  and  should  be  analyzed  ty 
plate  theory.  Plate  theory,  however,  utilizes  a  fourth-order  partial 
differential  equation,  solutions  to  which  are  not  readily  obtained, 
especially  for  plates  of  variable  thickness  under  arbitrary  loading. 

In  a  paper  being  prepared  by  ^he  NACA  there  is  presented  an  approximate 
plate  theory  based  on  a  method  sviggested  by  Professor  Eric  Reissner  for 
analyzing  the  distortion  of  such  wings.  (Essentially  the  same  theory 
derived  in  a  different  manner  has  appeared  in  a  very  recent  Swiss  paper 
(see  reference  1).)  The  present  paper  describes  this  theory  and  its 
applif  ation  to  delta  wings . 

Figure  1  shows  a  thin  cantilever  delta  wing  of  arbitrary  thickness 
variation.  A  natural  assumption  from  the  point  of  view  of  structural 
deformation  and  one  that  is  convenient  from  the  point  of  view  of  calcu¬ 
lating  aerodynamic  loading  due  to  structural  deformation  is  that  the 
deflection  of  any  cross  section  parallel  to  the  root  consists  of  only  a 
translation  and  a  rotation.  Thus,  if  the  spanwise  variation  of  the 
translation  of  the  trailing  edge  is  denoted  by  W(x)  and  the  spanwise 
variation  of  the  rotation  is  denoted  by  0(x),  the  deflection  at  ary 
point  in  the  wing  is  given  by 

w  »  W(x)  +  y  0(x)*~^  (1) 

The  right-hand  side  of  equation  (1)  can  be  regarded  as  the  first  two 
terms  of  an  expansion  of  w  in  powers  of  y.  For  greater  accuracy  a 
third  term,  y2a(x),  which  represents  a  parabolic  camber  deformation 
may  be  added  to  the  series.  It  might  be  mentioned  that  in  assuming  the 
displacement  of  any  cross  section  to  be  made  up  of  a  translation  and  a 
rotation  nothing  new  is  being  done.  What  is  new,  however,  is  that  these 
quantities  are  being  calcxilated  more  accurately  by  regarding  the  wing 
as  a  plate  rather  than  as  a  beam. 

The  next  step  in  the  analysis  is  to  choose  the  functions  W  and  6 
so  as  to  minimize  the  potential  energy  of  the  system.  This  results  in 
two  ordinary  differential  equations  for  W  and  0,  which  are  given  as 
follows: 

(aiW»)«  ♦  (a20")''  -  p^  (2) 
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where 


(a2W")"  ♦  (a30")"  -  2(1  -  |i)(ai6‘)'  •  Pg 


(3) 


and 

•»v3 

D  ■  - r  ■  Plate  stiffness 

12(l  -  1.2) 

p  ■  P(x,  y)  ■  Lateral  pressure 
h  «  h(x,  y)  "  Local  thickness 

The  primes  denote  differentiation  with  respect  to  x.  The  a's  are 
functions  only  of  the  plate  material  and  configuration  and  are  defined 
by  the  integrals  shown.  The  right-hand  sides  of  these  equations,  p^ 

and  pg,  are  the  running  lateral  load  and  running  torque,  respectively, 

on  sections  parallel  to  the  root  and  are  given  by  chordwise  Integrals 
as  shown  here. 

With  equations  (2)  and  (3)  solved  for  W  and  9  the  deflections 
of  the  wing  are  known.  The  stresses  can  then  be  determined  by  suitable 
differentiation  of  the  deflections.  The  theory  presented  herein  so  far 
applies  to  solid  or  nearly  solid  thin  wings  under  any  loading  of  any 
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plan  form,  not  only  delta  wings.  In  what  follows  some  specific  problems 
involving  delta  wings  for  which  these  eqiiations  have  been  solved  exactly 
are  discussed  and  coiq>arisons  ajre  made  with  experiment. 

Figure  2  shows  a  coiqsarison  between  theory  and  experiment  for  the 
deflections  of  a  triangular  plate  of  uniform  thickness  under  uniform 
load.  The  circles  show  the  experimentally  measured  deflections  along 
the  chord  at  the  0.2,  O.li,  0.6,  0.8,  and  tip  spanwise  stations.  The 
dashed  lines  show  the  deflections  at  these  stations  calculated  by  means 
of  equations  (2)  and  (3).  If  the  theoiretlcal  analysis  is  refined  by  the 
addltlcm  of  the  parabolic'  camber  term  previously  mentioned,  the  solid 
curves  az^  obtained  vdilch  are  in  somevdiat  better  agreement  with  the  test 
results. 

Figure  3  shows  a  similar  comparison  for  a  60^  triangular  plate. 

Again  the  circles  are  experimentally  measured  deflections,  the  dashed 
curves  were  obtained  by  use  of  equations  (2)  and  (3)  and  the  solid  curves 
by  the  analysis  refined  by  the  addition  of  the  parabolic  camber  term. 

Now  let  us  look  at  the  stresses.  Figure  U  shows  the  con^rison 
between  theory  and  experiment  for  a  triangular  plate  of  uniform 
thickness  under  \inlform  load.  The  circles  give  the  experimentally 
determined  maximum  principal  stresses  along  the  chord  very  near  the  root 
and  at  the  0.2,  O.U,  and  0.6  spanwise  stations.  Again  the  dashed  curves 
give  the  results  of  calculations  based  on  linear  chordwise  deflections 
and  the  solid  curves  Include  the  effect  of  the  parabolic  camber  term. 
Vfhlle  the  linear  approximation  gave  fairly  accurate  deflections,  it 
appears  from  this  comparison  that  the  inclusion  of  the  parabolic  camber 
term  is  necessary  for  satisfactory  stress  calculations.  The  maximum 
stress  occurs  at  the  root  at  the  trailing  edge.  It  may  be  of  interest 
to  remark  that  this  stress  is  55  percent  higher  than  the  stress  given 
by  the  usual  Nc/I  formula  of  beam  theory.  Thus,  in  a  delta  plate 
there  is  a  stress  concentration  analogous  to  the  stress  build-up  at  the 
root  of  the  rear  spar  in  conventional  sweptback  wings. 

For  the  purposes  of  examining  theoretically  the  effects  of  change 
in  cross-section  shape  and  spanwlse-thlckness  taper,  an  exact  theoretical 
solution  to  the  equations  presented  has  been  obtained  for  a  triangular 
plate  of  diamond  chordwise  cross  section  and  linear  spanwise  taper  under 
uniform  loading.  The  results  for  deflections  along  the  leading  and 
trailing  edges  are  shown  in  flgtire  5  uid  for  compaHson  there  are  shown 
the  deflections  along  the  leading  edge  and  along  the  trailing  edge  of  a 
uniformly  thick  triangular  plate  of  the  same  weight  under  uniform 
loading.  The  tapered  plate  deflects  considerably  less,  but  the  tip 
tends  to  curl  up  uxMier  the  uniform  loading. 

In  figure  6  theoretical  stress  curves  are  presented  for  the  two 
li$o  triangular  plates  of  the  same  weight  under  uniform  loading  as  was 
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done  for  deflections.  The  ciirves  at  the  left  show  the  stresses  along 
the  root  chord.  The  curves  at  the  right  show  stresses  along  the  chord 
at  60-percent  span.  The  maxlmiim  stress  In  the  diamond  cross-section 
plate  Is  well  below  that  of  the  uniform  plate.  Values  of  stress  close 
to  the  maximum  value  appear  all  along  the  span  of  the  diamond  cross- 
section  plate. 

Figure  7  presents  the  results  of  an  experimental  Investigation  and 
a  theoretical  study  of  the  natural  vibration  of  a  triangular  canti¬ 

lever  plate.  This  figure  shows  the  shape  of  the  nodes  In  the  Interior 
of  the  plate  for  the  first  three  natural  modes  of  vibration  as  determined 
by  experiment.  The  upper  nun^rs  give  the  experimentally  determined 
values  of  the  frequencies.  The  lower  numbers  are  the  frequencies  esti¬ 
mated  by  an  approximate  solution  to  the  equations  of  the  present  theory 
by  an  energy  method  assuming  only  six  degrees  of  freedom. 

In  general  it  may  not  be  possible  to  obtain  closed-form  solutions 
of  the  differential  equations  presented  for  wings  navlng  certain  desired 
airfoil  section^  taper^  and  loading.  In  such  cases  numerical  methods 
may  be  used.  To  check  the  accuracy  of  the  numerical  methods  for  delta 
wings  using  the  present  theory,  the  case  of  a  hS^  triangular  plate  of 
uniform  thickness  under  \inlform  load  was  solved  by  a  numerical  method 
which  considers  five  stations  along  the  span  and  which  utilizes  the 
sinq)le  quadratic  differentiation  approximation.  A  comparison  between 
this  numerical  method  and  the  closed-form  solution  of  the  differential 
equations  is  shown  in  figure  8  where  the  curves  are  the  results  of 
closed-form  solution  of  differential  equations  (2)  and  (3)  and  the 
sqxiares  are  the  results  of  nuzaerlcal  solution  of  the  same  equations. 

The  results  are  as  good  as  can  be  expected.  The  solution  involves  five 
simultaneous  linear  algebraic  equations  for  the  linear  chordwise  deflec¬ 
tion  theory.  For  the  parabolic  camber  theory  it  would  be  necessary  to 
solve  10  equations. 

A  method  of  analyzing  thin  solid  or  nearly  solid  wings  has  been 
presented.  This  paper  has  shown  that  accurate  results  can  be  obtained 
by  this  method  for  triangular  plates  of  uniform  thickness  under  uniform 
load.  It  has  also  been  indicated  that  the  method  is  readily  applicable 
to  triangular  plates  of  other  than  uniform  thickness  and  that  the  equa¬ 
tions  of  this  theory  may  be  readily  solved  by  numerical  methods. 
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Figure  1.-  Coordinate  ayatem  and  aaaumed  deflection  function 
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Figure  2.-  Defleetlona  of  a  triangular  plate  of  unlfom 
thiekMaa  unAer  unlfoim  load.  ^ 
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Figure  7*-  First  three  modes  of  vibration  of  a  uniform 

triangular  plate. 
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Figure  8.-  Campari  son  of  deflections  of  a  triangular  plate 
obtained  by  nmwrieal  solution  and  by  exact  solution. 
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STABILIZATION  OF  COMPRESSION  COVERS 
OF  SHELL  WINQS  BY  POSTS 
By  Paul  Selde  and  Paul  F.  Barrett 
Langley  Aeronautical  Laboratory 


Some  Interest  has  been  displayed  In  the  \ise  of  vertical  posts  as  a 
means  of  stabilizing  the  thick  compression  cover  of  thin  wings.  In 
order  to  obtain  Information  that  mlc^t  serve  as  a  basis  for  design,  a 
theoretical  buckling  analysis  of  this  type  of  construction  was  made. 

The  theoretical  analysis  was  confined  to  the  type  of  structure 
shown  In  figure  1.  The  stziicture  considered  Is  a  long  box  beam  sub¬ 
jected  to  pure  bending.  The  box  beam  consists  of  top  and  bottom  covers 
(not  necessarily  of  the  same  thickness),  two  side  webs  which  are  not 
shown,  and  vertical  posts  connecting  the  top  and  bottom  covers.  The 
chordwlse  and  spanwlse  spaclngs  of  the  posts  are  uniform;  however  they 
may  differ  from,  each  other.  Althovi^pi  only  two  spanwlse  rows  of  posts 
are  shown  here,  the  analysis  applies  to  a  box  beam' with  any  nximber  of 
rows.  In  analyzing  this  structure.  It  was  euBsumed  that  the  covers 
remained  flat  and  that  the  side  webs  provided  them  with  simple  support. 
The  posts  themselves  were  assumed  not  to  restrict  the  rotation  of  the 
covers  at  the  points  of  attachment. 

In  the  computations  attention  was  restricted  to  the  problem  of 
detenalnlng  the  poet  axial  stiffness  reqvilred  for  the  poets  to  8u:t  as 
ribs  or  shear  webs  idien  buckling  occurs;  so  that  either  transverse  or 
longitudinal  nodes  throu^^  the  posts  could  be  obtained.  With  this 
axial  stlffhess  the  posts  act  as  thou^  they  are  rigid  and  no  advante^ 
results  from  Increasing  the  stiffness.  The  results  of  the  analysis 
and  coiqiutatlons  that  apply  to  the  partlciilar  case  of  covers  of  equal 
thickness  are  shown  In  figure  2.  The  ordinate  Is  the  required  post 
axial  stiffness  and  the  abscissa  Is  the  ratio  of  longitudinal  distance 
between  posts  to  cover  width  L/t,  To  see  the  meaning  of  the  curves, 
consider,  for  example,  the  curve  which  applies  for  a  box  with  one  span- 
wise  row  of  posts.  For  any  combination  of  L/b  and  post  axial  stiffness 
that  falls  on  or  to  the  right  of  the  curve,  buckling  will  occur  with 
nodes  throtigh  the  posts.  In  this  case  with  transverse  nodes,  where  a 
chordwlse  line  of  posts  acts  as  a  rib.  The  buckling  stress  coefficient 
for  transverse  nodes  corresponding  to  the  particular  value  of  L/b  Is 
given  at  the  top  of  the  chart.  For  any  combination  falling  to  the  left 
of  the  vertical  asymptote  to  this  curve  It  Is  not  possible  to  obtain 
rib  action,  and  buckling  will  always  occur  with  vertical  movement  of 
the  posts.  Rib  action  la  obtainable  at  higher  buckling  stress  coeffi¬ 
cients  If  more  rows  of  posts  are  used.  Similar  sets  of  curves  were 
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obtained  for  various  values  of  the  ratio  of  tension-cover  to  compression- 
cover  flexural  stiffness  and  are  presented  in  reference  1  where  numerical 
values  of  post  stiffhess  are  given. 

A  sunmary  of  the  buckle  patterns  obtainable  by  using  one  spanvlse 
row  of  rigid  or  effectively  rigid  posts  is  given  in  figure  3*  Similar 
charts  are  available  in  reference  1  for  two  and  three  rows  of  posts. 

The  three  regions  of  the  chart  correspond  to  three  different  modes  of 
buckling.  For  values  of  L/b  and  the  ratio  of  the  flexural 

stiffness  of  the  tension  cover  to  that  of  the  compression  cover,  that 
fall  in  region  I,  buckling  occurs  with  rib  action  along  every  chord- 
wise  line  of  posts;  in  region  II ^  with  spar  action  along  every  sx>anwlse 
row;  and  in  region  III,  with  vertical  movement  of  the  posts.  The  dashed 

Dp 

lines  represent  the  nodes  that  occur.  Taking  the  value  of  —  =  1,  that 

Oc 

is,  covers  of  eqxial  thickness,  it  can  be  seen  that  buckling  with  spar 
action  cannot  occur  and  that  either  buckling  with  translation  of  the 
poets  or  rib  action,  depending  on  the  post  spacing,  is  obtainable.  It 
is  apparent  that,  for  ratios  of  cover  flexural  stiffness  usually  encoun¬ 
tered  in  aircraft  construction,  spar  action  with  posts  cannot  be  obtained. 

In  order  to  obtain  some  idea  of  the  behavior  under  load  of  an 
actual  box  beam  stiffened  by  posts,  the  specimen  shown  in  figure  4  was 
designed  and  tested.  The  nominal  thickness  of  the  covers  was  3/8  inch, 
the  box  depth  was  about  6  inches,  and  the  box  width  was  28  inches. 

Three  rows  of  ^-inch-diameter  posts,  spaced  6  Inches  in  the  spanwise 

direction,  were  uised.  The  weight  of  the  posts,  including  all  attach¬ 
ment  material,  was  about  6  percent  of  that  of  the  compression  cover. 
According  to  the  theory,  the  box  would  develop  a  cover  buckling  stress 
of  U0,0C)0  psl,  and  the  compression  cover  would  buckle  with  transverse 
nodes  throu^  the  posts. 

The  box  beam  was  mounted  in  the  combined  load  testing  machine  of 
the  Lsngley  structures  research  laboratory  and  subjected  to  a  pure 
bending  moment  which  put  the  upper  sTirface  in  compre'sslon.  Strains 
were  measured  at  three  locations  on  the  upper  surface  of  the  top  cover 
to  serve  as  a  basis  for  detecting  buckling.  The  growth  of  the  strains 
with  applied  bending  moment  is  shown  in  figure  If  the  structure 
behaved  as  an  ordinary  beam,  the  strains  would  be  expected  to  increase 
linearly  with  moment.  The  noticeable  departure  of  these  curves  from 
linearity  is  due  to  secondary  bending  effects  which  included  sn  over¬ 
all  downward  dishing  of  the  two  covers  relative  to  their  edges.  This 
over-all  dishing  continued  to  grow  with  load  until  it  became  of  the 
order  of  the  cover  thickness  and  then,  rather  sviddenly,  there  was  super- 
iaqposed  on  it  the  bxickle  pattern  Indicated  by  the  theoretical  analysis, 
namely,  a  short-wave  pattern  with  transverse  nodes  through  the  posts. 
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When  this  happened,  the  strains  Indicated  by  gages  1  and  3  reversed. 

The  biding  moment  at  vhlch  this  strain  reversal  took  place  vas  taken 
as  the  experimental  buckling  load  and  vas  about  2  percent  below  the 
theoretically  ccmiputed  moment.  The  structure  continued  to  take  load 
until  failure  occurred  thro\igh  crippling  of  several  posts  at  a  point 
8  percent  above  the  strain-reversal  buckling  load. 

The  results  of  the  test  showed  that  despite  differences  between  the 
test  specimen  and  the  idealized  structure  used  in  the  analysis,  the  mode 
of  buckling  predicted  .•/  the  theory  does  occur  at  the  predicted  load. 
However,  the  test  also  showed  that  significant  prebuckling  defomaticnis 
may  occur,  namely,  cross-sectional  deformations  due  to  downward  dishing 
of  the  covers  during  bending.  In  actual  wings,  such  cross-sectional 
deformations  could  very  likely  produce  vmdeslrable  aerodynamic  effects, 
and  it  would  therefore  be  necessary  to  use  suj^lementary  stiffening  to 
maintain  the  airfoil  shape.  Possibly  for  proportions  similar  to  those 
tested  a  proper  combination  of  i)osts  and  ribs  would  make  a  satisfactory 
design  both  from  the  viewpoint  of  deformations  and  strength.  Possibly 
for  some  other  proportions  posts  alone  would  be  adequate. 
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Figure. 2.-  Minimum  poet  axial  etlffness  required  for 
coDg>reaBlon  cover  to  buckle  with  tranaveree  nodes 
throu^  tbe  posts. 
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Figure  3<-  Buckling  phenomena  attalxiahle  for  box  beans  with 
rigid  or  effectively  rigid  posts. 


Figure  U.-  Box  beam  used  in  experimental  Investigation 
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Figure  Load- strain  curves  for  box  beam  stiffened  by  posts 
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STABILITY  OF  NULTIWEB  BOX  BEAMS 

By  Richard  A.  Pride,  Aldle  E.  Johnson,  Jr. , 
and  Roger  A.  Anderson 

Langley  Aeronautical  Laboratory 


The  thinness  of  vlngs  required  for  high-speed  aircraft  leads.  In 
SQiiie  cases,  to  the  use  of  nultlveb  construction,  (fig.  1).  The  skin  In 
the  central  portion  of  the  airfoil  section  may  be  slightly  ciirved  or 
flat  as  dictated  by  aerodynamics.  The  webs  may  be  Integrally  attached 
to  one  of  the  skin  surfaces  or  may  be  Joined  by  attctchment  angles  and 
rliwtlng. 

The  buckling  strength  of  the  resulting  multlveb  beam  has  been 
Investigated  theoretically  when  subjected  to  a  pure  bending  moment 
(reference  1).  It  has  been  found  that  for  a  beam  vlth  integral  Joints 
and  many  vebs,  beam  buckling  can  be  predicted  by  the  usual  plate  buck¬ 
ling  fonmila 


where  t|  Is  a  plasticity  correction  factor  and  the  value  of  k^  can 
be  obtained  from  figure  2.  In  this  figure  k^  Is  plotted  against  the 

ratio  of  web  depth  b^  to  the  veb  spacing  bp  for  several  ratios  of 

veb  thickness  to  cover-plate  thickness.  Buckling  coefficients  greater 
than  k  Indicate  beam  proportions  for  which  the  webs  provide  a  stabi¬ 
lizing  effect  to  the  covers.  Coefficients  below  k  indicate  beam  pro¬ 
portions  where  the  webs  have  a  destabllzlng  effect  on  the  cover. 

Crushing  stress  induced  in  the  webs  by  bending  of  the  covers  has  been 

neglected.  In  the  course  of  planning  experimental  research,  charts 
similar  to  figure  2  were  also  prepared  for  the  case  of  a  beam  with  only 
two  shear  webs  and  also  for  the  case  of  three  shear  webs,  whereas  this 
chart  Is  for  the  case  of  many  shear  webs.  When  the  k^  values  given  in 

figure  2  for  ^  «  1.0  for  multiwebs  are  conpared  with  the  corre¬ 
sponding  curves  for  the  two-  and  three-web  construction  (fig.  3)/  It 
can  be  seen  that  the  curves  do  not  greatly  differ.  Accordingly,  the 
simpler  two-web  section  was  selected  for  the  experimental  study. 

Results  of  pure  bending  teste  on  two-web  sections  are  shown  In 
figure  k.  The  experimental  buckling  stress  Is  plotted  against  calcu¬ 
lated  buckling  stress  and  the  closeness  of  the  points  to  a  k^®  line 


ccniFiDBrnAL 


56 


CONPIDaiTIAL 


indicates  the  adequacy  of  the  pure  bending  theory.  Both  elastic  and 
plastic  buckling  tests  are  included  in  this  data.  Various  airplane 
compcuiles  and  other  laboratories  have  conducted  tests  on  fabricated 
beams,  and  in  those  Instances  vhere  the  type  of  fabrication  approxi¬ 
mated  Integral  Joints  betveen  vebs  and  covers  (for  example,  refer¬ 
ences  2  and  3)»  the  agreement  with  theory  vas  satisfactory.  In  other 
tests  in  vhlch  beams  vere  fabricated  vlth  thin  vebs  and  heavy  covers 
Jointed  by  substantial  attachment  angles,  buckling  of  vebs  and  covers 
occurred  independently  and  the  present  integral-plate  theory  vas  not 
applicable.  Some  theoretical  and  experimental  research  is  being  done 
to  study  the  effect  of  attachment  angles  on  the  behavior  of  fabricated 
beams.  The  remainder  of  this  discussion  %rlll  be  confined,  however,  to 
beasts  in  vhlch  a  definite  Interaction  betveen  veb  and  cover  buckling 
exists. 

The  buckling  produced  by  a  spanvlse  variation  in  bending  moment 
vlll  now  be  considered.  Since  multlveb  beams  used  in  aircraft  wings 
are  subjected  to  distributed  loads,  the  webs  of  the  beam  carry  shear 
stress  and  the  bending  stress  varies  along  the  span.  A  buckling  ana¬ 
lysis  of  a  multlweb  beam  under  this  genial  kind  of  loading  has  not 
been  made.  However,  analyses  have  recently  been  made  of  the  webs  alone 
and  compression  covers  alone,  which  indicate  significant  qualitative 
effects. 

Figure  5  shows  the  theoretical  interaction  curve  between  pure 
bending  and  p\ire  shear  in  a  single  singly  suppoirted  web  plate.  The 
figure  is  plotted  in  terms  of  the  ratios  of  partlculcur  stresses  present 
under  combined  loads  to  the  critical  values  of  those  stresses  when  they 
are  acting  alone.  Gradient  in  the  bending  moment  has  been  neglected, 
but  the  effect  of  crushing  stress  of  the  type  produced  by  bending  of 
the  covers  has  been  included.  A  different  interaction  curve  is  shown 
for  each  of  several  valvies  of  the  crushing- stress  ratio.  The  flatness 
of  the  curves  near  their  upper  ends  indicates  that  the  ability  of  the 
webs  to  carry  bending  stress  is  not  reduced  substantially  until  the 
shear  stress  approaches  its  critical  value.  The  closeness  of  the 
curves  near  their  upper  ends  shows  that  for  moderate  values  of  shear 
stress  the  crushing  stress  has  little  effect  on  the  bending  strength  of 
the  veb.  It  is  only  vhen  the  shear  is  large  that  the  crushing  stress 
need  be  considered  as  a  contributing  factor. 

Figure  6  shovs  the  results  of  theoretical  calculations  (refer¬ 
ence  4)  for  buckling  in  compression  of  a  simply  supported  plate  repre¬ 
senting  a  portion  of  the  conpresslon  cover  bounded  by  two  vebs  and  tvo 
ribs.  The  effect  of  a  linear  stress  gradient  on  plates  of  three  length- 
to-vldth  ratios  is  shown.  Beneath  the  plates  are  the  corresponding 
nondlmenslonal  buckllng-stress  diagrams.  The  horizontal  line  is  the 
uniform  stress  required  to  buckle  the  plate  and  is  shown  as  100  percent 
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stress  at  tuckllng.  The  sloping  lines  represent  different  degrees  of 
gradient  stress  at  huckllng.  Here  It  Is  seen  that  the  stress  at  the 
■ore  highly  stressed  end  of  the  plate  Is  appreciably  higher  at  buckling 
than  the  unlfon  cospresslve  stress  required  to  buckle  the  plate.  This 
effect  Is  especially  pronounced  In  the  smaller  aspect  ratios. 

In  a  qualitative  vay,  the  theoretical  vork  Just  described  for  the 
Individual  components  Indicates  that  beam  buckling  strength  of  a  multl- 
veb  beam  as  predicted  by  pure  bending  theory  should  be  adversely 
affected  by  the  presence  of  shear  and  crushing  stresses  In  the  veba  but 
should  be  Increased  by  the  presence  of  stress  gradients  In  the  covers 
as  veil  as  the  vebs.  In  order  to  study  the  over-all  effect  of  web  shear 
and  moment  gradient  on  a  beam,  HACA  has  proceeded  vlth  an  experimental 
program. 

Figure  7  shows  the  type  of  tqst  being  conducted.  The  tvo-veb 
section  Is  supported  as  a  cantilever  beam  vlth  a  concentrated  shear  load 
and  end  moment  applied  at  one  end  of  the  beam.  This  produces  a  constant 
shear  diagram  and  a  moment  diagram  which  varies  linearly  from  a  miniiBuiB 
at  the  outboard  end  to  a  wiAviimmi  at  the  Inboard  end  of  the  test  section. 
By  use  of  the  combined  load  testing  machine  of  the  Langley  structxires 
research  laboratory,  the  values  of  the  minimum  moment  and  the 
moment  and  the  shear  load  at  buckling  can  be  changed  In  the  test  section. 
Ibe  test  section  amy  be  thought  of  as  representing  any  one  bay  between 
bulkheads  of  a  wing. 

The  results  of  tests  are  presented  In  figure  8,  as  a  series  of 
moment  diagrams  at  buckling  for  a  beam  subjected  to  various  loadings. 

The  length  between  bulkheads  of  this  particular  beam  Is  9.6  times  the 
spacing  of  shear  webs.  The  horizontal  line  is  the  pure  bending  mosient 
required  to  buckle  the  beam  and  Is  shown  as  100  percent  pure  bending 
strength.  The  lines  of  various  slopes  represent  moment  diagrams 
resulting  from  various  amounts  of  shear  loading.  The  maximum  moment 
developed  at  buckling  was  greater  than  100  percent,  and  the  length  over 
which  the  unlfonn  buckling  stress  was  exceeded  was  about  twice  the 
’  spacing  of  shear  vebs.  For  this  particular  beam  this  result  Is  not 
surprising,  since  for  the  greatest  shear  loading  the  web  shear  stress 
was  only  about  6  percent  of  Its  critical  value.  Also,  the  crushing 
stress  was  only  about  1  percent  of  its  critical  value.  These  small 
values  result  In  an  extreisely  small  reduction  of  pure  bending  strength 
for  a  web  plate  (see  fig.  ^).  The  corresponding  reduction  for  the  beam 
as  a  whole  Is  negligible. 

NACA  Is  proceeding  vlth  tests  on  beams  vlth  vebs  thinner  than  the 
cover  so  that  the  web  shear  stress  will  be  larger  In  relation  to  its 
critical  value  than  In  the  test  results  presented  here.  On  the  basis 
of  the  tests  that  have  been  completed.  It  can  be  concluded  that  for 
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long  teams  under  gradient  tending  In  idilch  the  veb  shear  and  crushing 
stresses  are  not  large  In  coaqparlson  to  their  critical  values,  the  pure 
tending  strength  can  be  exceeded  over  a  length  equal  to  one  to  tvo 
times  the  spacing  of  shear  vebs  at  the  oiore  highly  stressed  end  of  the 
beam. 
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Figure  C(niblDatlo&8  of  shear,  bending,  and  crushing  stress 
vhlch  produce  buckling  of  an  Infinitely  long  flat  plate  with 
sliqply  supported  edges. 
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Figure  6.-  Theoretical  stress  dlagraas  at  buckling  for  slaqply 
supported  plates  under  a  eo8Q>ressive  stress  that  varies 
linearly  in  tbs  direction  of  loading. 
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Figure  7.-  Loading,  shear,  and  nonent  diagrams  for  test  specimen. 


Figure  6.-  Experimental  moment  diagrams  at  buckling  for  test 
specimen  under  gradient  bending* 
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DESICaV  FOR  STIFFNESS  AS  WELL  AS  STRENGTH 
By  Norris  F.  Dow 
Langley  Aeronautical  Laboratory 


INTRODUCTION 


While  there  have  been  a  number  of  studies  of  minimum  weight  design 
of  aircraft  structures,  most  of  them  have  been  concerned  primarily  with 
the  strength -to -weight  ratio  of  the  structure  rather  than  the  stlffhess- 
to-welght  ratio.  The  reason  that  the  major  emphasis  has  been  upon 
strength  rather  than  stlffhess  Is  that  until  recent  years  a  structure 
of  adequate  strength  would  usually  achieve  adequate  stiffness  automat¬ 
ically.  Now  stlffhess  Is  becoming  more  of  a  problem,  partly  because 
higher  speeds  are  Imposing  more  severe  stiffness  requirements  and  the 
accompanying  thinner  wings  are  making  these  requirements  harder  to  meet, 
and  partly  because  Increases  have  been  made  In  the  strength-to-welght 
ratio  without  corresponding  Increases  in  the  stiffness -to -weight  ratio, 
so  that  the  more  efficient  structure  from  a  strength  standpoint  has 
become  less  likely  to  be  adequately  stiff. 

The  present  paper  Is  a  study  to  determine  what  are  the  possibilities 
of  increases  in  stiffness -to -weight  ratio.  This  study  Is  made  along 
two  lines  to  determine  the  possibilities  associated  with 

(1)  The  use  of  the  most  effective  distribution  of  a  given  material 

(2)  The  use  of  the  most  effective  material 

A  few  detailed  studies  of  both  stiffness  and  strength  are  made  along 
^th  these  lines,  and  current  and  possible  future  developments  are 
discussed . 


most  effective  DISmiBUTION  OF  MATERIAL 


First  let  us  consider  the  possibilities  associated  with  the  most 
effective  distribution  of  material.  Ihe  principles  Involved  In  the 
distribution  of  material  for  maximum  stiffness  are  extremely  elementary, 
consisting  basically  of  putting  the  material  where  It  will  carry  the 
load  most  directly.  If  you  have  a  simple  tension  load,  you  don't  carry 
It  with  a  helical  spring;  you  use  a  straight  rod.  If  you  have  a  bending 
load,  you  don't  use  a  plus-sign  section;  you  use  an  I-beam.  (See  fig.  1. ) 
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By  the  same  token,  if  you  are  primarily  concerned  with  torsional  stiff¬ 
ness  of  the  wing,  you  would  prefer  not  to  put  a  hlg  percentage  Of 
material  in  longitudinal  stiffeners  but  to  put  as  much  as  possible  of 
the  material  in  the  skin  where  it  will  work  for  you,  as  suggested 
schematically  in  figure  2. 

In  more  technical  language,  the  principles  for  obtaining  maximum 
stlffhess  can  be  stated  thus: 

(1)  Put  the  material  at  the  point  of  meixlmum  stress,  in  the  direc¬ 
tion  of  maximum  stress. 

(2)  Keep  the  distance  over  which  the  high  stress  acts  to  a  mlnimiun. 

(3)  Avoid  local  or  general  bending  or  twisting  distortions  which 
reduce  the  effectiveness  of  the  material. 

(4)  Avoid  slippages  in  Joints  by  making  the  structure  as  integral 
as  possible. 

Because  of  the  desirability  of  making  everything  as  Integral  as 
possible  and  because  of  general  Interest  in  Integral  stiffening,  a 
detailed  study  has  been  undertaken  of  the  characteristics  of  sheet  with 
various  configurations  of  integral  stiffening.  Integrally  stiffened 
sheet  has  been  proposed  in  a  number  of  configurations  (references  1  to  4), 
varying  from  fairly  simple  ribbing  to  complex  grldworks  tailored  for' 
specific  designs.  In  order  to  make  our  study  cover  as  many  of  these 
configurations  as  could  reasonably  be  done,  we  considered  first  simple 
longitudinal  or  transverse  ribbing  (see  fig.  3) t  and  then  the- case  for 
which  the  ribs  are  skewed  at  an  angle  9  to  the  longitudinal  axis.  Ihe 
resulting  waffle-like  plates  have  densities  and  properties  that  can  be 
varied  by  varying  size,  spacing,  and  orientation  of  ribs.  The  present 
paper  is  concerned  primarily  with  the  stlffhess  properties  of  waffle-  i 
like  plates,  si>eclflcally  their  axial  and  shear  stiffnesses,  which  are 
related  to  the  stlffhesses  which  waffles  can  provide  to  resist  bending 
or  twisting  of  the  structure  as  a  whole,  and  the  bending  and  twisting 
stiffnesses  of  the  waffles  themselves,  which  are  related  to  their  own 
resistance  to  local  bending  or  buckling  distortions.  After  all,  the 
main  reason  for  having  stiffening  of  any  kind  is  to  Increase  the  resist¬ 
ance  to  local  deformations.  If  it  were  not  for  such  local  deformations 
or  buckling,  slsqple  solid  sheet  could  be  used,  and  such  solid  sheet 
would  always  pzt>vlde  a  structure  of  maximum  stiffness  for  a  given  weight 
of  a  material.  The  problem  is  to  increase  the  buckling  resistance  of 
sheet  with  the  least  possible  loss  in  its  compressive  and  shear  stiff¬ 
ness  properties. 
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Let  us  examine  first  the  compressive  and  shear  stiffness  properties 
of  a  typical  vaffle.  With  the  vaffle  as  the  upper  surface  of  the  wing, 
the  con^ressive  stiffness  measures  the  stlffhess  provided  against 
handing  of  the  wing,  and  the  shear  stiffness  measures  the  stiffness 
provided  against  twisting  of  the  wing.  Curves  of  compression  and  shear 
stiffnesses  for  a  typical  waffle  proportion  are  plotted  in  figure  4. 

In  this  figure  the  ordinate  is  the  ratio  of  stiffnesses  of  waffle  and 
solid  sheet  of  equal  weight,  and  the  abscissa  Is  the  angle  d  at  which 
the  ribs  are  skewed.  The  curves  were  calculated  assvunlng  pin  Joints  at 
the  comers  of  the  grids.  With  such  pin  Joints  the  compressive  stiff¬ 
ness  decreases  to  a  minimum  with  the  ribs  at  approximately  60°,  and 
the  shear  stiffness  increases  to  a  maximum  with  the  ribs  at  45° • 
Accordingly,  from  the  standpoint  of  stiffness  of  the  wing  as  a  whole, 
the  choice  of  the  best  angle  for  the  stiffening  elements  depends  upon 
the  relative  importance  of  the  bending  stiffness  and  the  torsional 
stiffness  of  the  particular  wing  in  question. 

As  might  be  expected,  skewing  the  stiffening  elements  also  affects 
the  stiffness  of  the  waffle  to  resist  local  deformation  or  buckling. 
Evidently,  the  choice  of  the  best  angle  for  the  stiffening  for  a  partic¬ 
ular  application  can  hardly  be  made  on  the  basis  of  axial  or  shear 
stiffness  properties  alone.  Calculations,  similar  to  those  for  the 
axial  and  shear  stiffnesses,  of  the  bending  stiffness  and  the  twisting 
stiffhess  of  the  waffle  show  that  these  properties  vary  even  more 
substantially  than  the  conpression  and  shear  properties  as  the  emgle 
of  ribbing  is  varied.  The  curves  of  bending  and  twisting  stiffhesses 
are  plotted  in  figure  5  as  the  ratio  of  stiffnesses  of  waffle  and  equal 
weight  solid  against  the  angle  of  ribbing.  While  the  compression  and 
shear  stiffnesses  were  never  greater  than  those  of  equal  weight  solid 
sheet,  the  bending  euid  twisting  stiffnesses  can  be  substantially  greater 
and,  hence,  the  buckling  strength  can  be  increased. 

In  order  to  determine  approximately  what  are  the  effects  of  skewed 
stiffening  on  the  buckling  strength,  buckling  loads  were  estimated  from 
the  bending  and  twisting  stiffnesses  using  orthotropic  plate  theory. 

The  resulting  values  for  the  buckling  loads  showed  that  for  a  given 
waffle  proportion  in  the  elastic  range  the  critical  compressive  load 
should  increase  as  the  ribs  are  skewed  because  of  the  profound  effect  on 
buckling  of  the  twisting  stiffness  which  is  a  maximum  at  45°.  Although 
this  was  only  a  preliminary  analysis  it  was  qualitatively  confirmed  by 
elastic  buckling  tests  of  long  square  tubes  made  up  from  waffles .  In  fig¬ 
ure  6  are  plotted  the  results  of  the  buckling  tests  on  waffles  having 
ribs  at  0°,  30°,  and  45°.  In  each  case  the  circle  represents  the 
measured  critical  load  divided  by  the  critical  load  for  a  solid  sheet 
of  equal  width  and  weight.  The  results  confirm  that  the  stifi'ening 
elements  can  be  skewed  to  increase  the  shear  stiffness,  and  at  the  same 
time  the  elastic  compressive  buckling  strength  is  actually  increased. 
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If  the  waffle  Is  compared  with  conventionally  stiffened  (that  Is, 
riveted  up)  construction.  It  appears  to  have  several  favorable  features 
from  the  stiffness  standpoint  quite  apart  from  those  associated  with 
favorable  angular  orlentr.  tlon .  First,  because  the  waffle  cem  readily 
be  provided  with  high  twisting  as  well  as  longitudinal  and  transverse 
bending  stiffnesses.  It  can  be  expected  to  resist  buckling  with  less 
over -all  depth  of  section  than  the  conventional  panel.  Accordlnj^ly, 

It  can  be  effectively  placed  out  closer  to  outside  fibers  to  contribute 
greater  bending  or  twisting  stiffnesses  to  the  structure  as  a  whole. 

Second,  proper  proportioning  of  skin  and  grid  permits  Ineffectivenesses 
due  to  local  buckling  distortions  of  the  skin  to  be  readily  avoided. 

And  third,  any  losses  In  stiffness  due  to  loosenesses  In  riveted  Joints 
between  skin  and  stringers  are  obviously  completely  eliminated. 

!nie  foregoing  discussion  of  waffle-type  stiffening  has  been  con¬ 
cerned  primarily  with  the  stiffness  properties,  euid  the  evaluation  of 
stiffness  properties  are  only  one  small  phase  of  the  structural  evalua¬ 
tion  of  such  integrally  stiffened  construction.  Studies  of  other 
properties  such  as  ultimate  strength  characteristics  are  proceeding, 
and  Langley  Is  negotiating  with  the  Wyman-Ciordon  Company  In  order  to  obtain 
some  waffles  fabricated  as  press -forgings  to  aid  In  these  studies.  The 
test  specimens  which  were  made  to  check  the  buckling  strengths  were 
castings  and  hence  are  of  no  value  for  other  than  such  tests  entirely 
In  the  elastic  stress  range. 


liK)ST  EFFECTIVE  MATERIAL 


Let  us  now  turn  to  the  consideration  of  the  possibilities  asso¬ 
ciated  with  the  use  of  the  most  effective  material.  Perhaps  the  first 
Idea  that  comes  to  the  mind  of  the  aircraft  designer  confronted  with 
the  problem  of  making  the  structure  stiffer  is  to  use  steel  in  place 

of  aluminum.  Steel  has  a  Young's  modulus  2^  to  3  times  that  of  aluminum 

and  its  smller  volume  for  the  same  wei^t  permits  it  to  be  crowded 
out  closer  to  outside  fibers.  This  smaller  volume,  however,  is  usually 
the  stumbling  block,  for  associated  with  it  are  the  well-known  buckling 
difficulties  always  associated  with  thin  sheets,  and  buckling  is  of 
necessity  accompanied  by  loss  of  effective  widths  of  skin  and  hence 
losses  in  both  stiffhess  and  strength.  Buckling  difficulties  can  of 
course  be  reduced  by  the  standard  procedures  of  adding  stiffeners,  or 
possibly  integrally  stiffened  steel  plates  similar  to  the  various  types 
Just  described  could  be  fabricated. 

Another  possibility  to  be  considered  is  a  solid  steel-faced 
meigne slum-core  sandwich.  Let  us  examine  the  stiffness  properties  of 
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such  a  sandwich  In  the  same  way  that  we  did  for  the  waffle.  In  this 
case  the  variable  Is  the  relative  thickness  of  faces  and  core.  Indicated 
plctorlally  In  figure  7  by  the  curve  and  quantitatively  by  the  value 
of  the  abscissa.  In  this  figure  Is  plotted  the  bending  stlffhess  of 
the  sandwich  conpared  to  solid  aluminum  of  eqxial  weight.  One  end  of 
the  curve  represents  zero  core,  that  Is  solid  steel;  the  other  end 
represents  solid  magnesium.  At  the  same  time  that  the  relative  propor¬ 
tions  of  face  and  core  are  changed,  the  over-«ll  thickness  of  the  semd- 
wlch  Is  changed,  as  represented  by  the  width  of  the  curve.  In  order  to 
keep  the  weight  constant.  For  the  point  on  the  curve  corresponding  to 
the  proportions  for  which  the  sandwich  has  the  same  over -all  thickness 
as  well  as  the  same  weight  as  alvimlnum,  the  sandwich  has  approximately 
40  percent  greater  bending  stiffness  than  aluminum,  or  If  the  faces 
are  made  slightly  thicker  and  the  core  correspondingly  thinner  so  that 
the  weight  Is  still  the  same,  the  sandwich  can  be  thinner  and  yet  still 
somewhat  stlffer  In  bending  than  aluminum. 

This  curve  represents  only  the  elastic  bending  stlffhess  and 
correspondingly  the  resistance  to  local  deformation  or  buckling  In  the 
elastic  stress  range.  The  properties  In  the  plastic  range  depend  upon 
the  stress -strain  properties  of  the  steel  and  magnesium  alloy  used, 
and  both  very  high-strength  steel  faces  and  a  high-strength  magnesium- 
alloy  core  are  needed  If  the  sandwich  Is  to  compete  with  a  hl^-strength 
aluminum  alloy  such  as  7^S-T6  at  high  stresses  beyond  the  elastic  range. 
For  the  case  of  75S-T6  plate  and  a  sandwich  of  equal  weight  and  thick¬ 
ness,  for  exanple,  these  curves  of  figure  8  show  the  relative  stress - 
strain  characteristics  when  the  sandwich  Is  made  with  4340  steel  faces 
and  ZK6QA  magnesium-alloy  core.  While  a  comparison  on  the  basis  of 
efficiency  to  resist  plate  buckling  would  be  much  more  favorable  to  the 
sandwich  than  this,  still  these  curves  do  Indicate  the  necessity  of 
using  the  best  possible  materials  If  the  sandwich  Is  to  compete  with 
7^S  In  the  high  stress  region.  Except  for  this  region  the  sandwich  cein 
provide,  with  less  thickness,  resistance  to  buckling  equal  to  that  for 
aluminum  of  equal  weight .  Accordingly,  for  thin  wings  for  which  stiff¬ 
ness  Is  critical  a  steel-clad  magnesium  skin  can  provide  greater  bending 
or  torsional  stlffhess  to  the  wing  as  a  whole  than  aluminum  skin  of  the 
same  wel|^t  because  the  sandwich  skin  can  be  made  thinner  without 
reducing  the  reslstemce  to  local  deformation  or  buckling  and  because 
being  thinner  It  can  be  placed  effectively  nearer  the  extreme  fibers 
of  the  wing. 

At  the  present  time  the  General  Plate  Company  Is  attenqptlng  to 
fabricate  for  us  steel-clad  magnesium  sandwiches  for  structural  evalua¬ 
tion.  General  Plate  has  successfully  bonded  a  number  of  dissimilar 
materials  (see  reference  auid  has  had  acme  success  bonding  steel  on 
magnesium  on  small  pieces,  but  as  yet  nothing  of  any  size  has  been 
produced. 
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One  feature  to  remember  particularly  about  the  curves  for  the 
steel-clad  magnesium  sandwich  Is  that  while  the  sandwich  can  have,  for 
the  same  weight  and  thickness,  greater  elastic  bending  stlffhess  than 
aluminum,  or  for  the  same  bending  stiffness  and  weight  It  can  have  less 
thickness  than  aluminum  and  hence  can  provide  greater  bending  or  twisting 
stiffness  to  a  thin  structure  as  a  whole,  nevertheless,  even  In  the 
elastic  range.  It  has  no  greater  axial  stlffhess  than  aluminum;  that  Is, 
the  Initial  slopes  of  the  stress-strain  curves  for  the  sandwich  and  for 
equal  weight  aluminum  are  the  sasie.  Improve  the  axial  stiffness 
characteristics  (and  correspondingly  the  shear  stiffness  characteristics) 
a  material  with  a  higher  ratio  of  modulus  of  elasticity  to  density  than 
magnesium,  aluminum,  or  steel  Is  needed. 

The  fact  that  such  metals  as  magnesium,  aluminum, 

Iron  all  have  approximately  the  same  modulus -to -density  ratio,  as  noted 
In  the  following  table,  has  given  rise  to  the  Isqpresslon  that  there  Is 
little  hope  of  any  Improvement  In  this  ratio: 


Material 

Modulus 
ClO^  psl) 

Modulus 

Density 

(lO^  In.') 

Ifeigneslum 

6.5 

103 

Aluminum 

10 

102 

Titanium 

16 

100 

Iron 

28 

100 

Higher  values  of  the  modulus -to •density  ratio,  however,  are  not  contrary 
to  nature.  A  number  of  materials  having  high  modulus -to -density  ratios 
are  listed  In  the  following  table: 


Modulus 

Modulus 

Material 

(lO^  psl) 

Density 
(iC^  in.') 

Aluminum 

10 

102 

Chromium 

•  26 

108 

Molybdenum 

46 

125 

Tungsten  carbide 

96 

177 

Titanium  carbide 

55 

260 

Berylllvun 

42 

640 
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Chroalum  Is  a  little  better  than  aluatlnum;  nolybdenw  Is  even  better. 

In  fact  Bolybdenum  mlj^t  veil  be  suitable  for  some  highly  specialized 
applications  for  which  high  stlffhess  Is  a  primary  requisite.  If  a 
face  material  similar  to  molybdenum  were  available  for  cladding  magne¬ 
sium,  for  example.  It  vould  make  a  very  Interesting  sandwich.  Molyb¬ 
denum  Itself  Is  probably  too  scarce  to  make  a  molybdenum-clad  magnesium- 
core  sandwich  useful  for  other  than  limited  applications.  For  possibly 
even  more  limited  applications  there  are  the  carbides.  For  exeuiqple, 
there  Is  tungsten  carbide  which  has  a  Young's  modulus  of  96,000,000  psl 
and  a  modulus -to-denslty  ratio  of  ITT,  and  there  Is  titanium  carbide 
with  a  modulus  of  ^^,000,000  psl  and  a  modulus -to-denslty  ratio  which 
Is  nearly  three  times  that  of  aluminum.  Finally  there  Is  beryllium. 
Included  here  primarily  for  completeness  to  show  that  a  high  modulus - 
to-denslty  ratio  Is  not  contrary  to  nature,  for  even  If  berylllvun 
should  become  available  In  quantity  Its  poisonous  characteristics  might 
preclude  Its  use. 

The  carbides  are  becoming  available  in  somewhat  bigger  pieces. 
Possibly  carbides  should  be  considered  as  presently  applicable  to  such 
things  as  thin  wings  for  very  small  missiles  or  to  small  wind-tunnel 
models  for  which  stiffness  Is  the  primary  requisite,  but  even  for  such 
applications  only  If  adequate  allowance  for  the  limitations  of  the 
carbides  such  as  the  fact  that  they  are  very  brittle  can  be  made. 

Further  study  of  their  limitations  is  proceeding,  amd  the  Kennametal 
Corporation  is  under  contract  to  supply  us  with  carbide  test  specimens 
to  aid  In  this  study. 

Despite  such  limitations  as  scarcity,  brittleness,  or  toxicity, 
molybdenum,  the  carbides,  and  beryllium  do  show  that  a  high  modulus - 
to-denslty  ratio  Is  not  contrary  to  nature.  Since  a  high  modulus -to- 
denslty  ratio  Is  not  contrary  to  nature.  Is  there  any  possibility  of 
liiq>rovement  Ip  the  modulus -to -density  ratio  of  common  materials?  Most 
authorities  are  of  the  opinion  that  there  Is  no  possibility  of  Improve¬ 
ment  In  either  modulus  or  modulus -to -density  ratio.  There  are,  however, 
a  few  dissenters.  Ibe  chief  dissenters  are  Koster  in  Germany  (refer¬ 
ence  6)  and  Dudzlnskl  In  Great  Britain  (references  T  and  8),  but  there 
have  also  been  a  few  Investigators  in  this  country  (references  9  to  11) 
who  have  suggested  that  Young's  modulus  Is  susceptible  to  memlpulatlon 
by  alloying  In  much  the  same  way  that  other  mechanical  properties  of 
metals  can  be  manipulated.  Dudzlnskl  has  si^ggested  that  the  modulus - 
to-denslty  ratio  as  well  can  be  lnq>roved. 

Because  this  statement  Is  a  somewhat  controversial  one,  let  us  review 
some  of  the  data  on  which  It  Is  based.  First  consider  the  data  based  on 
Koster 's  work,  presented  by  Zener  as  evidence  that  small  amounts  of 
alloying  elements  which  form  solid  solutions  may  be  expected  to  decrease 
the  modulus.  These  data  are  reproduced  In  figure  9  which  shows  the 
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effect  <Hi  Young's  aodulus  of  solid  solutions  of  aluBlnua,  aanganese, 
beryllluB,  and  aagneslua  In  copper.  lot Ice  particularly  that  both 
■anganese  and  berylllua  decrease  the  aodulus  when  alloyed  vlth 
capper,  even  thouc^  their  aodull  Individually  -  given  by  the  nuad)ers 
In  brackets  -  are  both  greater  than  that  of  copper.  The  fact  that  the 
aodulus  of  the  alloy  Is  not  simply  the  weighted  average  of  the  moduli 
of  the  elements  suggests  the  possibility  that  for  some  cOid)lnatlons 
iaproveaents  might  be  realized.  Actually,  Just  a  little  farther  along 
the  curve  for  aluminum  alloyed  vlth  copper,  beyond  the  point  at  which 
simple  solid  solutions  are  formed,  the  modulus  Increases  to  a  value 
substantially  higher  than  that  for  either  aluminum  or  copper  alone,  as 
shown  in  figure  10.  This  plot  Is  also  based  on  KSster's  work.  Here 
again  Is  shown  the  effect  on  Young's  modulus,  this  time  for  various 
percentages  of  aluminum  alloyed  with  copper.  Finally,  let  us  look  at 
possible  variations  In  the  modulus -to-denslty  ratio  (fig.  11).  Here 
the  data  are  derived  from  Dudzlnski's  work  In  which  he  shows  how  the 
modulus -to-denslty  ratio  of  aluminum  may  be  changed  by  adding  alloying 
elements.  Dudzlnskl  has  gone  so  far  as  to  state  that  he  believes  a 
practical  aluminum-alloy  with  a  Young's  modulus  of  12,000,000  psl  can 
be  developed.  This  prediction,  depending  as  it  does  upon  the  Interpre¬ 
tation  of  the  word  'practical”,  does  not  seem  so  extravagant  when  it 
Is  considered  that  the  Aluminum  Cosqpany  of  America  lists  its  alloy  32S 
as  having  a  Young's  modulus  of  11,500,000  psl. 

Despite  data  such  as  these,  the  prevailing  Impression  seems  to  be 
that  there  Is  no  hope  of  Improvement  In  the  modulus -to -density  ratio. 

In  fact  the  possibility  has  been  suggested  that  KSster's  and  Dudzlnski's 
data  may  be  drastically  in  error. 

In  an  effort  to  help  clarify  this  situation,  Langley  requested  the 
Aluminum  Company  of  America  to  make  some  aluminum-alloy  sheet  having 
a  Young's  modulus  of  12,000,000  psl  and  a  density  approximately  the 
same  as  that  of  pure  aluminum.  The  Aluminum  Company  of  America  assured 
us  that  the  sheet  could  be  issde  but  that  It  would  be  a  far  cry  from 
the  practical  alloy  that  Dudzlnskl  was  trying  to  develop.  That  Is 
unquestionably  how  It  turned  out  because  among  other  things,  the  imaterlal 
Is  very  brittle.  Ibdenlable,  however,  as  the  stress -strain  curves 
given  In  figure  12  show,  by  adding  alloying  eleawnts  to  pure  aluminum 
they  have  both 

(1)  Increased  the  mechanical  properties 

(2)  Increased  the  Young's  aK>dulus 

niey  did  not  quite  achieve  the  12,000,000  they  were  shooting  for 
In  the  vlth^rain  direction  but  they  exceeded  It  In  the  cross -grain 
direction.  However,  since  the  average  elongation  In  the  vlth-graln 
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direction  was  between  2  and  3  percent  over  a  2-inch-eage  length,  at 
the  sane  tine  that  they  were  raising  the  modulus,  they  were  also 
increasing  the  brittleness . 

At  present  it  appears  that  most  hl£^  modulus -to-denslty-ratio 
materials  are  brittle.  The  problem  that  evidently  must  be  solved  if 
higher  modulus-to-denslty -ratio  structural  materials  are  to  be  achieved 
appears  to  be  that  of  providing  ductility.  Whether  this  goal  can  be 
achieved  is  still  a  matter  for  speculation. 


CONCLUDING  RafABKS 


Thken  as  a  whole,  the  problem  of  designing  aircraft  structures  for 
stiffness  as  well  as  strength  should  be  considered  to  be  in  its  Infancy, 
and  so  considered  It  seems  to  be  a  healthy  Infant.  It  is  a  problem 
that  promises  to  grow  In  Importance,  and  its  future  certainly  does  not 
appear  hopeless.  Integrally  stiffened  construction,  new  materials,  and 
new  combinations  of  materials  all  offer  possibilities  6md  certainly  in 
time  at  least  some  of  the  possibilities  will  be  realized. 
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Figure  S**  Configurations  of  Integrally  atlffened  sbMt  considered. 


Figure  4.-  CoHpressl've  and  shear  stiffnesses  of  waffles. 
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Figure  Bendlog  and  twisting  stiffnesses  of  waffles 
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Flgiire  6.-  Saqperlaestal  Iraclillng  loads  for  waffles 
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BENDING  STIFFNESS 
FOR  EQUAL  WEIGHT 
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THICKNESS  RATIO, 

Figure  7«-  Bending  etlffneeees  of  steel-clad  magnesium  sandwiches 
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BINARY  ALLOYS  OF  COPPgR 


ATOMIC  PERCENT  ALLOY  ELEMENT 

Figure  9**  Effect  on  modulus  of  adding  alloying  elements  In 
solid  solution  with  copper. 


Figure  10. •  Effect  modulus  of  alloying  aluminum  vlth  copper. 
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BINARY  ALLOYS  OF  ALUMINUM 


Figure  11.-  Effect  on  noduluB-to-denslty  ratio  of  adding 
alloying  elements  to  aluminum. 


Figure  12.-  Stress-stnln  curves  for  special  hl^  modulus 

aluminum  alloy. 
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THEORIES  OF  PIASTICITY 
By  S.  B.  Batdorf  and  Bernard  Budlansky 
Langley  Aeronautical  Laboratory 


Before  the  occurrence  of  structural  failure,  plastic  flow  usually 
derelops.  This  plastic  flow  reduces  the  stiffness  of  the  material, 
changes  the  stress  distribution,  gives  rise  to  residual  stresses,  and 
produces  permanent  distortions.  In  addition,  it  plays  a  role  in  the 
fracture  of  ductile  materials  and  probably  is  iiwolved  also  in  cieep 
and  fatigue.  A  portion  of  our  research  effort  has  therefore  been 
devoted  to  the  development  of  an  adequate  understanding  of  plasticity. 
Although  this  goal  has  not  yet  been  achieved,  it  is  felt  that  substan¬ 
tial  progress  has  been  made,  and  this  paper  will  constitute  a  report  on 
our  activities  in  this  field.  As  such  it  will  be  primarily  an  account 
of  a  seeurch  for  an  acceptable  theory.  Practical  applications  will  be 
the  concern  of  the  next  paper. 

The  general  character  of  the  stress-strain  behavior  of  materials 
in  simple  tension  or  compression  has  been  known  for  a  long  time  and 
this  has  permitted  adequate  solutions  to  be  found  for  the  stress  dis¬ 
tribution,  buckling  loads,  and  meocimum  strength  of  beams  and  columns. 
The  solution  of  the  corresponding  problems  for  plates  and  shells  such 
as  go  to  make  up  the  bulk  of  the  airframe,  however,  requires  a  knowl¬ 
edge  of  the  plastic  strain  laws  for  arbitrarily  varying  combined- stress 
loading  conditions,  and  Just  what  form  these  laws  should  take  has  not 
as  yet  been  established. 

Figure  1  shows  a  plot  of  a  typical  stress-strain  curve  for  a 
strain  hardening  material  in  simple  tension.  There  are  two  equations 
in  general  use  for  describing  the  stress-strain  relations  beyond  the 
elastic  range  of  the  material.  One  of  these,  a  =  EgS,  relates  the 
total  stress  to  strain  by  means  of  a  induced  modulus  called  the  secant 
modulus  which  is  defined  in  the  figure.  The  other,  da  =  I^dc,  relates 
Increment  of  stress  to  Increment  of  strain  by  using  the  slope  of  the 
curve  or  the  tangent  modulus.  When  these  laws  are  generalized  for 
application  to  combined-stress  loading  conditions,  they  retain  their 
character  of  having  either  finite  or  incremental  form.  The  general¬ 
ization  of  the  finite  one  has  been  guided  in  part  by  the  form  of  the 
elastic  equations  for  three-dimensional  states  of  stress,  while  the 
generalization  of  the  Incremental  one  has  been  influenced  by  the  equa¬ 
tions  for  fluid  flow.  Althou^  the  two  types  of  law  are  equivalent  in 
the  case  of  simple  tension,  in  their  general  fonnulations  they  are  not. 
One  of  the  most  conspicuous  differences  is  that  finite  laws  imply  that 
stress  and  strain  are  uniquely  associated,  whereas  according  to  incre¬ 
mental  theories  the  relation  between  stress  and  strain  depends  upon  the 
loading  path,  or  stress  history. 
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Analyses  based  on  plastic  stress-strain  lavs  are  much  more  compli¬ 
cated  than  those  based  on  elastic  lavs  because,  although  the  equations 
of  elasticity  are  linear,  those  of  plasticity  are  nonlinear.  Of  the 
tvo  tyi)es  of  plastic  stress-strain  lav,  hovever,  the  finite  lav  is  much 
easier  to  use  because  vlth  it  the  stress  distribution  can  be  obtained 
directly  in  terms  of  the  final  load.  The  incremental  lav,  on  the  other 
hand,  requires  that  once  plastic  action  has  started  the  stress  distri¬ 
bution  must  be  found  for  each  increment  of  load,  this  process  being  con¬ 
tinued  until  the  final  load  is  reached.  One  vould  therefore  prefer  to 
•  use  the  finite  type  of  stress-strain  lav.  More  important,  hovever,  is 
the  question  of  vhich  is  right.  Both  types  of  equation  can  be  vritten 
vith  adjustable  parameters  vhich  can  be  so  chosen  as  to  agree  vlth 
experimental  results  for  proportional  loading  -  that  is,  loading  in 
vhich  all  stress  components  increase  at  the  same  rate.  When  this  is 
done,  hovever,  the  tvo  types  of  lav  vlll  disagree  as  to  the  stress- 
strain  relationships  for  all  other  types  of  loading  and  the  disagree¬ 
ment  is  sharpest,  generally  speaking,  vhen  the  loading  deviates  the 
most  strongly  from  proportional  loading.  One  of  the  cases  in  vhich  the 
loading  departs  strongly  from  proportionality  is  that  of  buckling.  So 
let  us  nov  see  vhat  kind  of  results  are  obtained  here  from  the  tvo  types 
of  theory. 

During  the  var,  an  extensive  series  of  tests  on  plastic  buckling 
of  plates  vas  conducted  at  the  Langley  Laboratory.  On  the  basis  of 
these  and  other  data,  Gerard,  then  of  the  Republic  Aviation  Company, 
proposed  the  secant  modulus  as  the  effective  modulus  for  plate  buckling 
(reference  1).  (The  same  proposal  vas  made  lndei>endently  by  McDonald 
of  Dov  Chemical  Company.)  Shortly  after  this,  Handelman  and  Prager  of 
Brovn  University  gave  a  theoretical  solution  to  the  problem  of  plastic 
buckling  of  plates  on  the  basis  of  the  incremental  theory  of  plasticity 
(reference  2)  and  Stovell  of  Langley  solved  it  on  the  basis  of  the 
finite  theory  of  plasticity  (reference  3).  (Incidentally,  both  solu¬ 
tions  had  been  obtained  by  Bijlaard  in  1938,  but  his  vork  vas  not  gener¬ 
ally  knovn  in  this  country;  a  reviev  of  his  vork  appears  in  reference  U. ) 
In  figure  2  the  results  of  the  tvo  theories  are  compared  vith  experi¬ 
mental  data  on  the  buckling  of  long  hinged  flanges.  In  this  figure 
buckling  stress  is  plotted  against  the  elastic  buckling  strain.  The 
ratio  of  the  heights  of  the  buckling  curve  and  the  elastic  line  ma^  be 
regarded  as  a  knock-dovn  factor  to  be  used  above  the  elastic  limit.  At 
higher  stresses  the  finite  theory  says  there  is  a  knock-dovn  factor  and 
the  incremental  theory  says  there  is  not.  In  fact,  in  this  particular 
case  the  finite  theory  predicts  a  curve  coinciding  vith  the  stress- 
strain  curve  for  the  material,  vhich  meems  the  secant  modulus  aiiplles. 
Test  data  obtained  from  extruded  cruciform  section  columns  that  buckled 
by  tvlstlng  are  in  good  agreement  vlth  the  finite  theory.  (A  cruciform 
column  is  equivalent  to  four  flanges  vlth  a  common  hinge  line). 
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Figure  3  shows  similar  results  for  simply  supported  plates  buck¬ 
ling  in  compression.  Square  tubes  were  used  as  test  specimens.  The 
dotted  curve  is  the  incremental-theory  result  computed  by  Handelman  and 
Prager.  Although  it  is  reduced  below  the  elastic  value  at  high  values 
of  stress,  it  is  not  reduced  nearly  enough.  The  finite-theory  result 
falls  below  the  stress-strain  curve  in  this  case  euid  is  again  in  good 
agreement  with  experiment. 

These  ccmqparlsons  tend  to  confirm  the  finite  type  of  law  for  plas¬ 
ticity  and  rule  out  the  incremental  type.  However,  many  plasticity 
experts  reject  the  finite  type  of  law  in  favor  of  the  incremental  type 
and  it  may  now  be  of  interest  to  inquire  into  the  reasons  for  this. 


There  is  first  of  all  a  theoretical  objection  to  finite  laws  based 
on  continuity  considerations  (fig.  U).  The  closed  curve  in  figure  U 
represents  the  original  elastic  limit  of  a  material  being  subjected  to 
compression  emd  shear.  This  means  that  stasrtlng  from  stress  0  no 
plastic  action  occurs  as  long  as  the  combinations  of  shear  and  com¬ 
pression  stays  inside  the  closed  cuirve.  If  now  a  material  is  loaded  in 
pure  compression  up  to  point  A  beyond  the  original  elastic  limit, 
then  the  new  elastic  limit  in  compression  is  A  emd  the  dotted  cui*ve 
indicates  the  supposed  location  of  the  new  elastic  boundary.  If  em 
Increment  of  stress  is  then  added,  bringing  the  stress  to  the  state 
represented  by  point  B  in  the  plastic  region,  the  finite  theory  pre¬ 
dicts  that  an  Increment  of  plastic  shear  strain  dy^  will  be  produced, 
given  approximately  by  the  formula 


where  dr  is  the  increment  of  shear  stress  and  G  is  the  elastic  sheeu* 
modulus.  If  the  increment  of  stress  is  such  as  to  bring  the  material 
to  the  state  represented  by  point  C  in  the  elastic  range,  no  plastic 
action  is  involved  so  the  plastic  shear  strain  is  zero.  One  would 
expect  on  continuity  considerations  that,  as  B  and  C  are  brought 
closer  together  until  they  actually  converge  on  the  boundary,  the  same 
answer  would  be  obtained  whether  we  used  B  or  C.  This,  however, 
does  not  occur  according  to  the  finite  theory  of  plasticity.  Incre¬ 
mental  theories  of  the  types  that  have  been  considered  up  to  the  present 
are  free  from  this  discontinuity  in  behavior  because  according  to  these 
theories  the  plastic  shear  strains  associated  with  AB  and  AC  are 
both  zero. 


As  a  matter  of  fact,  it  is  Just  this  difference  between  the  two 
theories  which  is  responsible  for  the  success  of  the  finite  type  of 
theory  in  accounting  for  the  buckling  of  long  hinged  flanges,  because 
these  are,  for  practical  purposes,  colxmms  ^Ich  fall  by  twisting  and  the 
buckling  stress  of  these  flanges  is  determined  almost  entirely  by  the 
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effective  shear  modulus.  The  plastic  shear  strain  associated  in  the 
finite  theory  with  stress  change  AB,  which  was  considered  theoreti¬ 
cally  objectionable,  gives  us  the  knock-down  factor  needed  to  produce 
agreement  between  theory  and  experiment  for  buckling. 

At  this  point  it  was  decided  to  investigate  the  problem  of  the 
shear  stiffness  of  a  plastically  compressed  material  by  direct  experi¬ 
ment.  A  forged  billet  of  14S-TU  aluminum  alloy  was  machined  into  a 
thin  cylinder  which  was  then  compressed  into  the  plastic  region  8uad 
then  twisted,  with  the  compressive  stress  held  constant.  It  turned  out 
that  the  initial  modulus  at  the  beginning  of  twisting  was  substantially 
the  elastic  modulus  as  predicted  by  the  incremental  theories  6tnd  in 
contradiction  to  the  finite  theories. 

This  posed  quite  a  dilemma.  The  finite  type  of  theory  was  in  good 
agreement  with  plastic-buckling  experiments,  but  disagreed  with  direct 
tests  \inder  controlled  loading  conditions  using  cylinders.  The  incre¬ 
mental  type  of  theory  on  the  other  hand  was  in  good  agreement  with  the 
cylinder  tests  but  in  disagreement  with  buckling  experiments.  There 
appear  to  be  three  possible  explanations  for  this  situation:  (l)  Incre¬ 
mental  type  theories  might  be  right,  but  the  measured  buckling  stresses 
might  be  far  below  the  theoretical  because  of  the  eccentricities  in  the 
specimensj  (2)  the  Handelman  and  Prager  results  for  incremental  theory 
were  obtained  before  Shanley's  concept  of  plastic  buckling  developing 
with  Increasing  load  was  generally  accepted  and  did  not  incorporate 
this  principle,  whereas  the  results  shown  for  finite  theory  do  incorpo¬ 
rate  it;  and  (3)  it  is  possible  that  neither  type  of  plasticity  theory 
is  right  since  each  of  them  appears  to  disagree  with  a  certain  class  of 
experiments.  Let  us  now  examine  each  of  these  possibilities. 

It  can  be  noted  in  figure  2  that  the  buckling  stress  of  long 
hinged  flanges  falls  in  a  relatively  small  scatter  band.  One  would 
expect  that  if  initial ‘eccentricity  were  responsible  for  the  dtscrep-' 
eincy,  the  scatter  band  would  be  wide,  with  at  least  a  few  points  close 
to  the  theoretical  value.  This  makes  it  appear  quite  unlikely  that 
initial  eccentricities  are  responsible  for  the  discrepancy  between  the 
experimental  results  and  the  predictions  of  the  incremental  theory. 
Besides,  there  is  another  explanation  for  the  scatter.  Coupons  cut  from 
different  parts  of  the  cross  section  were  found  to  have  slightly*  dif¬ 
ferent  yield  strengths.  If  the  highest  and  lowest  theoretical  curves 
corresponding  to  these  variations  in  yield  stress  are  plotted  here  on 
the  basis  of  finite  theory,  the  experimental  points  are  all  included 
within  the  two  limiting  curves,  which  are  shown  dashed  in  figure  2. 

The  second  possibility  was  that  the  error  in  the  incremental- 
theory  solution  is  due  to  the  neglect  of  the  Shanley  principle.  But 
incorporation  of  the  Shanley  principle  will  make  no  change  so  far  as 
long-hlnged-flange  results  are  concerned.  The  reason  for  this  is  that 
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the  buckling  stress  of  a  long  hinged  flange  Is  determined  by  shear 
stiffness  alone,  and  according  to  Incremental  theory  the  effective  mod¬ 
ulus  In  shear  of  a  plastically  compressed  material  Is  the  elastic  value 
whether  the  material  continues  to  load  in  compression  or  not.  In  the 
case  of  the  simply  supported  plate  the  incorporation  of  the  Shanley 
principle  in  the  buckling  calculation  does  change  the  results  but  not 
enough  to  bring  the  theory  into  good  agreement  with  the  experiment. 
However,  to  clinch  the  matter,  an  experiment  was  conducted  in  which  a 
simply  supported  plate  was  annealed  in  order  to  give  it  a  very  low  buck¬ 
ling  stress.  The  results  are  shown  in  figure  5  in  which  all  theoretical 
curves  incorporate  the  Shanley  principle.  The  curves  up  near  the  top 
show  the  results  before  annealing.  The  test  data  are  in  excellent 
agreement  with  finite  theory  but  also  in  fair  agreement  with  incre¬ 
mental  theory.  The  bottom  curves  show  the  results  after  annealing.  The 
test  data  are  still  in  fairly  satisfactory  agreement  with  finite  theory, 
but  disagree  seriously  with  the  incremental  theory.  Also,  the  scatter 
is  so  very  small  in  relation  to  the  large  discrepancy  that  it  can  hardly 
be  seriously  supposed  that  the  discrepancy  is  due  to  initial  eccentrici¬ 
ties.  In  further  support  of  this  view,  there  is  the  fact  that  light 
would  not  pass  between  a  straight  edge  and  the  specimen  below  the  buck¬ 
ling  load. 

If  one  concludes  that  laws  of  the  incremental  type  «u:«  wrong  for 
buckling  and  that  finite  laws  are  wrong  for  a  right-angle  loading  path, 
then  one  is  in  the  uncomfortable  position  of  having  no  generally  reli¬ 
able  theory  of  plasticity  on  the  basis  of  which  to  make  structural  anal¬ 
yses,  In  an  effort  to  remedy  this  situation,  a  fresh  start  has  been 
made  and  a  theory  has  been  developed  for  the  stress-strain  relationships 
in  the  plastic  range  on  a  new  basis  (reference  5)«  As  was  mentioned 
before,  finite  and  incremental  theories  in  their  simpler  forms  can  be 
considered  to  be  patteraed  after  the  laws  of  elastic  solids  and  viscous 
liquids,  which  are  treated  as  uniform  homogeneous  mediums.  In  refer¬ 
ence  5  a  metal  was  considered  to  be  an  inhomogeneous  medium,  specifi¬ 
cally  a  polycrystal  ccxaposed  of  a  large  number  of  randomly  oriented 
grains.  Plastic  deformation  was  considered  to  be  due  to  slip  in  those 
grains  for  which  slip  planes  were  most  highly  stressed  in  shear.  The 
resulting  theory  is  a  little  like  the  kinetic  theory  of  gases  in  that 
the  macroscopic  quantities  that  are  of  interest  for  engineering  purposes 
are  obtained  by  summing  the  contributions  of  a  large  nvimber  of  elemen¬ 
tary  physical  processes. 

This  theory  has  been  applied  to  proportional  loading  of  cylinders 
in  biaxial  tension  and  found  to  be  in  excellent  agreement  with  experi¬ 
ment  (reference  6).  All  proposed  plastic  laws  are  reasonably  adequate 
when  applied  to  proportional  loading,  however,  and  our  principal  con¬ 
cern  must  be  how  the  theory  works  out  when  the  stress  ratios  vary 
during  the  loading.  Before  discussing  this  problem,  some  fundamental 
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differences  between  this  and  previous  types  of  plasticity  theory  must 
be  explained. 

Figure  6  shows  a  schematic  representation  of  the  elastic-plastic 
boundary  of  a  materiel,  loaded  in  compression  and  shear  according  to 
conventloned  theory  and  according  to  what  will  be  called  herein  the 
slip  theory.  The  solid  ellipse  Is  In  each  case  the  original  elastic 
' Imlt  of  a  material  without  strain  hardening.  Suppose  the  material  Is 
loeuled  In  compression  to  the  state  A.  As  a  result  of  this  loading, 
conventional  theories  contemplate  a  new  elastic  boundary  something  like 
the  dashed  line  which  separates  the  region  of  loading  from  the  region 
of  unloading.  In  the  ireglon  of  loading  the  plastic  stress-strain  law 
applies,  which  Is  of  either  finite  or  incremental  type.  According  to 
slip  theory,  when  the  material  Is  loaded  to  point  A  there  are  three 
regions:  Regions  (l)  and  (2)  correspond  to  continued  plastic  deform¬ 
ation,  while  region  (3)  Is  the  new  elastic  region  corresponding'  to  the 
unloading  region.  A  small  stress  AB  In  region  (1)  Induces  plastic 
deformation  In  only  those  grains  of  the  material  which  already  experi¬ 
enced  plastic  deformation  while  going  from  0  to  A,  so  that  all  plas¬ 
tically  loaded  grains  continue  loading  or  at  least  do  not  mload.  The 
equations  of  the  slip  theory  Imply  that  a  finite  stress-strain  law 
should  be  valid  In  this  region.  In  region  (2)  some  plastically  defoimed 
grains  are  loading  and  some  are  unloading.  In  this  region  the  slip 
theory  Is  equivalent  to  a  mixture  of  finite  and  Incremental  types  of 
theory,  being  mainly  finite  at  the  upper  end  and  mainly  Incremental  at 
the  lower  end. 

As  a  result  of  this  transition.  In  some  respect  slip  theory  agrees 
with  finite  theory.  In  others  with  Incremental  theory,  while  In  still 
other  respects  It  agrees  with  neither  one.  For  example,  according  to 
finite  theory,  the  plastic  strain  In  the  stress  state  A  Is  Independent 
of  the  loading  path.  According  to  Incremental  theories  the  plastic 
strain  depends  on  the  loading  path,  and  In  particular  a  loading  path 
which  Is  always  on  the  right-hand  side  of  the  compression  axis  would 
lead  to  a  different  strain  from  a  loading  path  which  is  always  on  the 
left.  According  to  the  slip  theory  the  strain  at  point  A  is  Inde¬ 
pendent  of  the  loading  path  provided  the  loEullng  path  stays  within 
region  (3)  but  will  depend  on  loading  path  If  any  part  of  It  traverses 
region  (2).  As  a  second  example,  let  us  consider  the  Initial  shear 
stiffness  at  stress  A  associated  with  stress  Increment  AB.  According 
to  finite  theory  this  Is  below  the  elastic  value  by  an  amount  which 
depends  on  point  A,  but  not  on  the  slope  of  the  line  AB,  provided  AB 
points  upwards  or  corresponds  to  continued  locMllng.  The  Incremental 
theory  gives  the  elastic  modulus  G  under  all  circumstances.  The  slip 
theory  agrees  with  the  finite  theory  In  region  (1),  with  the  Incremental 
theory  In  region  (3);  and  makes  a  transition  between  these  results  In 
region  (2). 
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This  particular  feature  of  the  slip  theory  suggests  a  possible 
explanation  of  the  dllemsa  encountered  earlier  In  this  paper.  According 
to  the  Shanley  principle,  plastic  buchllng  can  be  expected  to  occur  with 
Increasing  load.  One  may  expect  the  corresponding  stress  Increment  to 
lie  In  region  (l),  so  that  the  stiffness  and,  therefore,  the  buckling 
stresses  are  those  predicted  by  the  finite  lavs  of  plasticity.  The 
eiqwrlments  on  plastic  cylinders,  however,  were  carried  out  In  region  (2) 
In  which,  according  to  slip  theory,  the  predictions  of  the  finite  lavs 
are  no  longer  valid. 

This  suggested  explanation  has  been  subjected  to  direct  test  In 
some  further  experiments  performed  at  the  Langley  Laboratory.  A  group 
of  cylinders  vere  compressed  to  point  A  and  then  while  the  compressive 
load  continued  to  Increase  they  were  loaded  In  torsion.  In  this  way 
stress-strain  data  were  taken  with  the  stress  Increment  AB  having  a 
number  of  different  slopes,  some  In  region  (2)  and  some  In  region  (1). 

The  Instrumentation  was  carefully  worked  oyer  to  get  maximum  sensitivity 
In  order  to  get  a  good  measurement  or  the  effective  shear  stiffness  at 
the  very  beginning  of  twisting.  It  was  found  that  the  Initial  shear 
stiffness  was  the  elastic  value  G  for  all  cases  In  region  (2)  and  also 
for  a  case  near  the  border  here  In  region  (1).  In  the  case  of  the  very 
steepest  stress  increment,  a  rather  remarkable  thing  happened;  the  effec¬ 
tive  shear  stiffness  was  the  elastic  value  for  the  first  50  or  100  pounds 
per  square  inch  of  applied  shear  stress  and  then  very  rapidly  drppped  off 
to  the  reduced  value  predicted  by  finite  theory.  This  suggests  that  the 
slip  theory  may  be  qualitatively  coirect  about  the  effective  shear  stiff¬ 
ness  but  that  quantitatively  It  leaves  something  to  be  desired. 

The  slip  theory  has  also  been  tested  by  Investigating  the  validity 
of  Its  fundamental  physical  assumptions  by  observing  a  polished  and 
etched  specimen  of  aluminum  under  a  microscope  during  tensile  testing. 

This  study  conflnned  some  of  the  basic  concepts  of  the  theory  and 
sxiggests  that  others  may  need  refinement  if  accurate  results  are  to  be 
obtained. 

SumBlng  up  what  has  been  said  about  theories,  one  may  say  that 
finite  theories  are  not  satisfactory  when  very  abrupt  changes  occxu*  In 
the  loading  path  such  as  the  right-angle  turns  investigated  In  the 
plastic-cylinder  experiments.  Because  of  this  many  experts  have  tended 
to  reject  finite  theories  In  favor  of  Incremental  theories  except  for 
proportional  loading  where  the  two  types  of  theory  agree.  However,  Incre¬ 
mental  theories  of  the  types  proposed  up  to  this  time  are  unsatisfactory 
for  buckling.  A  new  type  of  theory  based  on  a  physical  model  of  a  plas¬ 
tically  deforming  metal  suggests  that  finite  and  incremental  theories  each 
contain  limited  range  of  application.  There  are  some  prospects  that  the 
new  approach  might  lead  to  the  generally  valid  fonnulation  of  stress- 
strain  lavs  In  the  plastic  range  that  Is  at  present  lacking.  But  this 
would  require  some  refinements  In  the  theory  and  we  cannot  be  sure  that 
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they  can  he  successfully  incorporated.  In  the  meantime  it  appears  that 
there  is  now  theoretical  as  veil  as  experimental  Justification  for  using 
a  finite  law  in  a  wider  variety  of  problems  than  has  been  generally  con¬ 
ceded  up  to  now,  which  is  very  fortunate  since  this  is  the  simplest 
type  of  lav  to  use  in  a  structural  ancdysis. 
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Flguw  1.-  Stress- strain  curve  for  sin5>le  tension, 


ELASTIC  BUCKLING  STRAIN 

Figure  2,-  Buckling  stresses  of  long  hinged  flanges 
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Figure  B***  Buckling  stresfles  of  square  simply  supported  plates. 
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Figure  Losdlxig  and  tinloading  In  finite  and  Incronental  theories 
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Figure  Buckling  of  heat-treated  and  of  annealed  1^  alumln\im 

alloy  plates. 
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Figure  6.<-!  CoiBparlson  of  slip  and  conventional  theories  vlth 
respect  to  elastic  and  plastic  regions. 
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SOME  APPLICATIONS  OF  PLASTICITI  TO  SIIRUCTURAL  ANALYSIS 

By  Elbrldge  Z.  Stowell,  George  E.  Griffith, 
George  J.  Heimerl,  and  John  E.  Duberg 

Langley  Aeronautical  Laboratory 


Three  applications  of  plasticity  theory  are  discussed:  the  first 
concerns  stress  distribution,  the  second  concerns  maximum  strength,  6uid 
the  third  is  concerned  with  the  selection  of  efficient  materials. 

Stress  distribution.-  It  is  well  known  from  both  theory  and  experi¬ 
ment  that  stresses  are  Increased  by  the  presence  of  a  hole  or  notch.  In 
case  of  a  circular  hole  in  a  vide  plate  under  tension,  the  stresses  are 
undisturbed  some  distance  away  from  the  hole  but  increase  rapidly  as 
the  edge  of  the  hole  is  approached,  finally  reaching  a  value  three  times 
as  large  as  the  stress  applied  to  the  sheet  (fig.  1  and  reference  1). 

It  has  been  known  for  a  long  time  that  this  stress  concentration  factor 
of  3  is  reduced  in  the  plastic  range  of  the  material,  but  Just  what 
numerical  value  it  takes  at  any  applied  stress  level  amd  what  the  stress 
distribution  is  have  not  been  known.  Experiments  have  been  conducted 
here  on  24S-T3  aluminum-alloy  sheet  to  help  throw  light  on  these  ques¬ 
tions  (fig.  2  and  reference  2).  As  the  stress  level  increases  so  that 
the  material  at  the  edge  of  the  hole  becomes  plastic,  the  stress  peak 
rounds  off  and  the  stress  concentration  diminishes. 

An  approximate  theoretical  solution  has  been  obtained  by  using 
plasticity  theory  for  the  stress  distribution  in  this  case  by  using  a 
finite  stress-strain  law  (reference  3) .  One  or  two  radii  from  the  edge 
of  the  hole  the  theory  is  not  in  too  good  agreement  with  experiment,  but 
a  long  way  from  the  hole  and  ri^t  at  the  edge,  where  the  stress  is  of 
greatest  interest,  the  agreement  is  very  close.  Figure  3  shows  a  plot 
of  stress  concentration  factor  as  a  function  of  applied  stress  far 
removed  from  the  hole.  The  factor  remains  3  until  the  stress  at  the 
hole  reaches  the  elastic  limit.  As  the  applied  stress  Increases,  the 
stress  concentration  factor  drops  off  toward  unity. 

Fortunately  the  theory  gives  a  very  simple  formula  for  the  stress 
concentration  factor  for  a  circular  hole,  neunely 


In  this  formula  Eg  is  the  secant  modulus  corresponding  to  the  stress 
at  the  edge  of  the  hole  and  Egg  is  the  secemt  modulus  at  the  applied 
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stress.  At  sufficiently  low  stresses,  both  of  these  moduli  are  equal 
to  the  elastic  modulus  so  that  the  vell-Jmovn  value  of  3  is  obtained  In 
the  elastic  range  of  the  material.  By  use  of  this  formula,  stress  con¬ 
centration  factors  for  materials  other  than  24S-T3  can  be  found  If  the 
stress -strain  curve  for  the  material  is  known. 

The  above  formula  for  a  circular  hole  can  be  generalized  in  a 
rather  natural  way  in  order  to  apply  it  to  other  types  of  stress  raiser, 
as  follows: 


where  Kj  is  the  stress  concentration  factor  obtained  in  the  usual  way 

from  the  theory  of  elasticity.  A  test  has  been  made  of  this  formula  by 
experiments  on  large  sheets,  4  feet  by  12  feet,  of  24S-T3  material  with 
semicircular  notches  and  fillets  (fig.  4).  The  plot,  which  is  similar 
to  that  of  figure  3>  is  for  the  sheets  with  the  semicircular  notches  for 
which  the  elastic  stress  concentration  factors  were  about  2,  4,  and  6. 
The  curves  were  computed  from  the  generalized  formula  just  given  and  the 
points  are  experimental.  A  plot  entirely  similar  to  this  plot  was 
obtained  from  the  tests  on  the  sheets  with  fillets. 

Since  a  method  for  finding  the  maximum  stress  at  the  edge  of  a 
hole  is  available,  the  question  may  be  what  happens  when  that  stress 
reaches  the  very  highest  stress  ever  reached  by  the  stress -strain  curve 
for  the  material?  Will  the  specimen  fracture?  Some  preliminary  experi¬ 
ments  made  at  Langley  Laboratory  indicate  that  this  is  so,  at  least  for 
the  light  metal  alloys.  If  further  experiments  confirm  this  result,  a 
method  for  computing  the  fracture  stress  for  perforated  sh^et  for 
aluminum  and  magnesium  alloys,  which  should  be  more  accurate  than 
present  methods  of  calculation,  will  be  available. 

Ifayimuiit  strength.-  The  second  application  of  plasticity  is  to 
problems  of  maximum  strength. 

No  discussion  of  this  subject  would  be  complete  without  mention 
of  the  subject  of  columns.  From  the  point  of  view  of  the  plasticity 
involved,  the  subject  is  much  simpler  than  the  corresponding  problem 
in  the  strength  of  plates  in  that  it  is  a  one -dimensional 

problem,  whereas  plates  are  a  three-,  or  perhaps  more  accurately,  a 
two-and-one-half-dimenslonal  problem.  Although  this  problem  of  column 
strength  has  been  discussed  in  the  literature  for  about  fifty  years, 
it  has  not  been  until  the  last  few  years  that  the  mechanics  of  column 
behavior  €uid  the  plasticity  have  been  correctly  combined  to  show  the 
existence  of  a  maximum  loeid.  Shanley  (reference  4)  was  the  first  to 
point  out  the  exlst«ice  of  a  difference  between  the  buckling  (first 
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bending)  load  and  the  naxlmum  load  in  19^7 and  the  results  of  the 
calculations  (reference  which  show  how  much  the  maximum  load  can 
exceed  the  tangent wnodulus  buckling  load  are  shown  in  figure  The 
factor  by  which  the  tangent Hoodulus  load  must  be  multiplied  to  get  the 
maximum  load  is  plotted  against  the  ratio  of  tangent -modulus  lo6id  to 
the  load  at  the  yield  stress.  The  different  curves  represent  different 
materials  as  defined  by  the  value  of  n  in  the  well-known  Bamberg- 
Osgood  formula.  Ohe  maximum  load  la,  usually,  never  more  than  a  few 
percent  above  the  buckling  load. 

In  the  case  of  a  plate,  the  maximum  load  in  some  cases  may  be  very 
much  more  than  the  buckling  load  because  of  the  additional  support  pro¬ 
vided  by  the  adjacent  structure.  The  simplest  case  of  a  plate  is  that 
of  a  long  hinged  flange  and  this  case  was  studied  first.  Figure  6  . 
shows  the  test  setup;  the  specimen  was  an  extruded  cruciform-section 
column  consisting  of  four  identical  flEmges  mutually  stabilizing  each 
other . 

If  one  analyzes  the  buckling  and  post -buckling  behavior  of  a  plate 
on  the  basis  of  purely  elastic  considerations,  no  maximum  load  Is  ever 
reached.  Figure  7  shows  the  curve  of  load  carried  by  the  flcmge  al  a 
function  of  the  rotation  of  the  midpoint.  Elastic  considerations  lead 
to  the  upper  curve  which  has  no  maximum.  Actually,  the  experimental 
curve  (fig.  7)  shows  a  very  definite  maximum  so  that  a  detailed  consid¬ 
eration  of  plastic  effects  is  required  to  obtain  a  satisfactory  solution 
for  maximum  strength. 

Such  a  detailed  study  has  been  made  (reference  6)  and  the  behavior 
of  the  plate  is  found  to  be  as  indicated  in  figure  ti.  At  successively 
higher  loads,  the  strain  distribution  across  one  flange  will  be  shown 
at  the  left  of  figure  8.  Application  of  the  theory  of  plasticity  will 
then  give  the  corresponding  stress  distribution  shown  at  the  right  of 
figure  8.  Ihe  significant  thing  to  note  is  that,  although  the  strain 
along  the  hinge -line  Increases  indefinitely  as  the  load  Increases,  the 
corresponding  stress  does  not;  it  reaches  a  maximum  and  then  diminishes. 

In  fact,  the  reversal  in  the  hinge -line  stress  occurs  at  the  same  time 
(or  very  nearly  so)  that  the  load  on  the  flange  reaches  a  maximum.  !nie 
hinge -line  stress  falls  off  because  it  is  given  by  the  product  of  the 
hinge-line  strain  and  the  effective  modulus,  and  the  effective  modulus, 
which  depends  on  the  shear  stress  as  well  as  the  conpresslve  stress,  is 
decreasing  more  rapidly  than  the  hinge-line  strain  is  increasing. 

The  strain  emd  stress  distribution  here  correspond  to  a  flange  of 
one  particular  width-thickness  ratio.  From  a  number  of  such  solutions 
for  different  values  of  the  width -thickness  ratio,  a  general  curve  for 
the  maximum  strength  of  a  flange  can  be  drawn  (fig.  9) •  (In  the  ordinate, 
the  critical  stress  has  been  divided  by  the  compressive  yield  stress  in 
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order  to  make  the  curve  nearly  independent  of  material.)  The  points 
represent  experimental  data  obtained  from  a  series  of  24S-T4  cruciform 
section  columns.  An  interesting  feature  of  the  calculations  is  that 
no  matter  what  cruciform  column  in  the  series  is  considered,  the 
combined-stress  intensity  along  the  hinge  line  was  a  constant,  to  about 
one  percent. 

A  similar  analysis  has  been  carried  out  on  the  maximum  strength  of 
a  square  simply  supported  plate.  It  was  assumed  again  that  the  strain 
distribution  was  the  same  as  if  the  material  remained  elastic  and  the 
plastic  stress-strain  relations  were  used  to  compute  maximum  loads. 
Results  for  the  plate  are  compared  with  those  for  the  flange  in  fig¬ 
ure  10.  In  the  lower  region  the  curves  essentially  coincide,  with 
differences  appearing  in  the  upper  region.  Test  points  are  shown  for 
flanges,  Z-,  I-,  and  channel-section  columns.  Since  the  Z-,  I-, 

and  channel-sections  are  con^osed  of  combinations  of  flanges  and  webs, 
one  would  expect  the  test  results  to  be  between  the  curves  for  the 
flange  alone  and  for  the  plate  alone,  and  this  proves  to  be  the  case. 
Again  here,  in  calculating  the  plate  curve,  the  stress  Intensity  along 
the  side  edges  proved  to  be  a  constemt  at  failure. 

'The  applications  of  plasticity  discussed  so  far  have  been  at  room 
temperature.  In  order  to  show  that  the  basic  achievements  at  room 
temperature  also  apply  at  elevated  temperatures,  some  test  on  buckling 
and  maximum  strength  have  been  carried  out.  Figure  11  (taken  from 
reference  7)  shows  test  points  corresponding  to  the  indicated  tempera¬ 
tures  and  superimposed  upon  the  same  theoretical  curves  of  figure  10. 

Selection  of  efficient  structural  material.-  An  interesting  problem 
in  structural  design  is  the  selection  of  the  material  which  will  permit 
the  design  of  the  structure  with  minimum  weight.  The  criterion  to  be 
eD!5>loyed  depends  on  the  design  condition  -  whether  the  structure  is 
being  designed  by  stiffness  considerations,  strength  considerations,  or 
by  considerations  of  buckling  strength  -  and  these  conditions  therefore 
require  separate  treatments. 

Figure  12  shows  stress -strain  curves  for  several  available  struc¬ 
tural  materials,  including  one  magnesium,  two  aluminum,  and  one  titanium 
alloy,  and  stainless  V  steel.  Since  stress  is  load  per  unit  area,  a 
division  of  stress  by  density  will  be  a  measure  of  load  per  unit  weight. 
Figure  13  (taken  from  reference  8)  shows  the  same  stress -strain  curves 
as  in  figure  12  but  with  stress  now  divided  by  density  so  that  the 
ordinate  is  a  measure  of  the  efficiency  with  which  a  load  is  carried. 

In  the  elastic  range,  all  the  materials  lie  nearly  on  the  same  line; 
this  result  is  a  consequence  of  the  familiar  fact  that  all  common  mate¬ 
rials  have  about  the  sane  ratio  of  Young's  modulus  to  density.  Beyond 
the  elastic  limit,  differences  appear. 
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If  yield  strength  is  the  sole  criterion  determining  design,  then 
the  most  efficient  material  is  found  hy  choosing  that  material  with  the 
highest  ratio  of  yield  stress  to  density.  This  quantity  is  indicated 
by  a  tab  on  each  curve  and  it  turns  out  that  titanium  alloy  is  the 
best . 


If  ultimate  strength  is  the  criterion,  ultimate  stress -density 
ratios  may  be  used  in  a  similar  manner. 

If,  however,  stiffness  is  the  sole  criterion  determining  design,  as 
aerodynamic  considerations  sometimes  require,  then  the  most  efficient 
material  is  found  by  choosing  that  material  with  the  highest  ratio  of 
stress  to  density  at  the  permissible  strain.  At  low  values  of  strain 
there  is  little  choice.  At  high  values  of  strain  there  is  a  spread, 
with  titanium  alloy  appearing  to  best  advantage  for  some  values  of 
strain  and  75S-T6  aluminum  alloy  for  other  values . 

If  buckling  is  the  criterion  determining  design,  the  most  efficient 
material  is  found  by  a  slightly  more  complicated  method.  The  ratio  of 
buckling  stress  to  density  is  plotted  against  an  appropriate  structural 
index.  For  column  buckling,  figure  lU  shows  as  abscissa  the  column 
index,  which  is  the  product  of  the  load  P  to  be  carried,  the  end 
fixity  c,  and  the  shape  factor  f,  divided  by  the  square  of  the  dis- 
temce  L  over  which  the  load  myst  be  transmitted.  With  this  index  a 
single  curve  represents  the  buckling  condition  for  any  one  material. 

On  this  basis,  the  choice  again  depends  on  the  value  of  the  index,  with 
titanium  alloy  again  best  at  high  values  of  the  index. 


A  similar  plot  for  plate  buckling  is  shown  in  figure  15,  in  which 
the  index  suitable  for  plates  is  the  product  of  the  load  P  to  be 
carried  and  the  square  root  of  the  buckling  coefficient  k,  divided  by 
the  square  of  the  width  b  of  the  plate.  At  low  values  of  the  index, 
magnesium  alloy  looks  best  and  at  higher  values  T5S-T6  and  then  titanium 
alloy  takes  over.  It  should  perhaps  be  mentioned  that  the  fact  that  a 
curve  for  ZK60-A  is  not  shown  for  structural  index  of  6,  for  example, 
does  not  mean  that  a  plate  could  not  be  designed  for  that  value  of  the 
index.  The  curve  was  discontinued  because  higher  values  of  the  struc¬ 
tural  index  would  correspond  to  large  and  presumably  unacceptable  plastic 
strains,  if  the  plate  were  designed  to  carry  the  load  at  the  buckling 
stress.  However,  any  point  to  the  right  of  a  curve  represents  a  non¬ 
buckling  condition  at  the  corresponding  stress  and  structural  index,  so 
that  for  practical  purposes  each  of  these  curves  may  be  thought  of  as 
terminating  in  a  horizontal  line  which  represents  the  highest  stress 
for  which  one  would  be  willing  to  design  the  plate. 

These  methods  can  be  used  to  select  materials  for  use  at  elevated 
temperatures.  As  mentioned  previously,  plate -strength  research  shows 
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that  methods  used  at  ordinary  tenqperatures  also  apply  at  elevated 
temperatures  provided  the  stress -strain  curve  Is  known  at  the  higher 
tenqperatures .  Comparisons  similar  to  these  can  be  made  at  elevated 
t«qperatures  provided  creep  Is  not  In^rtant  enough  to  change  the  picture 
seriously.  Figure  l6  shows  the  results  of  such  comparison  at  temper¬ 
atures  up  to  800^.  For  tension  or  compression  loading,  tltcmlum  alloy 
^s  best  up  to  about  400°,  with  hl£^-strength  steel  replacing  It 
above  400”.  For  column  buckling,  titanium  and  aluminum  alloys  are  about 
equally  good  up  to  400P,  then  titanium  alloy  takes  over  up  to  800^, 
being  finally  replaced  by  steel  farther  up.  For  plate  buckling  the 
materials.  In  succession,  are  magnesium  alloy,  aluminum  alloy,  titanium 
alloy,  and  steel. 

These  applications  of  plasticity  theory,  which  are  adapted  to  the 
use  of  a  finite  stress -strain  law,  show  that  there  Is  no  need  to  wait 
upon  the  discovery  of  the  most  general  stress -strain  laws  in  order  to 
obtain  Immediately  useful  results  to  some  practical  problems. 
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Figure  Relation  of  oaxtmum  load  for  columns  to  tangent 

modulus  load. 
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Cruciform- section  column  In  testing  machine 
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Figure  T.-  Comparlaon  of  observed  load  on  flange  with  that 
calculated  on  the  basis  of  elastic  theory. 
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Flfore  6.-  Strees  and  strain  distribution  across  a  flange  beyond 

the  buckling  load. 
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Figure  10.-  Relation  of  maxiram  streee  for  both  plate  and  fla^e 
to  critical  streee  cooqpared  with  room  temperature  experiments. 
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Figure  11.-  Relation  of  maximum  stress  for  both  plate  and  flange 
to  critical  stress  compared  with  elevated  temperature  experiments 
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Figure  12.-  Coopressive  stress-strain  curves  for  several  materials 
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Flgiire  13.-  Compressive  stress-density  ratios  plotted  against 
strain  for  the  same  materials  in  figure  12. 


Figure  lU.-  Structural  efficiencies  of  various  materials  for 
column  buckling  at  room  temperature. 
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Figure  15.-  Structijral  efficiencies  of  various  materials  for  plate 

buckling  at  room  temperature. 
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Figure  16.-  Comparison  of  structurally  efficient  materials  at 
elevated  ten^ratures  for  three  types  of  service. 
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FACTORS  GOVERNING  AERODYNAMIC  HEATING 
By  William  J.  O' Sullivan,  Jr. 
Langley  Aeronautical  Laboratory 


^e  first,  step  In  the  solution  of  structural  problems  arising  from 
aerod^amlc  heating  Is  the  detenoinatlon  of  the  aerodynamic  and  thermo- 
dynamic  factors  governing  the  exchange  of  heat  betveen'‘the  structure 
and  the  air.  The  NACA  Is  conducting  research  on  these  factors  at  super¬ 
sonic  Mach  numbers  and  full-scale  Reynolds  numbers  In  free  flight  by 
means  of  a  specially  designed  rocket-propelled  test  vehicle  known  as 
the  lOf-lO.  Already,  experimental  results  In  excellent  agreement  with 
theory  have  been  obtained  for  both  the  steady-state  case  of  thermal 
equilibrium  and  the  transient  case  of  aerodynamic  heating  and  cooling 
within  the  regime  of  immediate  Imjjortance  where  heat  transfer  by  radi¬ 
ation  may  be  neglected.  These  results  are  presented  In  detail  in  refer¬ 
ence  1.  Accordingly,  this  paper  briefly  outlines  the  aerodynamic  and 
thezmodynamlc  factors  involved  and  indicates  the  method  of  applying  the 
results. 

From  the  laws  of  the  conservation  of  energy  and  the  conservation  of 
mass  It  follows  that  when  a  particle  of  air  is  brought  to  rest  without 
energy  escaping  from  it.  Its  temperature  is  increased  because  of  the 
conversion  of  kinetic  energy  into  heat  energy.  The  temperature  attained 
Is  called  the  stagnation  temperature  Tg  and  is  given  by  the  equation 

Tg  =  T  +  T^.2LlJ:jM2  ’  (l) 

where  Tg  and  T  are  absolute  temperatures,  T  Is  the  temperature  of 
the  air  before  being  arrested,  7  Is  the  ratio  of  the  specific  heats 
of  the  air  usually  taken  as  1.4,  and  M  is  the  Mach  number  at  which 
the  air  was  moving  before  being  arrested.  The  second  term  on  the  right- 
hand  side  of  this  equation  Is  the  stagnation-temperature  rise.  In  fig¬ 
ure  1  is  plotted  the  stagnation  tempeirature  In  6igalnst  Mach  number 
for  sea- level  conditions  and  for  the  isothermal  stratosphere  extending 
from  about  35^000  feet  to  105,000  feet  altitude  In  the  NACA  Standard 
Atmosphere  (reference  2)  and  the  NACA  Tentative  Tables  for  the  Proper¬ 
ties  of  the  Upper  Atmosphere  (reference  3}«  To  the  right  of  this  graph 
are  marked  significant  temperatures  for  the  purpose  of  orientation. 

They  are  the  approximate  temperatures  at  which  alumlnxim  and  magnesium 
alloys  lose  90  percent  of  their  room  temperature  strength,  their  melting 
temperatures,  and  the  temperatures  at  which  a  grey  body  exhibits  Incip¬ 
ient,  dull-red,  and  brl^^t-red  glows.  It  Is  evident  that  aerodynamic 
heating  is  structurally  Important  at  Mach  numbers  currently  contem¬ 
plated  for  airplanes  and  already  achieved  by  guided  missiles.  This 
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paper,  however,  is  confined  to  the  lower  supersonic  Mach  numbers,  which 
are  of  Immediate  Importance  where  heat  transfer  by  radiation  may  be 
neglected. 

An  airplane  or  guided  missile  experiences  aerodynamic  heating 
because  the  air  in  contact  with  the  surface  is  brouj^t  to  rest  by  the 
action  of  the  viscosity  of  the  air.  llils  phenomenon  is  illustrated  in 
figure  2.  In  this  figure,  there  is  shown  on  the  left  a  typical  vari¬ 
ation  of  air  velocity  with  distance  from  the  skin  surface  through  the 
entire  thickness  of  the  boundary  layer.  As  the  surface  is  approached 
the  velocity  decreases  from  the  local  velocity  Vj^  Just  outside  the 
boundary  layer  to  zero  at  the  surface.  Accordingly,  the  air  temper¬ 
ature  Increases  as  the  surface  Is  approached,  as  shown  by  the  temper¬ 
ature  profile  on  the  ri^t.  It  la  Interesting  to  note  that  the  most  , 
rapid  decrease  In  velocity  and  Increase  of  temperature  occurs  close  to 
the  surface. 

Consider  first  the  steady-state  case  where  the  aircraft  skin  euid 
the  boundary  layer  have  reached  thermal  equilibrium;  that  is,  no  net 
flow  of  heat  occurs  from  the  boundary  layer  to  the  skin  nor  from  the 
skin  to  the  boundary  layer.  The  skin  and  the  air  in  contact  with  it 
are  not  at  the  stagnation  temperature  but  at  some  lower  teiiq>erature 
called  the  adiabatic  wall  temperature  T^y  because  the  temperature 
gradient  from  the  surface  outward  Into  the  stream  Is  causing  an  outwasrd 
flow  of  heat  within  the  boundary  layer.  The  fraction  of  the  local 
stagnation-temperature  rise  that  actually  occurs  under  this  equilibrium 
condition  is  called  the  boundary- layer  recovery  factor  3  and  is 
defined  by  the  equation 


"^AW  - 
Tg  -  T 


The  numerator  T^^y  -  Tj^  is  the  actual  temperature  rise  across  the 
boundary  layer.  The  denominator  Tg  -  Tj,  Is  the  local  stagnatlon- 
teiiq>erature  rise  or  the  maximum  possible  rise  that  could  occur. 


In  figure  3  la  shown  in  outline  the  aforementioned  RM-10  research 
vehicle.  The  value  of  the  recovery  factor  3  has  been  measured  experi¬ 
mentally  in  free  flight  at  supersonic  speeds  and  full-scale  Reynolds 
nimbers  in  axially  symmetric  flow  on  the  parabolic  body  of  revolution 
of  the  RM-10.  In  the  graph  of  figure  3  are  shown  time  histories  of  the 
Mach  number  M,  stagnation  temperature  Ts.  and  skin  temperature  at  a 

station  near  the  nose  T^  and  near  the  tall  T^  during  a  typical 

flight.  The  model  Is  propelled  from  rest  to  Its  maximum  Mach  number  In 
about  3*3  seconds  by  an  Internal  rocket  motor  and  thereafter  coasts 
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down  throui^  the  speed  range.  At  the  Instant  of  maxlmuai  skin  temper¬ 
ature  at  any  point  on  the  skin  there  Is  no  net  heat  flow  between  that 
point  of  the  skin  and  the  bovindary  layer  as  evidenced  by  zero  value  of 
the  rate  of  change  of  temperature  of  the  skin.  At  this  Instant,  the 
conditions  are  the  same  as  In  the  steady-state  case  and  provide  a  means 
of  deteimlnlng  the  recovery  factor  3.  Since  there  Is  no  net  heat  flow 
at  this  Instant  It  must  be  that  the  air  In  contact  vl1;h  the  skin  Is  at 
the  same  temperatuz^  as  the  skin.  Thus,  the  adiabatic  wall  temperature 

Is  known.  The  stagnation  tes^rature  Is  known  from  the  measured 
flight  conditions,  and  the  temperature  just  outside  the  boundary  layer 
can  be  fo\ind  either  by  calculating  the  flow  about  the  body  by  theory 
alone  or  by  theory  with  the  aid  of  experimental  pressure  distributions. 
This  method  of  determining  the  recovery  factor  assumes  that  the  longi¬ 
tudinal  flow  of  heat  throiigh  the  skin  Is  negligible,  which  point  was 
Investigated  and  found  to  be  the  case  In  these  tests. 

The  experimental  values  of  recovery  factor  so  obtained  are 
presented  in  figure  1^.  The  recovery  factor  3  is  plotted  against 
distemce  from  the  nose  of  the  model.  Comparison  Is  meuie  with  the  theo¬ 
retical  values  (references  4  and  5)  based  on  conditions  Just  outside 
the  boundary  layer.  The  results  are  in  excellent  agreement  with  theory. 

For  structural  purposes,  the  equilibrium  or  steady-state  temper¬ 
ature  ultimately  attained  In  steady  flight  can  readily  be  calculated 
for  airplanes  and  guided  missiles  having  conventional  structures  con¬ 
sisting  of  relatively  thin  reinforced  skins  that  offer  negligible  paths 
for  the  flow  of  heat  through  the  structure.  The  equation  defining  the 
recovery  factor  (equation  (2) )  need  only  be  solved  for  the  adiabatic  wall 
temperature  T^y,  which  is  the  equilibrium  temperature.  The  theoretical 
values  of  the  recovery  factor  may  confidently  be  employed  on  slender 
bodies  like  that  of  the  RM-10,  and  it  is  observed  that,  since  the  lami¬ 
nar  and  turbulent  values  differ  by  less  than  6  percent,  there  is  evi¬ 
dently  little  need  in  a  practical  problem  for  knowing  precisely  the 
point  of  transition  from  laminar  to  turbulent  boundary  layer. 

In  figure  5  is  considered  the  transient  case  where  the  aircraft 
skin  and  the  boundary  layer  are  not  in  thermal  equilibrium  so  that  the 
teiiq)erature  of  the  skin  is  either  increasing  or  decreasing.  Shown  in 
this  figure  is  an  element  of  skin.  This  skin  element  is  assumed  to  be 
sufficiently  distant  from  any  stiffeners  or  bulkheads  in  contact  with 
the  skin  so  that  no  flow  of  heat  takes  place  between  this  skin  element 
and  any  stiffeners  or  bulkheads.  Further,  let  it  be  assumed  that  the 
tenqperature  variation  along  the  skin  is  small  so  that  heat  flow  along 
the  skin  is  negligible.  In  short,  any  heat  flowing  into  or  out  of  this 
element  of  skin  is  assumed  to  come  from  or  return  to  the  boundary  layer 
in  contact  with  it.  Under  these  conditions,  the  rate  of  heat  flow  ^ 

into  or  out  of  this  element  of  skin  is  given  by 
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vhere  Is  the  temperatxire  of  the  skin  and  Is  the  adiabatic 

vail  temper  attire,  so  that  the  quantity  T^y  -  Ty  Is  the  tenqperature 

potential  that  Is  assumed  to  cause  heat  to  flow  between  the  boundaory 
layer  and  the  skin.  The  rate  of  heat  flow  Is  proportional  to  this 
temperature  difference,  to  the  area  A  of  the  element  of  skin  In  con¬ 
tact  with  the  boundary  layer,  and  to  a  coefficient  h^  called  the 
equivalent  heat-transfer  coefficient. 

Another  equation  can  also  be  written  for  the  rate  of  heat  flow 
Into  or  out  of  this  element  of  skin,  namely 


dt 


(4) 


In  this  equation,  A  Is  the  area,  r  Is  the  thickness,  and  d  Is  the 
mass  density  of  the  element,  so  that  Ard  Is  the  mass  of  the  element. 
This  mass  multiplied  by  the  specific  heat  of  the  skin  material  c,  and 

by  the  rate  of  change  of  the  skin  temperatxxre  gives  the  rate  of 

heat  accumulation  or  loss  of  the  skin  element. 

Equating  the  rates  of  heat  flow  given  by  equations  (3)  and  (4) 
and  solving  for  the  rate  of  change  of  skin  temperatures  gives 


Integrating  this  differential  equation  gives  the  equation  for  the  skin 
temperature 


vhere  Tq  Is  the  Initial  temperature  of  the  skin  before  the  start  of 
fllC^t.  The  specific  heat  of  the  skin  material  c  has  not  been  taken 
outside  of  the  Integral  because  its  value  may  change  somewhat  with 
temperature.  From  this  equation  It  is  evident  that  If  the  time  history 
of  the  trajectory  of  the  airplane  or  guided  missile  Is  available,  the 
tine  history  of  the  temperature  of  the  element  of  skin  can  be  readily 
determined  by  performing  a  stei>-by-step  Integration,  provided  that  the 
value  of  the  equivalent  heat-transfer  coefficient  h^  Is  known. 
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Ube  NACA  has  experimentally  determined  the  value  of  the  equivalent 
heat-tranefer  coefficient  h^  In  free  flight  at  full-scale  Reynolds 
numbers  and  supersonic  Mach  numbers  on  the  slender  pointed  body  of 
revolution  of  the  RM-10  research  vehicle  In  axially  symmetric  flow. 

These  values  were  obtained  by  use  of  equation  (^)  vhereln  all  quanti¬ 
ties  except  the  equivalent  heat-transfer  coefficient  vere  known  during 
those  portions  of  the  flight  in  which  the  skin  temperature  was  changing. 
The  adiabatic  wall  temperature  T^^y  was  found  by  use  of  the  experi¬ 
mentally  determined  recovery  factors  already  discussed,  which  were 
assximed  to  apply  In  the  transient  condition. 


The  correlated  heat-transfer  coefficient  results  are  presented  in 
figure  6.  It  may  readily  be  shown  by  dimensional  analysis  or  other 
theoretical  means  that  the  dimensionless  parameters  governing  the  value 
of  the  equivalent  heat-transfer  coefficient  in  aerodynamic  heating  in 
the  case  of  a  flat*  plate  alined  with  the  wind  are  the  Nusselt  number 
Nu,  the  Prandtl  number  Pr,  and  the  Reynolds  number  Re.  For  a  flat 
plate  at  subsonic  speeds  with  turbulent  boundary  layer  it  is  found  from 
the  research  of  Colburn  (references  6  euad  7)  that  these  parameters  are 
related  by  the  equation 


-  0.8 

Nu  Pr  ^  =  0.0296  Re 


(7) 


It  was  conjectured  that  on  a  slender  pointed  body  of  revolution  in 
axially  symmetric  flow  like  that  of  the  RM-10,  and  which  is  represent¬ 
ative  of  many  missiles  and  parts  of  airplanes,  this  flat-plate  formula 
would  closely  apply.  Accordingly,  the  Nusselt,  Prandtl,  and  Reynolds 
nvimbers,  the  formulas  for  which  are  shown  in  figure  6,  were  determined 
by  taking  the  value  of  the  length  I  as  the  axial  distance  from  the 
nose  to  the  element  of  skin  under  consideration  and  basing  the  remaining 
quantities  on  the  conditions  existing  Just  outside  the  boundary  layer 
at  the  element  of  skin  considered.  These  remaining  quantities  are  the 
coefficient  of  thermal  conductivity  of  the  air  k,  the  specific  heat  of 
the  air  at  constant  pressure  Cp,  the  absolute  viscosity  of  the  air  n, 
the  mass  density  of  the  air  p,  and  the  air  velocity  V.  The  results 

1 

-  7 

are  shown  in  figure  6  where  Nu  Pr  is  plotted  against  Re  for  nine 
stations  along  the  body  ranging  from  8. 9  inches  to  123. 5  inches  from 
the  nose,  obtained  from  the  flights  of  two  models.  For  comparison, 
Colburn's  flat  plate  equation  is  also  plotted  on  the  graph  and  is  the 
straight  line  shown. 

The  excellent  correlation  evidenced  by  the  experimental  values 
forming  a  systematic  variation  demonstrates  that  the  Nusselt,  Prandtl, 
and  Reynolds  numbers  are  Indeed  the  factors  governing  the  value  of  the 
equivalent  heat-transfer  coefficient  within  the  range  investigated. 
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The  close  agreement  hetveen  the  Colburn  equation  for  a  flat  plate  with 
turbulent  bovmdary  layer  and  the  experimental  values  for  the  slender 
body  of  the  RM-10,  which  was  1  foot  in  maximum  diameter  and  12.2  feet 
in  length,  demonstrates  that  for  practical  design  purposes  the  equiva¬ 
lent  heat-transfer  coefficient  can  be  computed  from  the  Colburn  equa¬ 
tion  for  slender  bodies  like  the  RM-10  in  axially  symmetric  flow  at 
supersonic  Mach  numbers  and  full-scale  Reynolds  numbers. 

The  fact  that  Colburn's  flat-plate  formula  is  in  such  close  agree¬ 
ment  with  the  experimental  values  obtained  in  three-dimensional,  axi¬ 
ally  symnetric  flow  on  the  slender  RM-10  body  is  strong  evidence  for 
believing  that,  in  two-dimensional  flow  about  thin,  smooth  airfoils  at 
low  lift  coefficients,  particularly  those  with  sharp  leading  edges,  the 
same  correlation  would  probably  hold.  This  supposition  is  yet  to  be 
verified  experimentally.  On  actual  wings  in  three-dimensional  flow  the 
IKJssibillty  of  important  shock  -  botindary- layer  interaction  is  present, 
the  effects  of  which  are  also  yet  to  be  investigated.  Even  on  a  body 
of  revolution  such  as  that  of  the  RM-10,  angle  of  attack  is  known  to 
influence  markedly  the  boundary  layer  and,  presumably,  also  the  aero¬ 
dynamic  heating.  Even  the  internal  structure  conceivably  can  play  a 
part  since  the  bulkheads  or  other  members  that  attach  to  the  skin  may 
act  as  heat  sinks  during  heating  and  as  heat  sources  during  cooling  so 
that  the  bovindary  layer  might  be  altered. 

Despite  these  foreseen  problems  that  are  yet  to  be  investigated, 
the  results  that  have  been  presented  herein  are  directly  applicable  to 
many  practical  cases  such  as  the  more  critically  heated  sections  of  the 
skin  of  ball! Stic -type  trajectory  missiles.  Further,  it  is  observed 
that  in  the  Nusselt,  Prandtl,  and  Reynolds  numbers  the  only  unknown 
quantity  is  the  proper  value  of  the  length  Z  to  employ  in  treating  a 
body  of  arbitrary  shape  at  any  given  value  of  angle  of  attack.  If  this 
theoretical  viewpoint  is  experimentally  proven  correct,  it  may  be  pos¬ 
sible  by  experimentally  determining  the  equivalent-flat-plate  lengths 
corresponding  to  various  points  on  any  given  body  to  be  able  to  calcu¬ 
late  the  aerodynamic  heating  and  cooling  of  the  body  under  transient 
conditions  for  any  given  trajectory. 

If  figure  3  is  again  examined,  it  is  observed  that  when  the  flight 
is  such  that  the  missile  is  propelled  up  to  the  maximum  Mach  number  in 
a  short  time  and  thereafter  coasts,  as  might  be  the  case  with  air-to- 
air  or  ground-to-air  type  guided  missiles,  the  maximum  skin  temper¬ 
atures  actually  attained  under  such  transient  conditions  are  far  less 
than  the  peak  stagnation  temperatvire.  These  lower  temperatures,  to¬ 
gether  with  the  failure  of  some  metals  to  lose  their  structural  strength 
properties  until  they  have  been  at  an  elevated  temperature  for  an  appre¬ 
ciable  time,  give  promise  of  affording  the  structural  designer  oppor¬ 
tunities  for  effecting  significant  weight  reduction  with  accompanying 
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incxvas*  in  perfomaacn  and  tha  poaalbillty  of  enploying  Baterlala  that 
would  be  ruled  out  on  the  basis  of  a  nontranelent  analysis  of  the  aero- 
dynaalc  beating. 
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Figure  1.-  Variation  of  stagnation  temperature  with  Mach  number. 


Figure  2.-  Illuatratlve  velocity  and  temperature  profiles  of 
a  turbulent  boundary  layer. 
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Figure  3«-  Time  histories  of  Mach  number,  stagnation  temperature 
and  akin  temperature  at  two  stations  during  a  flight  of  the 
NACA  RM-10  research  vehicle. 
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Figure  U,-  Experimental  values  of  recovery  factor  on 
HACA  RM-10  research  vehicle. 
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Figure  5.-  Illuatration  of  physical  concepts  involved  in  eqiiatlons 
governing  heat  transfer  in  aerodynamic  heating. 
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THEORETICAL  ARE  EXPSIIM£SITAL  DETERMIMTION  OF  THERMAL  STRESSES 
By  Richard  R.  Heldenfels  eind  William  M.  Roberts 
Langley  Aeronautical  Laboratojry 


DeslgQlng  an  aircraft  structure  \^lch  vlU  withstand  the  tempera¬ 
tures  resulting  from  aerodynamic  heating  presents  many  problems^  one 
of  ^Ich  Is  the  determination  of  the  thermal  stresses  Induced  In  a 
structure  by  a  given  ten^rature  distribution.  The  purpose  of  this 
paper  Is  to  describe  some  of  the  theoretical  and  experimental  thermal- 
stress  problems  stvidled  by  the  Langley  Structures  Research  Divlslcm. 

To  begin  with,  what  are  thermal  stresses?  Thermal  stresses  are 
those  stresses  which  are  Induced  In  a  structure  as  a  result  of  unequal 
thermal  expeinslon  of  Its  parts.  This  unequal  ejqpansion  may  result  from 
materials  with  different  coefficients  of  thermal  expansion  or  from  a 
temperature  distribution  that  Is  not  uniform.  Unequal  expansion  alone, 
however,  is  not  sijfficient  since  a  linear  distribution  across  the 
structvire  may  result  In  thermal  deformation  without  thermal  stresses; 
however,  such  simple  results  are  rare  in  practice.  In  any  case,  the 
thermal  stresses  depend  only  upon  the  temperature  distribution  and  the 
characteristics  of  the  structure  and  are  independent  of  any  applied 
loads. 


Under  what  conditions,  then,  will  aerodynamic  heating  induce  ther¬ 
mal  stresses  In  an  aircraft  struct\ire?  If  the  structure  is  made  of 
dissimilar  materials,  any  temperature  change,  and  thus  any  aerodynamic 
heating,  may  Induce  thermal  stresses.  If  the  structure  is  of  but  one 
material,  then  temperattire  differences  or  nonviniform  distributions  are 
necessary  to  produce  thermal  stresses.  Aerodynamic  heating  produces 
nonuniform  temperature  distributions  as  a  result  of  transient  condi¬ 
tions;  thiis,  thermal- stress  problems  may  eurlse  from  accelerating  and 
decelerating  supersonic  fll^t.  Although  the  temperature  problem  Is  a 
transient  one,  the  stress  distribution  at  any  Instant  of  time  can  be 
treated  as  a  steady-state  problem  because  stress  changes  6u:e  propagated 
much  faster  than  tenqperature  changes. 

Now  that  the  condltlans  vinder  which  thermal  stresses  are  likely  to 
occur  are  known,  the  next  question  is,  what  can  be  done  about  calcu¬ 
lating  them?  The  answer  is  that  for  a  given  temperature  distribution 
thermal  stresses  can,  in  principle,  be  determined  by  application  of  the 
theory  of  elasticity,  provided  the  material  properties  are  known  for 
the  temperattire  conditions  encoxmtered.  However,  the  problem  is  com¬ 
plicated  by  material  properties  which  vary  with  temi)erature  with  the 
result  that  exact  analytical  soluticxis  usually  ceumot  be  obtained  and 
approximate  numerical  solutions  are  very  tedious.  A  nxmiber  of  stress- 
analysis  procedures,  ranging  frcan  the  very  general  to  the  very  simple 
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and  Involving  various  assumptions  and  approximations^  have  been  estab¬ 
lished  and  this  material  Is  available  in  two  recently  p\d)llshed  papers, 
references  1  and  2. 

Even  if  there  are  theoretical  methods  for  determining  thermal 
stresses,  there  Is  a  need  for  experiiflental  methods  to  establish  the 
accuracy  of  approxlinate  theories  and  to  analyze  those  structures  for 
which  theoretical  analyses  are  Impractical.  Thermal  stresses  have  been 
measured  and  correlated  with  theory,  and,  since  this  work  has  not  yet 
been  published.  It  will  be  described  In  detail. 

How  does  one  measvire  thermal  stresses?  To  answer  this  question, 
commercially  available  wire  strain  gages  were  studied  and  their  charac¬ 
teristics  \mder  various  temperature  conditions  between  8o°  F  and  300°  F 
were  determined.  It  was  found  that  the  biggest  problem  was  to  devise 
a  satisfactory  means  of  differentiating  between  strain  readings  due  to 
thermal  expansion  and  those  d\ie  to  stresses.  Two  methods  were  available. 
One  was  to  provide  Individual  dTommy  gages  with  individual  temperature 
control;  however,  this  method  required  an  impractical  amount  of  appa¬ 
ratus  and  did  not  insure  a  high  order  of  accuracy  because  of  differences 
in  individual  gage  installations.  The  other  method  was  to  calibrate 
the  gages  under  all  combinations  of  stress  and  temperature  anticipated. 
The  inqjortant  resxilts  of  this  calibration  project  are  given  in  figure  1 
which  shows  the  effects  of  stress  and  temperatxrre  on  the  indicated 
strain  of  a  typical  wire  strain  gage  momted  on  a  75S-T6  aluminum-alloy 
specimen. 

Two  quantities  are  needed  to  describe  these  calibration  curves; 
the  slope  of  each  line  which  is  actua3J.y  the  modulus  of  elasticity  of 
the  specimen  and  the  horizontal  intercept,  or  zero  shift,  which  is  a 
measure  of  the  thermal  expansion  of  the  specimen,  the  strain  gage,  and 
the  bonding  material.  The  calibration  tests  revealed  that  the  strain 
gages  yielded  accurate  values  for  the  slopes  over  the  temperat^sre  range 
considered  and  that  there  was  negligible  variation  between  Individual 
gage  installations  on  the  seuae  material.  This  was  not  true  with  regard 
to  the  zero  shift  which  showed  considerable  variation  between  indi¬ 
vidual  gage  Installations.  Therefore,  the  only  calibration  required 
for  the  thermal -stress  investigation  was  the  determination  of  the  zero 
shift;  however,  the  zero  shift  had  to  be  determined  for  each  individxaal 
strain-^ge  installation  in  order  to  get  accurate  test  results.  This 
was  done  by  placing  the  entire  test  specimen,  with  gages  installed,  in 
a  furnace,  cycling  the  temperature  to  stabilize  the  gages,  and  then 
measTxring  the  zero  shift  for  the  various  temperatures  anticipated.  This 
zero-shift  calibration  of  each  stialn-gage  installation  is  a  distinct 
disadvantage  in  practical  axiplicatlons,  but  it  was  found  necessary  for 
accurate  re  stilts.  If  ezrors  of  1000  pounds  per  square  inch  are  per¬ 
missible,  however,  an  average  value  of  the  zero  shift  can  be  used. 
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Knowing  how  to  interpret  the  gage  readings  for  a  simple  state  of 
temperature  and  stress,  the  next-  step  was  to  see  how  satisfactory  this 
method  of  calibration  was  for  measuring  thermal  stresses.  Figure  2 
shows  the  specimen  used  for  such  a  thermal  stress  test.  It  was  a 
i  -inch  75S-T6  alluminm-alloy  plate  measuring  2h  by  36  inches  that  was 

eurranged  to  provide  a  simple  two-dimensional  distribution  of  thermal 
stress  under  steady-state  conditions.  It  was  heated  along  the  longi¬ 
tudinal  center  line  by  an  electrical  heating  element  and  cooled  along 
the  longitudinal  edges  by  water  flowing  through  plastic  tubes.  The 
sides  eind  edges  of  the  panel  were  insiilated  to  minimize  heat  loss 
through  these  surfaces  and  to  insure  that  the  strain  gages  were  at  the 
ssune  temperature  as  the  panel.  Strain-gage  rosettes  and  thermocouples 
were  mounted  on  one  quadrant  of  the  panel  and  the  flow  of  cooling  water 
was  controlled  so  that  the  temperature  distribution  was  symmetrical  in 
both  transverse  and  longitudinal  directions.  Since  the  thermal  stresses 
tended  to  buckle  the  center  of  the  plate,  it  was  mounted  on  a  wooden  jig 
to  restrain  this  buckling.  In  addition  the  strain  gages  were  mounted  on 
both  sides  of  the  plate  and  their  readings  averaged  so  that  the  test 
results  would  approximate  a  state  of  plane  stress.  (The  next  paper  by 
Myron  L.  Gossard  and  William  M.  Roberts  will  describe  an  analysis  of 
this  plate  ;dien  buckling  is  allowed, ) 

Figure  3  illustrates  the  temperature  increase  of  the  plate  over 
room  temperature,  that  is,  the  temperature  at  which  the  gages  were 
zeroed.  The  line  shown  is  the  temperature  distribution  assumed  for  the 
theoretical  calciilations  of  the  thermal  stresses  and  corresponds  to 
simple  conduction  between  a  concentrated  heat  source  along  the  longi¬ 
tudinal  center  line  and  a  concentrated  heat  sink  along  each  edge.  The 
temperature  difference  between  the  center  line  and  the  edges  is  150°  F. 
The  actual  distribution,  shown  by  the  test  points,  was  slightly  differ¬ 
ent  because  of  the  finite  area  of  the  source  and  sink  and  because  of 
heat  loss  through  the  sides  and  ends  of  the  plate,  but  it  agrees  within 
a  few  degrees  at  the  measured  points.  The  temperature  was  essentially 
constant  in  the  longitudinal  direction. 

Figure  4  Illustrates  the  longitudinal  direct  stress  induced  in  the 
plate  by  the  ten^rature  distribution.  The  squares  and  circles  repre¬ 
sent  test  data  and  the  lines  represent  theoretical  results.  The  stresses 
8tre  shown  at  the  two  stations  where  the  maximum  measured  valxies  were 
obtained  and  the  cc»5)lete  distribution  is  symmetrical  in  both  the  trans¬ 
verse  and  longitudinal  directions.  The  results  show  that  the  hot  parts 
of  the  plate  are  in  compression  emd  the  cold  parts  are  in  tension  as 
would  be  expected  since  the  cold  parts  restrain  the  greater  expansion 
of  the  hot  parts.  The  stresses  at  the  transverse  center  line  approxi¬ 
mate  those  given  by  an  elementary  analysis  of  a  very  long  plate;  how¬ 
ever,  this  peurticular  specimen  was  so  short  that  there  is  a  very  pro¬ 
nounced  secondary  effect  throu^out  the  res};  of  the  plate  because  the 
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stresses  must  go  to  zero  at  the  unrestrained  ends.  The  theoretical 
results  shovn  were  obtained  from  an  approximate  solution  based  on  the 
principle  of  minimum  complementazy  energy.  Althou^  this  Is  an  approxl 
mate  solution.  It  agrees  well  with  the  test  points. 

Figure  ^  Illustrates  the  Induced  shear  stresses.  The  shear 
stresses  are  associated  with  the  effect  of  the  stress-free  end  and 
depend  upon  the  rate  of  change  of  the  stresses  shown  In  figure  4;  how¬ 
ever,  the  agreement  between  theory  and  experiment  Is  still  good.  These 
stresses  are  Euatlsymmetrlcal  about  the  two  center  lines  and  are  thus 
zero  along  these  center  lines  as  well  as  along  the  unrestrained  edges. 

Figure  6  shows  the  transverse  direct  stress  induced  in  the  plate. 
These  stresses  also  arise  because  of  the  stress-free  end  and  are 
symmetrical  about  the  two  center  lines. 

In  conclusion,  fairly  accurate  measurements  of  thermal  stresses 
have  been  obtained  by  use  of  wire  strain  gages  and  a  special  cali¬ 
bration  technique.  However,  the  calibration  technique  has  the  dls- 
advsLntage  that  It  reqiilres  a  special  test  rxm  to  establish  a  zero  shift 
for  each  gagej  th\is  there  is  a  need  for  a  reliable,  self-ccanpensatlng 
strain  gage  for  the  measurement  of  thermal  stresses.  Additional  tests 
are  planned  to  check  the  accuracy  of  approximate  methods  for  the  stress 
analysis  of  built-up  structures. 
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Figure  3.-  Distribution  of  temperature  increase., 


Figure  U,-  Longitudinal  direct  stresses  Induced  by  the 
temperature  distribution. 
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Figure  5*“  Shear  8ti*e88e8  Induced  by  the  temperature 

dletrlbutlon. 


Figure  6.-  Tranaverse  direct  atreeaea  Induced  by  the 
temperature  dlatrlbutlon. 
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THERMIL  BUCKLING  OP  FIAT  PIATES 
Qy  Myron  L.  Gossard  and  William  M.  Roberts 
Langley  Aeronautical  Laboratory 


The  nonuniform  heating  of  supersonic  aircraft  may  cause  parts  of 
the  structure  to  buckle.  Since  buckling  of  the  aerodynamic  surfaces 
may  have  an  adverse  effect  on  the  performance  of  the  aircraft,  the 
designer  should  know  under  vhat  temperature  distributions  buckling  will 
occur.  In  some  cases  it  may  not  be  feasible  to  design  a  structure 
>^lch  will  not  buckle,  and  the  designer  should  know  the  depth  of  the 
buckles  in  order  to  estimate  their  effect  on  the  performance  of  the 
aircraft. 

Figure  1  shows  two  types  of  thermal  bulging  of  plate  panels  that 
nay  arise  from  nonuniform  temperatures.  The  first  type  of  bulging  is 
caused  by  nonuniformity  of  temperature  through  the  thickness  of  the 
plate.  If  the  temperature  varies  linearly  through  the  thickness  emd 
is  constant  over  the  two  surfaces,  an  unrestrained  plate  assumes  a 
spherical  shape  as  Illustrated  in  figure  1(b)  and  the  deformation  tedces 
place  without  the  appearance  of  thermal  stresses.  This  type  of  bulging 
has  the  same  8tructu3ral  significance  as  the  usxxal  initial  out-of-flatness 
that  occurs  in  sheet-metal  structures.  If,  however,  the  edges  of  such 
a  plate  are  restrained,  stresses  will  be  induced  in  the  plate. 

The  second  type  of  bulging  develops  when  the  teiq>erature  is  non- 
uniform  over  the  middle  sxirface  of  the  plate.  If  this  temperature 
variation  is  nonlinear,  for  example,  as  indicated  in  figure  1(c), 
thermal  stresses  will  be  Induced  in  the  plate  regardless  of  the  manner 
of  edge  restraint,  and  the  plate  will  buckle  somewhat  as  shown  in  fig¬ 
ure  1(d)  in  a  manner  analogous  to  that  which  occurs  under  externally 
applied  forces.  Combinations  of  these  two  types  of  bulging  will  gen¬ 
erally  occur  in  aircraft  structures  subjected  to  heating. 

The  temperature  distributions  and  buckle  patterns  in  various  types 
of  stiffened  panels,  under  thermal  conditions  designed  to  simulate  the 
initial  period  of  fllc^t,  have  been  experimentally  investigated  by  the 
Lockheed  Aircraft  Corporation  and  reported  in  Air  Force  Technical 
Report  ^77^.  RACA  in  this  research  has  been  concerned  thus  far 
with  a  theoretical  as  well  as  an  experimental  study  of  buckling  of  a 
flat  plate  subjected  to  a  simple  type  of  nonlinear  temperature  distri¬ 
bution  over  the  sxurfaces,  with  constant  temperature  through  the 
thickness. 

In  approaching  the  thermal  buckling  problem  theoretically,  the 
established  methods  of  the  theory  of  elastic  stability  can  be  applied. 
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provided  suitable  modifications  are  made  to  account  for  thermal  expan¬ 
sion.  The  thermal  conditions  under  vhich  buckling  could  begin  for  a 
flat  plate  can  be  determined  by  small-deflection  buckling  analysis; 
the  calculation  of  buckle  depth  requires  a  large- deflect ion  analysis, 
however.  The  difference  between  these  two  analyses  is  Illustrated  by 
the  following  two  eqviatlons: 


^  -  -BaV^  + 


£  +  _  2  \ 
°  dx^  dy^  ^  by  hx  by) 


The  large-deflection  equations,  first  derived  by  Von  Karman  except  for 
the  temijerature  term,  are  the  full  set  of  equations  as  given.  The 
small-deflection  equations  are  these  two  equations  if  in  the  first  eqxia- 
tion  the  terms  involving  deflection  w  are  removed.  These  terms  account 
for  the  stretching  of  the  midplane,  which  becomes  appreciable  when 
deflections  are  greater  than  about  one-third  the  plate  thickness. 
Although  the  terms  involving  w  in  the  first  equation  are  not  needed 
in  determining  the  critical  buckling  temperature  and  initial  buckle 
shape,  they  must  be  retained  even  for  small  deflections  if  the  magni¬ 
tude  of  the  buckle  is  to  be  determined. 

Exact  solutions  of  these  equations  are,  in  general,  difficult  to 
obtain,  but  approximate  solutions  can  be  made  by  several  methods. 

In  the  studies  of  thermal  buckling  the  energy  method  has  been 
used  to  analyze  a  simply  supported  flat  plate  of  the  dimensions  shown 
in  figure  2.  The  temperature  was  constant  through  the  plate  thickness 
and  distributed  nonllnearly  over  the  surfaces  in  the  tentlike  fashion 
shown  in  figure  2.  This  is  the  same  panel  as  that  Just  discussed  by 
Richard  R.  Heldenfels  and  William  N.  Roberts.  The  critical  buckling 
temperature  and  buckling  problem  were  detenolned  by  Introducing  the 
calculated  thermal  stresses  and  a  deflection  pattern  with  several 
undetermined  parameters  into  the  energy  expression  of  classical  theory 
Afid  by  minimizing  the  total  energy  with  respect  to  the  parameters. 
Deflections  for  temperatures  above  the  critical  (in  the  large-deflection 
range)  were  determined  by  assuming  that  the  shape  of  the  buckle  pattern 
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r«Miln*d  the  sane  as  at  the  critical  temperature  but  Increased  In  magni¬ 
tude  according  to  a  lav  vhlch>  on  a  weighted  average  over  the  plate, 
satisfies  equilibrium. 

The  results  of  this  study  are  shown  In  figure  3  where  the  maxlmim 
buckle  depth  Is  plotted  against  the  temperature  difference  between  the 
center  line  and  the  edges.  The  two  solid-line  curves  are  theoretical 
results,  with  deflection  starting  at  the  critical  temperatiire  given  by 
small-deflection  theory.  The  upper  solid  curve  would  apply  If  diirlng 
buckling  the  edges  of  the  plate  were  required  to  remain  straight  but 
were  othezvise  free  to  move  In  the  plane  of  the  plate;  whereas  the 
lover  solid  curve  applies  to  a  plate  which  Is  stress  free  at  the  edges. 
Test  data  shown  by  the  circular  symbols  were  obtained  from  a  panel 
mounted  so  as  to  provide  sliiq)le- support  and  stress-free  conditions  at 
the  edges;  therefore  the  lower  solid  curve  should  be  compared  with  the 
test  data.  However,  the  test  data  Indicate  an  initial  out-of-flatness 
of  the  plate  and,  therefore,  the  theoretical  result  should  be  adjusted 
to  account  for  the  effect  of  this  out-of-flatness.  When  this  is  done 
In  a  manner  suggested  by  P.  P.  Bijlaard  for  an  Initial  out-of -flatness 

at  the  plate  center  of  13^  percent  of  the  plate  thickness,  the  dashed- 

r 

line  theoretical  cvirve  Is  obtained  and  the  agreement  of  experiment  and 
theory  Is  seen  to  be  fairly  good. 

In  conclusion,  thermal  buckling  in  supersonic  aircraft  may  be 
analyzed  quite  accurately  If  the  temperature  distribution  Is  known. 

The  shape  and  magnitude  of  the  buckle  for  any  temperature  may  then  be 
determined.  This  determination  required  only  simple  modifications  of 
established  methods  but  the  calculation  procedures  required  are  apt  to 
be  conq>llcated  and  tedious. 
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-  Variation  of  buckle  depth  with  teniperature 
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SGMS  AIRCRAPT  STRUCTURAL  PROBLEMS  RESULTHC  FROM  CREEP 
By  Charle*  Libov«  and  George  J.  Helnerl 
Langley  Aeronautical  Laboratory 


At  teiiqieraturee  to  be  encountered  In  supersonic  flight,  tlsie 
bero—s  an  Influential  factor  In  determining  the  strength,  stresses, 
and  defoxmtlons  of  the  aircraft  structure.  This  occurs  throuifh  the 
■echanlsa  of  creep,  vhlch  Is  -a  phencaasnon  that  Is  highly  accentuated 
by  Increased  temperature.  Creep  Is  manifested  by  a  change  In  strain 
vlth  time  even  when  the  load  and  temperature  are  kept  constant. 

The  creep  characteristics  of  a  metal  are  usually  described  by  a 
set  of  creep  curves  of  the  kind  shoun  In  figure  1.  These  are  obtained 
by  testa  on  tension  specimens  vhlch  show  the  change  of  strain  with  time 
for  various  values  of  constant  tensile  load.  The  particular  curves  In 
figure  1  are  for  75S-T6  allumlnum  alloy  at  600^  F  and  are  baaed  on  an 
empirical  foxmula  suggested  In  reference  1.  The  material  vas  stabilised 
at  600^  F  for  100  hours  prior  to  test.  The  extreme  sensitivity  of  creep 
to  stress  should  be  noted.  Similar  sensitivity  exists  with  regard  to 
teaperaturd. 

One  obvious  effect  of  creep  Is  a  constant  change  In  geoamtry  as 
time  progresses.  In  an  aircraft  this  vould  mean,  for  example,  a  con¬ 
tinual  Increase  In  dihedral  of  vlnga  caused  by  ccmpresslve  creep  In  the 
upper  surface  and  tensile  creep  In  the  lover  surface.  Another  conse¬ 
quence  of  creep  Is  stress  relaxation,  vhlch  Is  manifested  by  a  gradual 
reduction  In  stress  vben  the  strain  remains  constant.  Stress  relax¬ 
ation  could  be  beneficial  In  the  aircraft  structure  by  relieving  high 
stress  peaks  resulting  from  thermal  eipanslon.  It  Is  harmful,  hovever. 
In  that  It  can  cause  loosening  of  rivets  and  other  connections  that  are 
sipposed  to  resMln  tl^t.  Creep  further  coipllcates  the  aircraft  design 
picture  by  virtue  of  the  fact  that  Its  strains  are  cumulative  and 
generally  result  In  fracture  In  some  finite  time. 

It  Is  clear,  then,  that  In  designing  the  aircraft  structure  for 
high  tesperatures,  consideration  vlll  have  to  be  given  to  the  factors 
of  pexmdsslble  deformation  and  desired  service  lifetime.  In  these 
respects  the  design  of  aircraft  to  prevent  creep  Involves  the  same 
problems  as  the  design  of  any  other  high- temperature  equipment.  For 
aircraft,  hovever,  one  additional  factor  vlll  have  to  be  considered, 
namely,  the  effect  of  creep  on  the  maximum  load-carrying  capacity  of 
msmibers  subject  to  Instability,  such  m  columns  and  flat  and  curved 
plates.  At  room  temperatures,  many  tests  and  theory  have  Indicated  a 
correlation  betveen  buckling  strength  and  the  stress- strain  curve. 

(For  exmmple,  see  references  2  and  3.)  The  same  sort  of  correlation 
has  been  shown  to  exist  at  high  toperatures  provided  there  Is  no 
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appreciable  creep,  or,  In  other  words,  provided  the  load  is  brought  up 
to  its  buchllng  value  in  a  fairly  short  tine  (reference  Ift 
ever,  the  load  Is  brought  up  to  a  value  below  the  buckling  strength  of 
the  structural  eleaent  and  kept  there  for  sone  time,  as  It  always  would 
be,  and  If  the  stress  Is  high  enough,  creep  will  occur,  and  there  arises 
the  possibility  that  the  eleaent  will  eventually  fall  even  though  the 
load  Is  below  the  short-tlae  buckling  strength. 

Initial  eccentricities  will  play  an  laportant  part  In  problesu  of 
this  kind.  Consider,  for  eaaaple,  a  sll^tly  crooked  column  carrying 
a  load  below  Its  short>tlne  buckling  load.  Because  of  the  crookedness, 
one  side  of  the  coluan  will  be  stressed  more  highly  than  the  other  side. 
Because  of  the  higher  stress,  creep  will  occur  aore  rapidly  on  one  side 
than  on  the  other  side.  This  will  Increase  the  crookedness,  which  In 
turn  will  Increase  the  stress  difference,  and  so  forth.  Instead  of 
resulting  In  a  final  equlllbrlua  position,  this  process  may  diverge  and 
cause  collapse  of  the  colinnn. 


In  problems  of  this  kind,  the  Internal  stress  In  the  structural 
eleaent  will  change  continually  with  tiae.  It  Is  clear  that  for  the 
analysis  of  such  problems  the  conventional  creep  curves  which  are 
obtained  at  constant  stress  are  Instifflclent.  Needed  are  tests  or 
theory  to  describe  what  happens  when  the  stress  changes  with  time. 
Shanley  has  proposed  such  a  theory  which  Is  adalttedlv  approximate  but 
offers  hope  for  engineering  applications  (reference 


Shanley 's  theory  for  uniaxial  tension  or  coopresslon  at  constant 
temperature  Is  Illustrated  In  figure  2.  If  the  material  Is  at  some 
condition  of  stress  and  strain  represented  by  the  point  labeled  (a,c) 
and  the  stress  Is  Instantaneously  changed  by  an  amount  do,  then  the 
strain  Is  assumed  to  change  elastically  by  an  amount  dc  according  to 
the  elastic  modulus  for  the  particular  temperature  T  of  the  material; 


the  strain  change  can  therefore  be  written  de 


and  It  Is 


represented  by  the  arrow  that  Is  parallel  to  the  elastic  line.  On  the 
other  band.  If  the  stress  is  kept  constant  for  a  short  Interval  of 
time  dt,  a  small  strain  change  dc  occurs  due  to  creep  and  Is  repre¬ 
sented  by  the  horizontal  arrow.  The  rate  at  which  this  strain  occxirs 
Is  assumed  to  be  a  function  only  of  the  Instcmtaneous  values  of  the 
stress,  the  strain,  and  the  temperature  and  not  of  how  these  conditions 
were  attained.  For  constant  stress,  therefore,  one  can  write 
dc  ■  f(o,c,T)dt,  where  dt  Is  a  small  time  Interval  and  f(o,c,T)  Is 
a  known  fuctlon. 


Shanley*  s  formulation  of  a  creep  theory  would  be  very  desirable 
from  an  engineering  viewpoint  If  It  proves  to  be  sufficiently  accurate 
because  any  loading  hletory  can  be  approximated  by  a  series  of  steps,  as 
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•hovn  In  figure  3*  Each  step  consists  of  an  instantaneous  stress 
change  by  an  aaount  da  followed  by  a  constant- stress  period  of 
length  dt.  The  strain  changes  occurring  during  each  of  the  two  parts 
of  the  step  nay  be  added  together  to  obtain  the  total  strain  change 

d«  «  do  +  f(a,e,T)dt 

Dividing  throu^  by  dt  gives  a  differential  equation  for  the  strain 
rate,  which  applies  in  the  limit  as  the  steps  become  infinitesimally 
small.  If  the  tesqperature  as  well  as  the  stress  were  changing,  another 
term  could  be  added  to  the  equation  to  express  the  strain  change  due  to 
thezmal  esqianslon. 

Application  of  the  theory  Jtist  described  to  a  given  structure  of  a 
given  material  at  a  given  tesperature  requires  a  knowledge  of  the  values 
of  the  creep  function  f(a,c,T)  for  many  coBd)lnatlona  of  a  and  c; 
these  values  may  be  obtained  from  conventional  creep  curves.  For 
exaaple,  the  equation  for  the  creep  curves  of  figure  1  is  given  in 
referezMe  1  as 


«-£  +  Ae®®t>^  (1) 

E 

where 

A  -  2.64  X  10"'^ 

B  -  1.92  X  10“3 
E  -  5.2  X  10^ 

K  -  0.66 

By  differentiating  c  with  respect  to  t,  regarding  a  as  constant, 
and  then  eliminating  t  in  terms  of  a  and  c,  the  following 
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•xjpreMlon  for  f(a,c,T)  for  T^~>T6  alumlnun  alloy  at  600°  F  is 
obtained: 


(2) 


In  order  to  Illustrate-  the  type  of  behavior  that  mi^t  be  expected 
of  structural  elements  subject  to  Instability  vhen  creep  enters  the 
plctxire,  an  analysis  has  been  made  of  a  simple  pin-ended  column  of 
7^-T6  aluminum  alloy  at  600°  F  by  \Mlng  Shanley's  creep  theory  In  con- 
Jtmctlon  vlth  equation  (2) .  The  probl«B  Is  that  of  an  H-sectlon  column 
vlth  a  veb  of  negligible  area  and  with  an  initial  crookedness  In  the 
shape  of  a  half  sine  vave.  (See  fig.  4.)  A  load  belov  the  buckling 
load  Is  Instantaneously  applied  and  produces  an  average  stress  a.  The 
load  Is  kept  constant  and  because  of  creep  the  column  continues  to 
deflect  laterally.  The  shape  Is  assumed  to  remain  a  half  sine  vave  and 
plane  sections  are  assumed  to  remain  plane.  The  problem  was  solved  by 
determining  the  time  history  of  the  flaisge  stresses  and  strains  at  the 
center  of  the  column  such  that  the  following  conditions  are  satisfied 
at  all  times:  stresses  In  equilibrium  with  external  loading}  strains 
ccmqmtible  vlth  the  assunqptlon  that  plane  sections  remain  plane  and 
that  the  shape  remains  a  half  sine  wave;  stresses  emd  strains  must 
satisfy  the  strain- rate  equation  of  figure  3*  (in  the  analysis  It  was 
necessary  to  substitute  a  numerical  value  only  for  the  constant  K. 

The  results  of  the  analysis  therefore  apply  to  other  materials  whose 
creep  curves  have  the  foxm  li^lled  by  equation  (l)  with  K  >  0.66.) 

The  growth  of  lateral  deflections  is  Illustrated  In  figure  where 
a  lateral-deflection  parameter  A  Is  plotted  against  a  time  parameter  r 
for  several  values  of  a  strale^tness  parameter  8.  The  main  thing  to 
notice  about  this  figure  Is  that  each  curve  has  a  vertical  asyqptote 
Indicating  that  the  lateral  deflections  approach  Infinity  (resulting  In 
collapse  of  the  column)  In  some  finite  time.  The  asynqptotlc  value  of  t 
will  be  referred  to  as  xcr* 
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The  straightness  parameter  S  is  defined  In  terms  of  the  flange 
stresses  at  the  middle  of  the  column  Immediately  after  the  load  Is 
Instantaneously  applied  and  before  any  creep  has  occtirred;  that  Is, 

- 

S  «  e 


where  Ci  Is  the -stress  on  the  Inside  or  concave  side  of  the  column 
and  Oq  Is  the  stress  on  the  outside  or  convex  side  of  the  column. 
These  stresses  are  obtainable  by  means  of  an  elastic  calculation.  If 
the  column  Is  perfectly  straight,  then  ffi  and  Oq  are  eqxial.  The 
stralfpitness  parameter  S  then  reduces  to  unity  and,  as  the  graph 
shows,  there  are  no  lateral  deflections  with  time  for  S  =  1.  All 
degrees  of  Initial  crookedness  or  straightness  are  included  within  the 
range  S  »  0  to  S  =  1. 

The  time  parameter  t  Is  proportional  to  the  time  in  hours  but 
also  Involves  other  quantities: 


The  parameter  Is  extremely  sensitive  to  a^,  which  appears  in  the 
e3q)onent  of  e,  and  therefore  to  the  Initial  eccentricity  and  the 
closeness  of  the  applied  stress  to  the  Euler  stress.  The  difference 
between  the  applied  stress  a  and  the  Euler  stress  Oq  also  appears 
explicitly  in  the  denominator. 

The  ratio  of  the  lateral  deflection  d  to  the  width  b  between 

flange  center  lines  is  expressed  as  the  sum  of  two  terms:  q* 

the  valiie  existing  immediately  after  the  load  Is  applied  and  obtainable 

by  an  elastic  calculation,  and  (  ,  which  Is  the  time- dependent 

\D/CREEP 

part.  It  Is  this  latter  term  that  is  given  by  the  parameter  A  through 

the  equation  j  = 

\b/c3fflEP  2Ba 

If  the  lateral  deflections  themselves  are  not  of  interest,  but  only 
the  service  lifetime  of  the  column  Is  important,  this  Infozmatlon  can  be 
given  In  another  plot  that  gives  Tq^.,  the  lifetime  parameter.  In  terms 
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of  S,  the  straightness  parameter.  This  curve  Is  shovn  In  figure  6. 

The  curve  goes  up  to  Infinity  asynqptotlcally  as  the  column  approaches 
straight. 

The  H-sectlon  column  In  Itself  Is  not  of  too  great  Interest. 
Nevertheless,  analyses  of  the  type  described  may  be  qualitatively 
Important  In  Indicating  the  significant  parameters  that  enter  Into  a 
problem.  With  the  significant  parameters  knovn,  the  number  of  tests 
required  to  obtain  design  Information  for  structures  not  easily  analyzed 
can  be  greatly  reduced,  ^e  results  of  a  limited  number  of  tests  on 
solid  rectangular-section  columns  have  been  published  In  reference  1. 

If,  for  the  solid  rectangular- section  column,  the  flange  stresses  a± 
and  oo  reinterpreted  to  be  extreme  fiber  stresses,  then  the 
results  of  reference  1  can  be  plotted  In  terms  of  the  parameters 
and  S,  which  were  found  to  be  significant  for  the  H-sectlon  column. 

The  results  are  shown  In  figure  7*  Some  correlation  Is  Indicated,  but 
the  number  of  tests  Is  too  few  to  conclude  that  these  parameters  aze 
Indeed  significant  for  the  solid  column.  The  dashed  curve  Is  taken 
from  the  previous  figure  and  replotted  to  the  present  scale.  The  trend 
of  four  of  the  tests  agrees  with  that  of  the  dashed  curve.  (The  solid 
bar  throTigh  one  of  the  test  points  Indicates  permissible  leeway  In 
plotting  this  point  because  of  an  Tincertalnty  In  the  Initial  crookedness 
of  the  colzmm.)  The  test  points  for  the  solid  column  fall  above  the 
theoretical  curve  for  the  H-sectlon  column,  as  Is  to  be  expected.  The 
two  test  points  on  the  left-hand  side  of  the  graph  which  definitely  do 
not  follow  the  trend  of  the  other  test  points  represent  tests  in  which 
the  applied  load  was  very  close  to  the  Euler  buckling  load.  The 
applied  stress  In  these  two  cases  was  the  same  and  was  9^  percent  of  the 
Euler  buckling  stress  as  contrasted  with  a  maxlxnum  of  84  percent  for  the 
other  tests.  Since  the  behavior  of  the  column  is  very  sensitive  to  the 
closeness  of  the  applied  load  to  the  Euler  buckling  load,  a  small  error 
In  the  calculated  Euler  buckling  load  could  cause  considerable  scatter. 
Furthermore,  the  extreme- fiber  stress  on  the  concave  side  of  the 
column,  as  computed  by  an  elastic  analysis,  turned  out  to  be  above  the 
prox>ortlonal  limit  for  these  two  cases.  The  finite  length  of  time  It 
takes  to  load  the  column  would  then  be  significant,  and  the  present 
analysis  which  assumes  elastic  behavior  during  Instantaztfous  loading’ 
would  not  apply.  Some  Indication  of  the  large  experimental  scatter  to 
be  expected  In  creep-buckling  tests  Is  afforded  by  the  fact  that  these 
two  test  colimms  were  Identical  except  that  one  was  more  crooked  than 
the  other,  yet  the  more  crooked  column  lasted  longer  than  the  stralghter 
one. 


The  column  problem  Just  described  Is  one  of  the  sliiq)lest  problems 
one  can  formulate  Involving  buckling  and  creep.  Further  experimental 
and  theoretical  work  lis  needed}  tests  of  plates  as  well  as  additional 
tests  of  coltmns  conteaqplated. 
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75S-T6  aluminum  alloy  at  600°  y 


Figure  2,-  Sketch  illustrating  basic  nmumi.t  t  ^  « 
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Figure  5.-  Growth  of  lateral  deflection  of  Initially  crooked 
column  under  constant  load. 


S  (STRAIGHT 

COLUMN) 

Figure  6.-  Chart  showing  yariatlon  of  lifetime  parameter  Tqj.  v 
straightness  parameter  S  for  initially  croolced  column  under 
c(»istaiit  load. 


COilFIUlTIAL 


o 


conmzRZAL 


U5 


20H 


o  EXPERIMENT  (RECTANGULAR  SECTION) 
- THEORY  (H-SECTION) 


'"Ter  I 


lOH 


o 


0 

-2ilO^  -10^ 


(|lns)o(«(oi,-oi) 


'  '  ' 

“I03 


Figure  7»-  CoogiariBon  of  H-sectlcm  theory  with  solid- rectangular 

section  tests. 
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SIZE  EFFECT  IN  NOTCH  FATIGUE 
By  Paul  Kuhn  and  Herbert  F.  Hardrath 
Langley  Aeronautical  Laboratory 


The  present  state  of  knowledge  on  fatigue  is  characterized  by  a 
practically  complete  lack  of  quantitative  lavs  or  rules;  no  one  quan¬ 
tity  can  be  calcvilated  from  any  other  quantity  with  reasonable  accuracy. 
Reliable  allowable  stresses  can  therefore  be  established  only  by  means 
of  rather  specific  tests  in  which  all  the  pertinent  variables  are  varied 
over  the  desired  range.  Together  with  the  scatter  found  in  fatigue 
tests,  this  consideration  requires  such  a  vast  amount  of  testing  that  it 
is  is^asible  to  fulfill  all  the  demands  for  information.  In  this 
situation,  any  reasonably  reliable  rule  that  permits  the  elimination  of 
one  parameter  from  the  test  schedules  will  effect  a  large  improvement. 

One  of  the  main  parameters  involved  is  the  factor  of  stress  con¬ 
centration.  For  the  sake  of  brevity,  the  word  "factors”  is  used 
throughout  rather  than  the  words  "stress  concentration. "  For  a  given 
stress  raiser,  such  as  a  notch,  a  theoretical  factor  can  be  obtained 
by  the  theory  of  elasticity  or  by  photoelastic  methods.  Comparative 
fatigue  tests  on  smooth  and  notched  specimens  give  the  actual  or  fatigue 
factor,  which  depends  on  the  stress  level  or  on  the  number  of  cycles  to 
failure;  in  the  special  case  when  the  specimen  fails  in  one- quarter  of 
a  cycle  -  the  ordinary  tensile  test,  it  is  called  the  static  factor. 

This  paper  is  concerned  with  only  the  opposite  limiting  case,  where 
many  millions  of  cycles  are  required  to  produce  failure,  which  implies 
that  no  plasticity  effects  are  involved. 

The  analysis  of  many  fatigue  tests  has  shown  that  the  fatigue  factor 
lies  somewhere  between  the  static  value  and  the  theoretical  one.  Tests 
on  geometrically  similar  specimens  of  differing  sizes  have  shown  that 
the  fatigue  factor  Increases  with  size,  and  it  is  more  or  less  generally 
believed  that  it  approaches  the  theoretical  value  as  sui  upper  limit. 

This  effect  of  size  might  be  called  the  geometric  size  effect,  because 
it  depends  more  on  geometry  than  on  the  material.  It  should  not  be 
confused  with  what  might  be  called  the  material  size  effect,  the  dif¬ 
ference  in  properties,  say,  between  a  sheet  l/l6  of  an  inch  thick  and 
a  plate  U  inches  thick,  eurlsing  from  the  fact  that  the  hot  and  cold 

work  necessary  to  convert  a  4- inch  plate  into  a  ^-inch  sheet  Improves 

16 

the  properties  of  the  material  (at  least  the  static  properties).  This 
material-size  effect  can  be  eliminated  from  the  fatigue  tests. 

The  existence  of  the  geometric  size  effect  has  been  definitely 
recognized  for  about  20  years.  The  most  plausible  explanation  afforded 
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l8  that  the  theory  of  elasticity  breaks  down  when  applied  to  small 
regions  of  certain  materials  such  as  structviral  metals.  The  theory  is 
based  on  tvo  a88uiq>tlons:  (a)  the  material  is  elastic  cmd  (b)  it  is 
homogeneous.  Hov,  metals  have  a  granular  structure.  Different  grains 
may  have  different  properties,  and  in  an  individual  grain,  the  properties 
may  be  directional.  The  properties  of  large  pieces  of  the  material  are 
determined  by  the  statistical  average  of  the  individual  properties. 
Hovever,  when  the  piece  under  consideration  includes  only  a  few  grains, 
it  is  necessary  to  consider  the  individual  properties.  At  this  level, 
the  assumption  that  the  material  is  homogeneous  is  no  longer  tenable, 
and  the  usual  theory  of  elasticity  is  no  longer  applicable. 

A  theoretical  solution  which  takes  into  account  simultaneously 
over-all  geometry  as  well  as  the  individual  properties  of  the  grains 
seems  rather  hopeless,  and  if  it  existed,  it  would  undoubtedly  be  too 
complex  for  practical  use.  An  interesting  suggestion  for  a  practical 
approximation,  however,  has  been  made  by  Neuber,  who  has  done  a  leurge 
amount  of  work  on  the  mathematical  theory  of  notches.  Neuber  assumed 
(reference  1)  that  the  specimen  acts  as  though  the  material  were  made 
up  of  blocks,  or,  say,  building  bricks,  which  have  the  property  that 
no  stress  gradient  can  develop  across  a  face  of  the  block.  By  applying 
this  concept  to  the  critical  region  at  the  bottom  of  a  notch,  Neuber 
arrived  at  the  general  formula  for  correcting  the  stress  concentration 
given  by  the  classical  theory  for  size  effect 


Kn  -  1 

%  -  1  +  - - - 

«  |r 

In  this  formula,  %  (with  the  subscript  N  for  Neuber)  denotes 
the  corrected  factor;  Kji,  the  theoretical  factor;  cd,  the  flank  angle 

of  the  notch;  R,  the  radius  of  the  notch;  and  R' ,  the  dimension  of  the 
Neuber  "building  block,”  which  is  dependent  to  some  extent  on  the 
material.  If  R'  approaches  zero,  the  Neuber  factor  becomes  equal  to 
the  theoretical  factor.  If  R'  becomes  very  large,  the  Neuber  factor 
becomes  equal  to  unity;  that  is,  the  notch  has  no  stress-concentration 
effect.  In  terms  of  the  widely  used  concept  of  "notch  sensitivity," 
a  small  valvie  of  R'  denotes  large  notch  sensitivity  emd  a  large  value 
of  R'  denotes  small  notch  sensitivity.  However,  the  formula  brings 
out  the  fact  that  notch  sensitivity  is  not  a  p\ire  material  property, 
but  depends  on  size. 

In  order  to  assess  the  applicability  of  the  Neuber  formula  to  the 
fatigue  problem,  a  collection  was  started  of  all  notch- fatigue  tests 
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that  could  be  fovind  In  the  literature;  particular  attention  was  paid 
to  tests  In  which  some  size  factor  was  varied  systematically.  The 
collection  to  date  Inclvides  references  2  to  1^.  As  a  first  approxl- 
matlon,  R'  was  assiuned  to  have  the  constant  value  0.02  Inch,  which 
was  calculated  by  Neuber  from  two  sets  of  static  strain  measurements 
on  steel  specimens.  The  use  of  a  constant  value  of  R*  Is  obviously 
only  a  crude  first  approximation;  nevertheless  the  results  were 
encouraging.  For  the  majority  of  the  tests,  the  use  of  the  Neuber 
correction  improved  the  prediction  of  the  fatigue  factor. 

A  stiidy  of  the  data  obtained  with  the  first  approximation  suggested 
that,  for  steel  specimens,  the  value  of  R'  might  be  taken  to  be  a 
function  of  the  tensile  strength  of  the  material.  The  tentative  cvirve 
is  shown  in  figure  1.  As  the  strength  of  the  steel  Increases,  R* 
decreases  -  that  is,  the  notch  sensitivity  Increases.  Judging  by  the 
data  available  now,  this  tentative  ctirve  appears  to  be  reasonably 
reliable  for  \iltimate  strengths  above  100  ksi;  for  lower  strengths,  it 
appears  to  be  less  reliable. 

The  next  two  figures  show  some  sample  comparisons,  chosen  at 
random,  between  calculated  Neuber  factors  and  experimental  fatigue  factors 
for  sets  of  tests  where  some  systematic  variation  of  size  was  made.  Fig¬ 
ure  2  shows  tests  on  arotatlng  beams  with  semicircular  grooves.  The 
materials  were  S.A.E.  23^5  and  S.A.E.  1020  steel  as  noted.  In  the  two  sets 
of  teats  shown  at  the  top,  the  ratio  of  notch  radius  to  specimen  diameter 
was  constant;  the  theoretical  factor,  shown  by  the  solid  line,  is  there¬ 
fore  constant.  The  Neuber  factor  is  shown  as  a  dashed  line.  The  dif¬ 
ference  between  the  solid-line  curve  and  the  dashed  curve  constitutes 
the  predicted  size  effect.  The  circles  denote  experimental  fatigue 
factors,  obtained  in  the  usvial  manner  from  the  S-N  curv«s  for  smooth 
and  notched  specimens.  In  the  two  te^ts  shown  at  the  bottom,  the  notch 
radius  was  constant;  the  theoretical  factor  therefore  varied.  The 
agreement  between  predicted  and  observed  size  effect  is  very  good  for 
all  cases. 

Figure  3  shows  tests  on  filleted  shafts.  The  agreement  is  very 
good  for  one  set  of  tests  on  mild  carbon  steel  (lower  right)  and  for  the 
nickel-molybdenum  steel.  For  the  other  tests  on  mild  carbon  steel  (upper 
left),  the  agreement  is  not  so  good.  However,  the  largest  discrepancy  is 

about  12^  percent,  which  is  of  the  same  order  as  the  uncertainty  of  the 

experimental  fatigue  factor  caused  by  uncertainty  of  the  S-N  curves  on 
which  the  experimental  factor  is  based. 

From  calculations  and  test  data  such  as  shown  in  figures  2  and  3> 
and  from  other  notch-fatigue  tests  where  no  size  variations  were  made, 
three  sunnary  plots  have  been  prepeured.  Figure  4  summarizes  all  tests 
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available  on  rotating  beams  with  grooves  and  V- notches.  The  ratio  of 
Neuber  factor  to  fatigue  factor  is  plotted  against  notch  radius. 

Because  the  number  of  points  is  rather  large,  some  groups  have  been 
averae^.  The  average  is  shown  by  a  circle,  the  number  above  the  circle 
is  the  nuatber  of  points  averaged,  and  tick  marks  indicate  the  highest 
and  the  lowest  points  in  the  group. 

The  plot  shows  that  the  Neuber  factor  agrees  with  the  fatigue 
factor  to  an  accuracy  of  ±10  percent  for  the  large  majority  of  the 
tests.  Outside  of  one  wild  point,  serious  discrepancies  appear  only 
for  radii  less  than  0.01  inch.  All  the  low  points  at  this  radius  are 
tests  run  by  one  investigator  in  one  investigation,  which  leaves  room 
for  the  suspicion  that  there  was  some  systematic  flaw  in  the  test 
technique. 

Figure  5  is  the  same  type  of  summary  plot  for  filleted  shafts. 

Again,  most  of  the  points  fall  within  the  ±10-percent  band.  There  are 
no  low,  or  vinconservative,  points  on  this  plot. 

Figure  6  is  the  svimmary  of  the  axial-load  tests.  The  three  types 
of  specimens  shown  at  the  top  are  the  ones  tested  at  the  Battelle 
Memorial  Institute  for  the  NACAj  their  material  is  S.A.E.  4l30  steel. 

The  radii  on  these  specimens  are  varied  to  obtain  different  theoretical 
factors.  The  points  fall  within  a  ±5-percent  band.  The  lower  specimen 
is  cylindrical  with  a  V-notch.  The  tests  of  these  cylindrical  specimens 
(reference  5)  were  run  on  a  number  of  alloy  steels  with  various  heat 
treatments.  Half  of  the  tests  were  run  at  room  temperature  and  half, 
at  -78°  Fj  but  there  was  no  significant  difference  due  to  temperature. 

The  width  of  the  scatter  band  for  these  tests  is  cons‘iderably  larger 
than  for  the  other  tests,  but  this  is  again  a  case  of  a  notch  radius 
less  than  0.01  inch,  for  which  the  scatter  was  larger  on  the  rotating 
beam  tests,  also. 

Figures  4,  5^  and  6  may  be  summarized  as  follows:  The  Neuber  for¬ 
mula,  used  in  conjunction  with  the  tentative  curve  of  R'  plotted 
against  tensile  strength,  predicts  the  fatigue  factor  from  the  theo¬ 
retical  factor  to  an  accuracy  of  ±10  percent  for  91  percent  of  the  tests 
if  test  specimens  with  radii  less  than  0.01  inch  are  excluded  and  for 
66  percent  if  no  test  specimens  are  excluded.  The  tests  Include  a  large 
variety  of  steels,  with  tensile  strengths  varying  from  60  to  225  ksi, 
very  different  configurations,  and  bending  tests  as  well  as  axial  loading 
tests. 

An  extreme  form  of  a  notch  is  a  crack.  For  this  discussion,  a 
crack  may  be  defined  as  a  notch  of  known  depth  and  very  small  but 
unknown  width.  From  this  definition,  it  follows  that  a  stress- 
concentration  factor  cannot  be  calculated  by  the  theory  of  elasticity, 
because  this  theory  requires  that  the  geometry  of  the  notch  be  com¬ 
pletely  defined.  With  the  Neuber  concept,  on  the  other  hand,  the  crack 
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is  the  limiting  case  of  a  notch,  and  the  factor  can  he  calculated 
readily  from  the  depth  of  the  crack  alone. 

Five  sets  of  tests  on  specimens  with  artificial  cracks  were  avail¬ 
able.  R,  E.  Peterson  produced  cracks  of  two  different  depths  by 
machining  a  narrow  groove  and  then  compressing  the  specimens  until  the 
gap  vas  closed  (reference  13).  In  German  tests  (reference  l6),  speci¬ 
mens  made  of  three  nitrlded  steels  were  stretched  ixntil  they  cracked. 

As  might  be  expected  for  this  type  of  specimens,  the  agreement  between 
calculated  and  experimental  fatigue  factors  is  not  so  good  as  for  the 
other  tests  shown.  The  discrepancies  were  of  the  order  of  ±25  percent, 
with  a  maximvim  discrepancy  of  50  percent.  For  the  German  tests,  a 
photomicrograph  made  it  possible  to  estimate  the  width  of  the  crack; 
with  this  information,  the  theory  of  elasticity  can  be  used  to  make  an 
order-of-magnitude  estimate  of  the  theoretical  factor.  The  error  was 
over  1000  percent;  the  Neuber  factor  was  therefore  definitely  a  better 
approxlmat ion. 

The  Neuber  block  concept  is  a  mathematical  device  for  expressing 
the  physical  idea  that  nature  Imposes  limitations  on  stress  gradients, 
at  least  in  materials  with  a  granular  structure.  Now,  stress  gradients 
exist  not  only  aroiind  notches;  a  beam  also  has  a  gradient  across  its 
depth.  This  suggests  that  the  Neuber  concept  might  be  applied  in  some 
way  to  estimate  the  size  effect  on  smooth  rotating  beams. 

It  has  long  teen  known  that  the  fatigue  strength  decreases  as  the 
diameter  of  the  rotating  beam  increases.  H.  F.  Moore,  for  instance, 
ran  a  series  of  tests  on  six  steels,  with  specimen  diameter  varied  from 
l/8  inch  to  about  2  inches,  and  found  that  the  smallest  specimens  were 
up  to  30  percent  stronger  than  the  leorgest  ones  (reference  8).  The 
lower  limit  of  strength,  which  might  be  called  the  true  strength  of 
the  material,  in  a  way,  would  be  shown  by  specimens  with  infinite 
diameter  or  zero  gradient.  It  is  not  possible  to  test  a  rotating  beam 
of  infinite  diameter,  but  it  is  relatively  easy  to  do  the  equivalent  - 
to  test  specimens  under  axial  loading.  This  was  done  in  some  recent 
German  tests  (reference  17)  on  two  kinds  of  steel;  the  rotating  beam 
specimens  in  these  tests  were  very  small  in  order  to  give  large  gra¬ 
dients.  They  varied  from  about  0.3  inch  down  to  0,0UO  inch  in  diameter. 
The  smallest  beam  specimens  had  36  percent  and  65  percent  more  strength 
than  the  axially  loaded  specimens. 

The  theoretical  method  of  applying  the  Neuber  concept  to  tests  of 
this  type  is  not  clear  at  the  moment.  Qnplrlcally,  it  was  found  that 
the  strength  of  the  rotating  beams  could  be  estimated  quite  well  by 
multiplying  the  axial  strength  by  the  factor 
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where  R  is  now  the  radius  of  the  specimen.  The  strengths  predicted 
in  this  manner  agreed  with  the  observed  ones  with  a  maximum  error  of 
about  6  or  7  percent  for  all  the  tests  mentioned  except  one  of  Moore's 
tests,  which  showed  a  20- percent  discrepancy.  (For  the  tests  made  by 
Moore,  the  strength  under  axial  loading  was  estimated  by  applying  the 
formula  to  the  rotatlng-beam  strengths  obtained  for  specimens  of  varying 
diameters.  In  effect,  then,  the  formula  checked  the  effect  of  varying 
the  diameter  in  the  test  range,  but  not  the  relation  of  rotating-beam 
strength  to  axial  strength. ) 

For  steel,  then,  there  seems  to  be  good  hope  that  the  geometrical 
size  effect  may  be  fairly  well  predictable.  Additional  work  remains  to 
be  done,  of  covirse,  to  clear  up  the  existing  discrepancies  if  possible, 
to  extend  the  method  to  other  cases,  and  to  define  its  limitations. 

On  aluminum  alloys,  little  can  be  said  at  present.  There  are  very 
few  systematic  investigations  of  size  effect,  and  the  results  from  other 
tests  are  very  conflicting.  A  first  approximation  using  R'  equal  to 
0.02  inch  as  for  steel  gives  good  agreement  for  sets  of  tests  of  some 
alloys,  but  even  this  first  approximation  does  not  look  as  promising 
generally  as  it  did  for  steel,  and  attempts  at  a  second  approximation 
have  been  xinsuccessful  so  far.  It  is  suspected  that  the  trouble  is 
caused  by  machining  effects.  This  suspicion,  voiced  by  various  inves¬ 
tigators,  is  supported  by  British  X-ray  measurements  (reference  l8)  of 
siurface  stresses  caused  by  machining  of  cast  aluminum  alloys,  which 
sxiggest  that  the  form  tools  used  to  produce  V-notches  may  produce  very 
high  resldvial  stresses. 
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Figure  1.-  Relation  betveen  Reuber  constant  R*  and  tensile 
strength  for  steel  specimens. 
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Figure  2.-  Coiiqparlaone  between  calculated  and  experimental 
fatigue  factors  for  rotating  beams  with  grooves. 
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Figure  Conqparlson  between  calculated  and  experimental 
fatigue  factors  for  filleted  shafts  (rotating  beam  tests) 


Figure  4.-  Correlation  between  calculated  and  eiqmrlmental 
fatigue  factors  for  rotating  beams  with  grooves  and 
V-notches. 
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Figure  Correlation  between  calculated  and  experimental 
fatigue  factors  for  filleted  shafts. 
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FlguK  6.-  Correlation  between  calculated  and  eaqperlmental 
fatigue  factors  for  axially  loaded  specimens. 
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INVESTIGATION  OF  FATIGUE  STRENGTH  OF  A 
FULL-SCALE  AIRPLANE 
Dwight  0.  Fearnow 
Langley  Aeronautical  Laboratoxy 


Recent  trends  in  airplane  design^  construction,  and  use  have  placed 
Increasing  emphasis  upon  the  problem  of  fatigue  in  airplane  structures. 
Unfortunately,  an  exact  theoretical  method  of  computing  fatigue  life  is 
not  available.  Moreover,  most  fatigue  tests  to  date  have  been  conducted 
on  small  polished  specimens  ;inder  closely  controlled  laboratory  condi¬ 
tions,  but,  although  such  basic  tests  are  essential  to  the  solution  of 
the  problem,  they  do  not  reproduce  all  the  practical  problems  of  a 
C0II^}lex  wing  structure. 

VAien  it  became  apparent  that  war  surplus  airplanes  could  be  made 
available  for  research  on  full-scale  test  specimens,  the  NACA  Special 
Committee  on  Surplus  Aircraft  Research  endorsed  the  proposal  that  an 
investigation  of  static  and  fatigue  strength  of  new  airplane  wing  struc¬ 
tures,  as  well  as  those  exposed  to  long  war  service,  be  carried  out 
through  the  testing  of  sxirplus  aircraft.  The  scope  of  the  program  was 
originally  quite  large.  It  covered  the  testing  of  some  600  airplanes 
of  15  different  types  and  was  to  be  carried  out  through  the  coordinated 
efforts  of  the  Air  Forces,  Navy,  aircraft  industiy,  and  NACA.  It  soon 
became  apparent,  however,  that  funds  and  personnel  would  not  be  avail¬ 
able  for  such  an  ambitious  program.  The  steering  committee  then  recom¬ 
mended  that  the  NACA  carzy  out  the  initial  phase  of  the  program  on  one 
type  of  transport  airplane.  The  NACA  proceeded  to  secure  23  of  the  most 
suitable  planes  available  from  the  War  Assets  Administration. 

The  original  objectives  of  the  program  had  to  be  modified  to  be 
consistent  with  the  drastic  reduction  in  the  number  of  airplanes  to  be 
tested.  Perhaps  the  most  inqjortant  objective  of  the  modified  program 
would  be  the  determination  of  some  correlation  between  the  results  of 
laboratoxy  tests  on  small  specimens  and  their  application  to  the  full- 
scale  wing  structure.  Other  objectives  would  be  the  determination  of 
the  spread  in  fatigue  life  in  wing  structures  which  are  as  nearly  iden¬ 
tical  as  assembly- line  production  methods  will  allow,  the  reduction  in 
static  strength  after  fatigue  failure,  the  effect  of  fatigue  damage  on 
such  parameters  as  natural  frequency  and  damping  characteristics,  and, 
if  possible,  the  loss  in  fatigue  life  associated  with  the  number  of 
hours  a  transport  airplane  has  flown. 

The  results  to  be  presented  here  are  the  first  available  from  this 
investigation  and  deal  entirely  with  tests  of  the  constant-level  type. 
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Subsequent  tests  should  yield  infomation  from  additional  constant-level 
tests  as  well  as  tests  in  which  the  stress  level  is  varied  in  a  random 
step  fashion. 

The  present  Investigation  is  being  conducted  on  the  wing  structure 
of  the  CuHiss-Wrlght  C-U6D  transport  airplane.  The  two  test  specimens 
troa.  which  data  have  been  obtained  were  subjected  to  about  600  hovirs  of 
flight  and  several  years  storage  in  an  open  depot.  A  test  specimen 
consists  of  a  left-  and  a  right-wing  outer  panel  and  a  center  section. 

The  outer  panels  are  of  the  two-spar  type  of  construction  while  the 
center  section  has  three  spars.  The  spars  are  made  of  2liS-T  aluminum 
alloy  but  all  skin  and  stiffeners  are  made  from  Alclad  2liS-T. 

Certain  modifications  are  made  to  each  airplane  to  prepare  it  for 
testing.  The  fuselage  is  cut  in  front  of  and  behind  the  wing,  and  this 
section  is  then  lnvez*ted  and  supported  between  structural  steel  back¬ 
stops  .  The  upper  portion  of  figure  1  shows  a  profile  view  of  half  the 
test  specimen  in  its  Inverted  position.  The  complete  fuselage  and  back¬ 
stops  are  omitted  for  clarity.  The  test  specimen  is  symmetrical  about 
the  center  line.  The  outer  wing  panels  are  cut  off  at  a  station 
U05  Inches  from  the  center  line  of  the  wing-fuselage  combination.  This 
reduced  the  total  span  by  about  35  percent.  Concentrated  masses  are 
attached  to  the  wing  at  the  point  where  it  was  cut  off  to  reproduce 
level-flight  stresses  at  station  21]^.  Station  21ii  was  designated  as 
critical  for  the  fatigue  tests  because  previous  tests,  utilizing  brittle 
lacquer,  indicated  numerous  points  of  local  high-stress  concentrations 
at  that  station.  Design  margins  of  safety  were  also  very  low  in  that 
region  of  the  wing.  (See  reference  1.)  Actually,  design  bending  moments 
were  reproduced  fairly  closely  over  a  considerable  portion  of  the  wing 
as  can  also  be  seen  in  the  lower  portion  of  flgvire  1,  where  bending 
moment  is  plotted  against  wing  semispan  station  in  inches  from  the  center 
line.  The  upper  dashed  curve  shows  that,  if  relieving  loads,  resulting 
from  nacelle  Inertia  effects,  were  not  introduced  during  the  laboratory 
tests,  design  bending  moments  (the  solid  line,  as  obtained  from  refer¬ 
ence  2)  would  be  reproduced  within  about  1  percent  over  a  considerable 
portion  of  the  test  specimen.  Since  it  appeared  desirable  to  force  the 
wing  to  fall,  if  possible,  at  station  21k,  nacelle  Inertia  effects  were 
introduced  by  means  of  a  spring  system  at  the  engine  mount,  and,  as 
shown  by  the  lower  dashed  line,  the  bending  moments  were  reduced  by 
about  10  percent  at  the  center  line. 

All  tests  were  conducted  by  the  resonant-frequency  method.  The 
.concentrated  masses  attached  to  the  test  specimen  provided  the  mean 
stress  condition  about  which  the  wing  vibrated.  The  fatigue  machine  now 
being  used  for  the  constant-level  tests  is  relatively  simple.  Basically, 
it  is  an  adjustable  eccentric  which  tranraits  vertical  force  through  a 
push  rod  and  a  spring  to  the  wing.  The  spring  is  necessary  to  permit 
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the  reqvdred  phase  shift  between  wing  and  forcing  fimction  so  that  the 
wing  resonance  characteristics  can  be  utilized. 

Figure  2  shows  a  plan  view  of  a  portion  of  the  right  semispan  of 
the  test  specimen.  In  this  figiire  the  fatigue  cracks  are  denoted 
arrows  drawn  to  the  points  vAiere  the  cracks  started.  The  cracks  are 
omitted  for  reasons  of  clarity.  The  straight  and  wavy  dashed  lines 
indicate  points  where  fatigue-crack-detector  wires  were  installed,  and 
the  squares  Indicate  the  nominal  locations  of  the  strain  gages. 

Although  this  figtire  shows  only  the  right  semispan,  the  points  of 
origin  of  all  eleven  fatigue  cracks,  from  both  the  left  and  right  semi¬ 
spans  of  the  first  two  specimens,  are  included.  For  purposes  of  identi¬ 
fication,  the  cracks  are  numbered  consecutively  in  or^r  of  occurrence. 

The  first  five  cracks  occurred  in  the  first  test  specimen,  whereas  the 
last  six  occurred  in  the  second  test  specimen.  The  R  and  L  beside 
each  crack  number  indicates  the  semispan,  right  or  left,  on  which  the 
crack  occurred.  In  three  Instances,  Involving  seven  of  the  eleven 
cracks,  duplicate  failures  were  obtained.  Cracks  2R,  IL,  and  IIL  occurred 
in  riveted  joints  which  were  symmetrical  about  the  center  line.  Cracks  6R 
and  7L  originated  at  identical  points  on  the  right  and  left  semi spans  at 
the  outboard  jiincture  of  the  wing  and  nacelle  group,  while  cracks  8R 
and  9L  originated  at  the  corners  of  Identical  Inspection  cut-outs. 

Cracks  3R  and  5L  occurred  at  the  comers  of  inspection  cut-outs  in  the 
outer  panels,  whereas  cracks  I4L  and  lOL  originated  in  the  joggle  in 
doubler  plates  located  on  the  outboard  and  inboard  sides  of  the  wing 
attachment  angle. 

The  strain  gages,  shown  by  the  squares,  were  located  so  as  to 
measure  only  nominal  stress  in  the  material  in  the  vicinity  of  the  points 
where  fatigue  cracks  originated  and  were  far  enough  away  from  the  point 
of  fatigue-crack  origin,  so  as  to  measure  no  stress  concentrations.  The 
fatigue-crack-detector  wires,  the  straight  and  wavy  dashed  lines,  were 
located  at  points  where  the  brittle-lacquer  tests  indicated  spots  of 
local  high  stresses. 

Although  it  was  intended  to  vibrate  the  test  specimens  about  the  1  g 
or  level-flight  condition  at  a  constant  ai^plitude  equivalent  to  an  incre¬ 
ment  of  load  factor  An  of  slightly  more  than  A0.6g,  this  could  be 
achieved  only  to  a  limited  degree.  In  the  starting  and  stopping  of  the 
machine,  a  few  cycles  of  load  were  applied  at  higher  values  of  An  as 
well  as  at  lower  values,  but  calculations  have  indicated  that  the  effect 
of  these  loads  was  negligible. 

Since  the  tests  were  not  exact  constant-level  tests  according  to 
standards  set  forth  in  fatigue  tests  of  small  specimens,  for  purposes 
of  analysis  the  assuiEq>tlon  was  made  that  only  those  cycles  of  load 
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applied  in  the  Interval  betvfeen  An  •  0.6g  and  An  •  0.75g  had  any 
effect  upon  the  fatigue  life  of  the  specimen.  It  was  further  asstimed 
that  those  loads  acted  at  the  midpoint  of  that  interval.  The  effec¬ 
tive  stress-concentration  factor  was  then  coiqmted  hy  simply 

dividing  the  stress  from  the  S-N  curve  of  imnotched  Alclad  2US-T 
aluminum-alloy  sheet,  as  found  in  reference  3,  hy  the  measured  stress. 
The  measured  stress,  as  mentioned  previously,  included  none  of  the 
stress  concentrations.  The  values  obtained  for  are  shown  in  the 
following  table: 


For  cracks  1,  2,  and  11,  all  of  which  originated  in  similar  riveted 
joints,  varies  from  2,6  to  3.2.  The  comer  of  an  inspection  cut¬ 
out,  where  crack  3  originated,  had  an  effective  stress-concentration 
factor  of  U>2.  The  table  also  shows  the  number  of  cycles  of  load 
required  to  cause  failure,  as  well  as  the  maximum  effective  stress. 
Maximum  effective  stress  is  defined  as  measured  mean  stress  plus  half 
of  the  alternating  stress  all  multiplied  hy  the  effective  stress- 
concentration  factor.  Although  it  may  be  mere  coincidence,  it  is  inter¬ 
esting  to  note  the  relatively  small  amount  of  spread  in  the  niunber  of 
cycles  required  to  cause  failure  for  each  of  the  eleven  cracks,  even 
though  all  cracks  did  not  originate  in  the  same  places  on  the  two  test 
specimens.  The  maximum  spread  is  in  the  order  of  a  factor  of  2;  the 
lowest  number  is  170,000  and  the  maximum,  308,000.  This  could  lead  to 
seme  speculation  on  the  possibility  of  determining  airplane  S-N  curves, 
and  a  word  of  caution  should  he  Injected  here  concerning  such  specu¬ 
lation.  The  relatively  small  amount  of  data  available  to  date  shovdd 
not  be  overlooked  and  subsequent  data  may  show  considerably  more  scatter. 

Figtire  3f  in  which  loss  of  tension  material  in  percent  is  plotted 
against  (^las,  shows  grai^cally  how  cracks  3  and  U  grew  after 
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Inception.  The  abscissa  is  shown  on  rectangular  coordinates  rather  than 
a  logarithmic  scale  to  enqphasize  the  fact  that  quite  a  number  of  cycles 
are  required  to  cause  appreciable  loss  of  tension  material.  It  is 
interesting  to  note  in  this  figure  that  the  rate  of  crack  growth  was 
quite  smsdl  until  between  5  and  10  percent  of  the  tension  material 
failed.  From  that  point  the  rate  of  growth,  as  shown  by  the  steepness 
of  the  curves,  was  quite  rapid.  It  was  Interesting  to  watch  the  growth 
of  the  crack  during  the  test.  In  each  case  the  crack  in  the  skin  would 
pass  over  a  stringer  and  proceed  to  about  halfway  to  the  next  stringer 
and  then  the  first  stringer  would  fail,  and  so  on  across  the  wing  in 
almost  a  straight  line. 

Figure  li  shows  the  reduction  in  natural  frequency  of  the  wing  as 
a  function  of  loss  of  tension  material  for  cracks  3  and  U.  The  natiiral 
frequency  of  the  test  specimen,  before  the  test  started,  was  about 
106  cycles  per  minute.  After  about  55  percent  of  the  tension  material 
had  failed,  this  value  had  only  been  reduced  by  about  2  percent  or  to 
about  103.5  cycles  per  minute.  No  change  in  natural  frequency  was 
apparent  until  after  the  fatigue  cracks  had  originated. 

Although  no  specific  damping  measurements  were  taken  during  the 
test,  the  very  nature  of  the  test  procediire,  that  is,  the  resonant- 
frequency  method,  was  such  that  any  change  in  damping  should  have 
resulted  in  a  change  in  amplitude  of  vibration.  No  such  change  in 
amplitude  was  apparent  until  after  a  fatigue  crack  had  originated. 


In  conclusion,  a  fatigue-strength  investigation  is  currently  being 
conducted  on  full-scale  airplane-wing  structures  by  the  resonant- 
frequency  method,  wherein  concentrated  masses  are  attached  to  the  wing 
to  reproduce  level-flight  stresses  over  approximately  h5  percent  of  the 
span.  The  res^llts  presented  Include  the  data  gathered  thus  far  on  the 
first  two  fatigue  specimens.  Although  the  results  cover  eleven  fatigue 
cracks,  there  are  not  yet  sxifficlent  data  on  which  to  base  definite  con¬ 
clusions.  The*  data  do  indicate  a  surprisingly  small  amount  of  spread 
in  fatigue  life  even  though  all  fatigue  cracks  did  not  originate  at  the 
same  points.  Computed  effective  stress-concentration  factors  for  similar 
riveted  Joints  varied  from  2.6  to  3*2,  whereas  the  one  value  computed 
for  the  comer  of  an  inspection  cut-out  was  U.2.  Neither  the  natural 
frequency  nor  the  danplng  of  the  two  test  specimens  appeared  to  be 
affected  by  fatigue  damage  until  after  a  fatigue  crack  had  originated. 

The  rate  of  fatigue-crack  growth  was  quite  small  until  approximately 


percent  of  tension  material  was  lost. 
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Figure  1.-  Comparison  of  design  and  test  applied  bending  moments 

over  span  of  test  specimen. 


Figttn  2.-  Moalnal  location  of  inatrunentatlon  and  points  of 

fatigue-crack  origin. 
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Figure  3«-  Typical  curves  of  fatigue-crack  growth 
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Figure  4.-  Effect  of  fatigue-crack  damage  on  natural  frequency 

of  test  specimens. 
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